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Preface 


This volume contains the results of the Workshop on Unmixing the 
SNCs: Chemical, Isotopic, and Petrologic Components of Martian Meteor- 
ites, held October 1 1-12, 2002, in Houston, Texas. Scientific organizers for 
the workshop were Allan Treiman and Christopher Herd (Lunar and Plan- 
etary Institute), assisted by co-conveners John Jones and David Mittlefehldt 
(NASA Johnson Space Center). 

Logistics, administrative, and publications support were provided by the 
Publications and Program Services Department of the Lunar and Planetary 
Institute. 


Contents 


Program 

Workshop Summary 
Abstracts — 

Incompatible Trace Elements in SNCs and Li Isotopes Systematics in a Shergottite 

J ’ A . Barrat, P> Beck, M. Chaussidon, Ph. GiUet, A. Jambon, C Gopel, and V. Sautter 

Exploring Trace Element and Isotopic Mixing Relationships in the Martian Meteorite Suite 

L. E. Borg 

What Should We Looking for in Martian Meteorites? Is Evidence of Crustal Process or Mantle Process 
More Important — and to Whom? 

J. S. Delaney and M. D. Dyar 

Crust-Mantle Reservoirs of Radiogenic Isotopes of Mars and Earth: Where Can We See a Mixing? 

G. Dreibus and E. Jagoutz 

Mossbauer Spectroscopy of SNC Meteorites 

M. D . Dyar 

Petrogenesis of Olivine-Phyric Shergottites Sayh A1 Uhaymir 005 and Elephant Moraine A79001 
Lithology A 

C. A. Goodrich 

Mars Meteorite Statistics and the Martian Uplands 
W. K. Hartmann 

Martian Basalt Oxygen Fugacity and Geochemistry: Implications for a Heterogeneous Martian Mantle 

C. D. K. Herd 

Origin of the Martian Crust and Mantle 
P. C. Hess 

New Aspects in the Isotope Systematics of Shergottites 
E. J. Jagoutz and G . D. Dreibus 

SNC Meteorites and Martian Reservoirs 
J. H. Jones 

Melting the Martian Mantle: Shergottite Formation and Implications for Present Day Mantle Convection 
on Mars 

W. S . Kiefer 

QUE 94201 : Reconsidering Its Origins as a Bulk Melt from a Volcanic Region of Mars 

D. A. Kring 

SNC Meteorites and Their Source Composition(s) 

J. Longhi 


Oxygen Fugacity Recorded in Pigeonite: Indications of a Heterogeneous Martian Magma Source Region? 
M. C. McCanta and M. J. Rutherford 

Geochemistry of Martian Meteorites and the Petrologic Evolution of Mars 

D. W. Mittlefehldt 

Silicon in Mars’ Core: A Prediction Based on Mars Model Using Nitrogen and Oxygen Isotopes in 
SNC Meteorites 

R . K. Mohapatra and S. V. S. Murty 

REE in Shergottite Augites and Whole Rocks 
£>. S. Musselwhite and M. Wadhwa 

An “Andesitic” Component in Shergottites with Restored LREE Abundances? 

L. E. Nyquist, C.-F. Shift, H . Wiesmann, and J. A. Barr at 

Evolution of Martian Atmospheric, Crustal and Mantle Xenon Components in Basaltic Shergottites 
K. D. Ocker, G. Holland, and J. D. Gilmour 

Noble Gases in the Terrestrial Planets, with Focus on what the SNC Meteorites Tell Us About Mars 

R. O. Pepin 

Shergottite Impact Melt Glasses Contain Soil from Martian Uplands 

M. K Rao and D. S. McKay 

Is Olivine Speedometer a Reliable Tool to Constrain Thermal Story of SNC Meteorites? 

V. Sautter, J. A . Barrat, and A. Jambon 

Nitrogen and Noble Gases in Mineral Separates from Zagami 

S. P. Schwenzer, R. K Mohapatra, S ' Herrmann, and U. Ott 

Chemical Composition of Newly Collected Antarctic Nakhlites, Y000593 and Y000749 

N. Shirai, Y. Oura, and M. Ebihara 

Some Puzzles About What Noble Gas Components Were Mixed into the Nakhlites, and How 

T. D. Swindle 

Element-Element Correlations Among Martian Meteorite Bulk Compositions: Peculiarities Explained(?) 
by Mixing, with Implications for the Composition of Mars 
A. H. Treiman 

Chemistry and Micro textures of Melt Pockets in the Los Angeles Basaltic Shergottite 

E. L. Walton and J. G. Spray 


List of Workshop Participants 


Program 


Friday, October 11, 2002 

RADIOGENIC ISOTOPES AND OXYGEN FUGACITY 
8:30 a.m. Berkner Rooms 


WELCOME AND ANNOUNCEMENTS 

Borg L. E.* [INVITED 20-MINUTE PRESENTATION] 

Exploring Trace Element and Isotopic Mixing Relationships in the Martian Meteorite Suite 

Lee C. T.* [INVITED 20-MINUTE PRESENTATION] 

A Brief Overview of Mantle Metasomatism 

Ripley E.* [INVITED 20-MINUTE PRESENTATION] 

Chemical Effects on Interaction Between Magma and Wall Rock 

Dreibus G.* Jagoutz E. 

Crust-Mantle Reservoirs of Radiogenic Isotopes of Mars and Earth: Where Can We See a Mixing? 

Nyquist L. E.* Shih C.-Y. Wiesmann H. Barrat J. A. 

An “Andesitic” Component in Shergottites with Restored LREE Abundances? 

Herd C. D. K * 

Martian Basalt Oxygen Fugacity and Geochemistry : Implications for a Heterogeneous Martian Mantle 

Musselwhite D. S.* Wadhwa M. 

REE in Shergottite Augites and Whole Rocks 

McCanta M. C* Rutherford M. J. 

Oxygen Fugacity Recorded in Pigeonite: Indications of a Heterogeneous Martian Magma Source Region? 


Friday, October 11, 2002 

BASALTS AND PHASE EQUILIBRIA (ON MARS) 
1 :30 p.m. Berkner Rooms 


Glazner A.* [INVITED 20-MINUTE PRESENTATION] 

Thermal Constraints on Magma Mixing and Assimilation 

Hess P.* [INVITED 20-MINUTE PRESENTATION] 

Origins of the Martian Crust and Mantle 

Longhi J.* 

SNC Meteorites and Their Source Composition(s) 

Kring D. A.* 

QUE 94201 : Reconsidering Its Origins as a Bulk Melt from a Volcanic Region of Mars 
Goodrich C. A.* 

Petrogenesis of Olivine-Phyric Shergottites Sayh A l Uhaymir 005 and Elephant Moraine A 79001 Lithology A 


* Denotes Speaker 




Friday, October 11, 2002 
POSTER SESSION 
5:00 - 7:00 p.m. Great Room 


Jagoutz E. J. Dreibus G. D. 

New Aspects in the Isotope Systematics of Shergottites 

Schwenzer S. P. Mohapatra R. K, Herrmann S. Ott U. 

Nitrogen and Noble Gases in Mineral Separates from Zagami 

Dyar M. D. 

Mossbauer Spectroscopy of SNC Meteorites 
Hartmann W. K. 

Mars Meteorite Statistics and the Martian Uplands 
Shirai N. Oura Y. Ebihara M. 

Chemical Composition of Newly Collected Antarctic Nakhlites, Y000593 and Y000749 

Sautter V. Barrat J. A. Jambon A. — POSTER CANCELED 

Is Olivine Speedometer a Reliable Tool to Constrain Thermal Story of SNC Meteorites ? 

Ocker K. D. Holland G. Gilmour J. D. 

Evolution of Martian Atmospheric , Crustal and Mantle Xenon Components in Basaltic Shergottites 
RaoM.N. McKay D. S. 

Shergottite Impact Melt Glasses Contain Soil from Martian Uplands 
Mohapatra R. K. Murty S. V. S. 

Silicon in Mars * Core: A Prediction Based on Mars Model Using Nitrogen and Oxygen Isotopes in SNC Meteorites 
Mittlefehldt D. W. 

Geochemistry of Martian Meteorites and the Petrologic Evolution of Mars 
Walton E. L. Spray J. G. 

Chemistry and Microtextures of Melt Pockets in the Los Angeles Basaltic Shergottite 



Saturday, October 12, 2002 
GEOPHYSICS AND TRACE ELEMENTS 
8:30 a.m. Berkner Rooms 


Parmentier M.* [ INVITED 20-MINUTE PRESENTATION] 

Convection and Layering in the Martian Mantle 

Kiefer W. S* 

Melting the Martian Mantle: Shergottite Formation and Implications for Present Day Mantle Convection on Mars 
Jones J. H.* 

SNC Meteorites and Martian Reservoirs 
Treiman A. H.* 

Element-Element Correlations Among Martian Meteorite Bulk Compositions: Peculiarities Explained(?) by Mixing, 
with Implications for the Composition of Mars 

BarratJ. A.* BeckP. Chaussidon M. GilletPh. JambonA. GopelC. Sautter V. — TALK CANCELED 
Incompatible Trace Elements in SNCs and Li Isotopes Systematics in a Shergottite 

Saturday, October 12, 2002 
NOBLE GASES AND SUMMATION 
1:30 p.m. Berkner Rooms 


Pepin R. O * [INVITED 20-MINUTE PRESENTATION] 

Noble Gases in the Terrestrial Planets, with Focus on what the SNC Meteorites Tell Us About Mars 

Swindle T. D.* 

Some Puzzles About What Noble Gas Components Were Mixed into the Nakhlites, and How 
Delaney J. S.* Dyar M. D. 

What Should We Looking for in Martian Meteorites? Is Evidence of Crustal Process or Mantle Process 
More Important — and to Whom? 


GENERAL DISCUSSION 


WORKSHOP ADJOURNS 


WORKSHOP SUMMARY 
Unmixing the SNCs 


1. INTRODUCTION 

Geochemical and petrologic studies of the Martian meteorites (nicknamed the “SNCs”) 
have proliferated in the past few years, from a wealth of new samples and the perfection of new 
analytical methods. An intriguing result from these studies is that the chemical and isotopic 
compositions of the Martian meteorites, all basalts or derived from basaltic magma, can be 
modeled as mixtures of a limited number of components. 

These mixing components were the focus of the workshop, and the attendees brought 
their expertise to explore several fundamental questions about the components. 

• How many mixing components can be recognized? 

• Can components defined in one type of system be correlated with components in another 
system? 

• Do the components represent physical masses of material, or chemical isotopic 
processes? 

• How and when did the components form, i.e., become separated from the average bulk 
composition of Mars? 

• How and when did the components come to be mixed into the SNCs? 

• What are the historical, tectonic, and geological implications of the components? 

Few definitive answers emerged, and finding them was not the workshop’s immediate 
purpose. Rather, the workshop served as a forum for active researchers to present and compare 
their results, and perhaps as a source of inspiration. A continuing theme was the comparison of 
Martian geochemistry with those of other planets, especially the Moon. The issues were set 
against the commonly agreed results that Mars differentiated early, -4.5 Ga, into a core and 
several silicate materials that could include mantle reservoirs, crust, and/or residua of a magma 
ocean. Several speakers discussed the idea of an early magma ocean, as many chemical systems 
in the Martian meteorites are closely analogous to those seen in lunar basalts. 

2. WORKSHOP ATTENDANCE, ORGANIZATION, AND FORMAT 

“Unmixing the SNCs” was held on October 1 1-12, 2002, at the Lunar and Planetary 
Institute in Houston, Texas. Fifty scientists and students attended, representing research groups 
from the U.S.A., Germany, Canada, Japan, and Hawai‘i. The scientific sessions included 23 oral 





presentations, seven of which were invited, and eight posters. The oral presentations were 
divided into four sessions, focused generally on radiogenic isotopes, petrologic and chemical 
systematics, geophysical constraints and trace element distributions, and noble gases. Ample 
time was provided for questions and answers, which ranged widely and engaged nearly all of the 
audience. 

3. OVERVIEW OF SCIENTIFIC PROGRAM 

The abstracts associated with the workshop were published in a hard-copy abstract 
volume, limited supplies of which are available upon request from the LPI 1 . The preliminary 
program and electronic versions of the abstracts are available via the LPI’s Web site at 
www.lpi.usra.edu/meetings/sncs2002/, and the abstracts are also included in this volume. In 
some cases, presentations differed considerably from the abstracts, which represents the rapid 
pace of investigations (and reinterpretation) of the Martian meteorites. 

3. 1. Terrestrial Analogs and Thermochemistry 

Several speakers from outside the Martian meteorite community were invited to the 
workshop to provide “real-world” analogs and constraints on the mixing processes. 

Dr. Edward Ripley (Indiana University, Bloomington) spoke about chemical interactions 
between basaltic magmas and surrounding rocks, illustrated by his work on massive sulfide ore 
bodies . Sulfur isotope ratios in these ores are characteristic of the surrounding metasediments 
(i.e., are characteristic of biogenic sulfate reduction) and not of mantle sulfur. Other isotopic and 
chemical systems are also consistent with extensive element mobility into the basalt magma, 
presumably mediated by aqueous fluids. Elsewhere in these deposits, fragments of the wall-rock 
metasediments were caught up in the magma. These fragments were partially melted by the 
basalts, and the melts (rich in incompatible elements) were removed from the fragments into the 
basalt magma. Traces of this assimilation are seen in trace element and isotopic signatures in the 
ore deposits and the basalt flows related to them. 

Dr. C.-T. Lee (Rice University) spoke about the many complications of mantle 
metasomatism on Earth, both as seen in mantle samples and in basalts melted from 
metasomatized mantle 3 . Mantle material in the Earth commonly has much greater abundances of 
highly incompatible trace elements (e.g., La, U, Th) than would be expected from the Mg/Fe 
ratios of their olivine and pyroxenes. This enrichment, the metasomatism, is caused by passage 
of low- volume partial melts from other parts of the mantle. These melts can be silicate-rich, 
carbonate-rich, and/or water-rich depending on their source region and its history. Thus, 
metasomatized mantle can have a wide range of compositions as can basalts melted from 
metasomatized mantle. 

Dr. A. Glazner (University of North Carolina) spoke about the thermal constraints on 
magma mixing 4 . Many interpretations of mixing and assimilation in the literature are based on 



incorrect interpretations from temperature-composition phase diagrams. Phase diagrams showing 
enthalpy (heat content) versus composition are much more appropriate for mixing and 
assimilation, as they permit one to derive the correct temperature, phases present, and relative 
proportions of phases after mixing. 

3.2. Isotopic Constraints 

Radioisotope systematics provide the strongest evidence for component mixing in the 
SNC meteorites. L. Borg led off the workshop with radioisotope data that strongly suggested 
mixing of isotopic components in the source regions of the Martian meteorites. These 
systematics are comparable to those of the Moon, i.e. cumulates from the lunar magma ocean, 
and KREEP (the late-stage silicate magma from the magma ocean). Borg introduced an e 143 Nd - 
s 142 Nd isochron diagram (the former derived from !47 Sm, half life 106 Ga; the latter from 146 Sm, 
half life 103Ma), on which most of the shergottite meteorites (“S” of SNC) are consistent with a 
single differentiation event at 4.51 Ga (and nothing between then and eruption); this is in 
agreement with Rb-Sr and Pb-Pb whole-rock isotopic compositions. G. Dreibus reviewed earlier 
isotopic and chemical evidence for the SNC meteorites, and related the meteorites to the Mars 
Pathfinder rock composition via their potassium contents. L. Nyquist introduced radioisotope 
data for the newly discovered NWA1068 shergottite, and used it and earlier data to show that the 
shergottites could be resolved into a mantle component and a crustal component. The crustal 
component is similar in some respects to the Mars Pathfinder “andesite” composition. J. Jones 
considered many radioisotope systems to show that the Martian meteorites and their source 
mantles are shockingly similar to those of the lunar basalts and their source regions; he favors a 
strongly layered mantle to preserve these chemical components over geological time. Treiman 
considered the bulk compositions of SNC meteorites as a backdrop for radioisotope systematics, 
and found that the effects of silicate fractionation and component mixing could (to some extent) 
be disentangled. At least three geochemical components can be recognized, and SNCs in a single 
age group seem to have similar proportions of the components. In a poster, E. Jagoutz raised the 
minority view that the isotopic characters of the SNCs do not indicate recent crystallization ages. 

3.3. Constraints from Redox 

Constraints on Martian basalt petrogenesis and models of the Martian interior from 
oxygen fugacity estimates are relatively recent 5 . C. Herd presented some new estimates of 
oxygen fugacity from mineral equilibria, and D. Musselwhite and M. McCanta presented 
estimates from the calibration of the Eu/Gd (or Eu/Sm) oxybarometer for Martian augites and 
pigeonite pyroxenes. All three methods show a large spread in oxidation states of the Martian 
magmas, although there are differences in the exact values of the oxygen fugacities. Herd 
reviewed the crustal assimilation models that have resulted from correlations between oxygen 
fugacity estimates and geochemical parameters, and explored in more detail implications of 



oxygen fugacity for a heterogeneous mantle (e.g., as presented by Borg). Oxygen fugacity 
estimates require that the enriched reservoir within the mantle be oxidizing, either due to greater 
amounts of ferric iron, or water (likely bound up in hydrous minerals). Isotopes require that the 
enriched reservoir was separated from the depleted reservoir at 4.5 Ga. Therefore, the 
crystallization of a hydrous magma ocean is a viable model to explain the characteristics of the 
SNCs. 

3.4. Constraints from Petrology 

P. Hess discussed two potential models for the Martian interior using analogies with the 
Earth and the Moon. Continental mantle lithosphere is a potential analog because it is depleted 
mantle yet still garnet-bearing, it has the same age as the continental crust, and there is a 
decoupling of major compatible and trace incompatible elements (due to variable 
metasomatism). The magma ocean model is attractive because it enables heat-producing 
elements to be segregated, and A1 can be sequestered into majorite within the lower part of the 
cumulate pile, thereby accounting for its depletion in the Martian basalts. J. Longhi presented 
many variations on partial melting models, which showed that trace elements are easy to twiddle 
to fit observations, but that major elements are much more difficult. Many variations on partial 
melting models cannot achieve the SNC parental magma compositions from the composition of 
the primitive upper mantle of the Earth. 

D. Kring questioned the evidence for QUE 94201 being a melt (i.e., containing no 
cumulus minerals), since its bulk composition plots in the olivine field. However, during 
discussion G. McKay showed results of experiments done using the QUE 94201 bulk 
composition in which olivine does not crystallize at the oxygen fugacity at which the rock likely 
formed (-IW+1). C. Goodrich presented results of studies of melt inclusions in SaU 005 and 
EET 79001 lithology A showing that the melt inclusions in chromite are more representative of 
parental melts, and that early olivine and chromite, as well as some low-Ca pyroxene are 
xenocrystic. Goodrich also showed that the La/Yb ratios of SaU 005 melt inclusions 
(representing xenolithic magma) and groundmass are similar, whereas the La/Yb ratios of EET 
79001 lithology A xenolith melt inclusions and groundmass are different. These observations are 
matched by differences in oxygen fugacity of the xenoliths and groundmass. Therefore it appears 
that melts of different mantle sources can co-exist in the same rock. 

3.5. Geophysical Constraints 

The workshop provided a unique opportunity to bring geophysicists into the discussion of 
the internal reservoirs of Mars. M. Parmentier presented the results of modeling involving the 
fractional crystallization of the entire silicate portion of Mars. Calculation of densities indicates 
that material at the top of the cumulate pile would be more dense than material at the bottom, and 
that overturn would occur. He addressed the question of the rate of overturn, which depends on 



the size of the denser slab and the viscosity of the mantle. Modeling indicates that overturn could 
occur as quickly as 10 Ma, especially if water was present. A result of overturn is that the colder 
material would collect at the core-mantle boundary, and that heat loss from the core might 
explain an early, brief period of magnetism that is postulated from orbital magnetic 
measurements. 

W. Kiefer, in a laudable interdisciplinary effort, used geophysical modeling of mantle 
convection in Mars to constrain the abundances and locations of Mars, heat-producing elements 
(K, Th, U). A dry Martian mantle can melt enough (and recently enough) to satisfy 
photogeo logic constraints only if >40% of heat-producing elements are in the mantle (not in a 
crust). This seems consistent with the enriched component residing in the mantle. If the Martian 
mantle were wet, the proportion of melting increases significantly, and more of the heat- 
producing elements could have been extracted to a crust. Also tying geophysics and 
geochemistry, S. Jacobsen used the spread of initial isotope ratios for SNC meteorites with a 
model of mantle convection to estimate the average age of the Martian crust as 3. 2-4.0 Ga (vs. 
~2.0 Ga for the Earth), and the timescale of mantle mixing as ~2 Ga (vs. ~0.5 Ga in the Earth). F. 
Singer discussed melting models for the interiors of terrestrial planets, and argued that spindown 
of a satellite captured into retrograde orbit early in a planet’s history could cause significant 
heating through tidal friction. 

3.6. Noble Gas Isotopes 

Noble gas abundances and isotope ratios in SNC meteorite have been studied intensely, 
but it is difficult to relate them to other radioisotope and geochemical systematics. R. Pepin gave 
an overview of processes that fractionate noble gas isotope ratios, emphasizing hydrodynamic 
escape. Pepin focused on Xe, and the difficulties in interpreting abundances of its nine isotopes. 
A serious problem is that Xe (even with nine isotopes) remains ambiguous — within their 
current precision, Xe isotope ratios can be matched by several different starting compositions 
with differing degrees of mass fractionation and addition of flssiogenic components. New Xe 
data muddles the story even further, and it appeared that Pepin’s main argument was that data 
from the SNCs was too equivocal, and that we need to get a direct sample of the Martian 
atmosphere! As outlined by T. Swindle, current theories also do not explain details of Xe isotope 
data in the nakhlites and ALHA84001 . Whereas most of the SNC data falls on a mixing line in 
132 Xe/ 129 Xe vs. Kr/Xe space between Martian atmosphere and Chassigny (thought to be 
representative of Martian interior), other data falls between Chassigny and the Nakhlites/ 
ALH84001, at lower Kr/Xe for a given 132 Xe/ 129 Xe ratio. Available explanations for these effects 
are inadequate, including involving aqueous alteration, adsorption of gases on soil and 
subsequent incorporation into the magma, shock effects, the involvement of clathrates, and even 
the season or latitude at which Nakhla and ALH84001 incorporated the gas! R. Mohapatra 
presented noble gas data on SNC meteorites found in hot deserts, which show an additional 


complication. Several of these SNCs contain a component of elementally fractionated Earth 
air — yet another pitfall in noble gas geochemistry. 

4. OVERALL COMMENTS AND LESSONS LEARNED 

The format of the workshop — short presentations followed by equal (or greater) time for 
discussion — worked well and led to probing, interdisciplinary discussions. These discussions 
would have been difficult at larger, more structured meetings. For example, geophysical models 
of early Mars differentiation provided mechanisms for petrologists and geochemists to explain 
the formation of different internal reservoirs, and to speculate on where these reservoirs might be 
sited. Conversely, geophysicists could constrain their models by the petrologic and geochemical 
observations. The participation of invited speakers from the non-SNC community provided some 
well needed “ground truth”, in the form of descriptions of processes that occur on the Earth that 
may be applicable to Mars. Such insight is important for a community that spends so much of its 
time analyzing meteoritic samples with no geologic context! 

The workshop served to improve the understanding of the various aspects of the SNCs 
within the community, and to raise questions that can guide future research. For example, the 
magma ocean model was invoked repeatedly as an explanation of the internal differentiation of 
Mars. How well does this model explain all of the trace element variations? Do the non- 
lithophile element distributions tell a different story? The nakhlite source appears to be distinct 
from the shergottites. Does this require an entirely separate mantle reservoir from the 
shergottites, a shergottite-like mantle that has been metasomatized, or assimilation of a different 
type of crustal reservoir? 

Finally, J. Delaney raised the issue of relevance to NASA’s Mars exploration program. 
The workshop considered geochemical reservoirs, focusing deeper and deeper into Mars’ mantle 
and into its earliest differentiation. Yet, future Mars missions will have access to, and investigate 
in detail, the Martian surface and near subsurface. How is our work relevant to the Mars 
exploration program? Delaney suggested a need to examine post-crystallization processes that 
affected the SNCs (including shock, metamorphism, alteration, and weathering) in the context of 
future spacecraft missions. If nothing else, workshops like “Unmixing the SNCs” provide a 
baseline for studies of Martian surface processes, be they in situ investigations or analyses of 
returned samples. 
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’CNRS UMR 6112 and Universite d’Angers, 2 bd Lavoisier, 49045 Angers Cedex, barrat@univ-angers.fr, 2 ENS 
Lyon, Dept, de Sciences de la Terre, 46 allee dTtalie, 69364 Lyon Cedex 07, 3 C.R.P.G.-CNRS, 15 rue ND des Pau- 
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In this paper, we report on two aspects of our re- 
cent investigations on Martian meteorites: 1/ Li and 
incompatible trace element abundances in shergottites 
and nakhlites and 2/ the Li isotopes systematics during 
the crystallization of the North West Africa 480 sher- 
gottite (NWA 480). Li is seen by many workers as a 
key element for tracing fluids during the genesis of 
nakhlites and shergottites [e.g. 1]. 

Li and incompatible trace-element abundances: 

In view of their isotopic compositions and their oxida- 
tion states, basaltic shergottites such as Zagami or Los 
Angeles are seen as being contaminated by a Martian 
crustal component [e.g. 2-6). The isotopic composi- 
tions of shergottites are consistent with mixing between 
depleted and enriched components [e.g. 2, 6, 7). Using 
the REE concentrations in pyroxenes, it has been 
shown that the patterns of parental melts of the Martian 
meteorites were parallel to those of the corresponding 
whole rocks [3, 8). Consequently, it can be assumed 
that the REE ratios of these meteorites and their paren- 
tal liquids are similar. This conclusion can be extended 
to other incompatible element ratios such as Th/Sm and 
Rb/Sm. The trends defined by the shergottites can be 
interpreted either as mixing lines between LREE- de- 
pleted and enriched melts or produced by contamina- 
tion of LREE-depleted melts by an enriched compo- 
nent (crust?) [9]. Nakhlites are far outside the trends 
defined by the shergottites, implying a distinct LREE 
component. This interpretation is in agreement with the 
Li abundances determined on whole rock samples. In a 
Li vs. Yb diagram (fig. 1), nakhlites+Chassigny and 
shergottites are well distinguished: shergottites display 
a significant range of Li/Yb ratios (between 1 .2 and 
3.6); Nakhla, NWA 817 and Chassigny display similar 
Li/Yb ratios (about 10). In a Li/Yb vs. Be/Yb plot (Fig. 
2), most of the available data for shergottites define a 
linear trend, indicating one more time depleted and 
enriched components. Two shergottites are clearly out- 
side this trend: DaG 476 and Dhofar 019. Two hy- 
potheses can explain these outliers: (1) these two mete- 
orites are hot desert finds and their (yet) anomalous 
Li/Yb ratios were possibly produced by terrestrial 
weathering; (2) alternatively, other (LREE-depleted) 
components are possibly involved during the petro- 


genesis of these two shergottites, as suggested by their 
unusual high Hf/Lu ratios [9,1 0). . . 

Li isotopes in the NWA 480: The Li abundances 
and isotopic compositions were determined in NWA 
480 pyroxenes by SIMS using a procedure already 
described [11]. This shergottite has been chosen be- 
cause: (1) its pyroxenes are large (typically 2 mm in 
length), and display a huge zoning (from En77 to En5); 
transects are easy to produce; (2) NWA 480 is more 
LREE-depleted than Shergotty and Zagami, which 
have been previously analyzed. Selected results are 
shown in figure 3. 

Li abundances'. About 30 analyses were achieved 
and all the results (but two) are in the range 3-4 ppm. 
Quite different results were obtained on Shergotty and 
Zagami where Li abundances decrease significantly 
from pyroxene cores to rims [1]. 

Li isotopes'. 8 7 Li increases systematically from core 
(about -9 per mil) to rim (+10 - +18 per mil). 

Terrestrial weathering, closed system fractional 
crystallization, mixing or some kind of metasomatism 
failed to explain such huge isotopic variations (30 per 
mil) and Li behavior. In agreement with previous re- 
sults obtained on Shergotty and Zagami pyroxenes [1], 
these results suggest a significant loss of Li (about 60 
%) by the residual magma. The involvement of aque- 
ous fluids is plausible but is not the only possibility. 

References: [1] Lentz R. et al. (2001) GCA 65, 
4551-4565, [2] Jones J. (1989) Proc. L.P.S.C. 19, 465- 
474. [3] Wadhwa M. et al. (1994) GCA 58, 4213-4229. 
[4] Wadhwa M. (2001) Science 291, 1527-1530. [5] 
Herd et al. (2001) Am. Miner. 86, 1015-1024. [6] 
Jagoutz E. (1991) Space Sci. Rev. 56, 13-22. [7] Borg 
L. et al. (1997) GCA 61, 4915-4931. [8] Wadhwa M. 
and Crozaz G. (1995) GCA 59, 3629-3647. [9] Barrat 
J.A. et al. (2002) GCA 66 (in press). [10] Barrat J.A. et 
al. (2001) MAPS 36, 23-29. [11] Decitre S. et al. 
(2002) G 3 , 3, 1, 20 p. 
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Fig. 1. Li vs. Yb for Martian meteorites. The sources of 
the data are too numerous to be given here. 



Be/Yb 


Fig. 2. Li/Yb vs. Be/Yb for Martian meteorites, lunar 
mare basalts and picritic glasses, non-cumulate 
eucrites, and terrestrial recent lavas. The sources of the 
data are too numerous to be given here. 



Fig. 3. Li abundances and 8 7 Li profiles across a NWA 
480 pyroxene. 
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EXPLORING TRACE ELEMENT AND ISOTOPIC MIXING RELATIONSHIPS IN THE MARTIAN 
METEORITE SUITE. Lars E. Borg, University of New Mexico, Albuquerque, NM 87131, lborg@umn.edu. 


Introduction: Variations in the incompatible trace 
element and isotopic systematics of the martian meteor- 
ites are consistent with mixing of compositionally dis- 
tinct materials e.g., [1-3]. One end of the mixing array 
is defined by a component characterized by a strongly 
light-rare-earth-element depleted pattern, a low 
87 Sr/ 86 Sr ratio, and a high e Nd 143 value. In many ways 
this component is comparable to the Earth's depleted 
mantle. The second component is characterized by a 
light-rare-earth-element enriched pattern, a high 
87 Sr/ 86 Sr ratio, and a low ENd 1 43 value. This component 
could be analogous to the Earth's crust [3], or a late 
stage liquid trapped in the cumulate pile of a magma 
ocean [2]. The nature of this "enriched" component 
ultimately underscores whether the compositional 
variations observed in the martian meteorites reflect 
differentiation in the crust (i.e. assimilation fractional 
crystallization) or are inherited from their source re- 
gions at the time of melting. Below, fundamental mix- 
ing relationships are reviewed and a mechanism to 
produce the compositional variations observed in the 
martian meteorites is proposed. 

Mixing relationships: The martian meteorites display 
a large range of incompatible trace element and iso- 
topic ratios. Furthermore, variations of trace element 
and isotopic ratios correlate and fall on calculated mix- 
ing curves defined using martian meteorites as compo- 
sitional end-members [Fig. 1]. In contrast. 


incompatible trace elements and initial isotopic compo- 
sitions of the shergottites are controlled by the same 
mixing process, whereas the major elements are not. 

The whole rock isotopic systematics of the martian 
meteorites are also consistent with mixing processes, 
and in some cases offer constraints on the age of for- 
mation of the end-member components. Figure 2 is a 


Rb-Sr whole rock mixing diagram 



87 Rb/ 86 Sr whole rock 

Figure 2. Whole rock Rb-Sr isochron plot of martian mete- 
orites. The mixing line (dashed) is calculated using 
QUE94201 and Shergotty as end-members and has a slope 
corresponding to 4.49 Ga. Crosses represent 10% intervals. 
Solid line regressed through data represents an age of 4.4 Ga. 
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Figure 1. La/Yb vs. ENd' 43 • Shergotty and QUE 94201 de- 
fine the ends of mixing line. Crosses represent 10% inter- 
vals. 

major element abundances do not correlate with either 
incompatible trace element or isotopic systematics of 
the martian meteorites. This implies that variations in 
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whole rock Rb-Sr isochron plot of the martian meteor- 
ites. The meteorites define a line with a slope corre- 
sponding to an age of 4.4 Ga. A mixing line between 
QUE94201 and Shergotty defines an age of 4.49 Ga, 
suggesting that planetary differentiation occurred near 
4.4 to 4.5 Ga. A more precise age of planetary differ- 
entiation can be defined using a two-stage Nd model 
isochron [Fig. 3]. This chronometer is based on the 
l46 Sm - 142 Nd (t Vl = 103 Ma) and the 147 Sm - l43 Nd (t, A = 
106 Ga) isotopic systems and is discussed in detail by 
Borg et al. [4]. Although the model isochron is also a 
mixing line, the slope corresponds to the age of source 
differentiation of 4.51 Ga. 

Sources of martian magmas: The 87 Rb/ 86 Sr, 
147 Sm/ I44 Nd, and ,76 Lu/ 177 Hf ratios of the martian me- 
teorite sources can be estimated from their ages and 
initial and whole rock isotopic compositions by assum- 
ing they formed at 4.51 Ga. To complete this calcula- 
tion, undifferentiated Mars is assumed to have chon- 
dritic I47 Sm/ 144 Nd and 176 Lu/ 177 Hf ratios, and a 
87 Rb/ 86 Sr ratio of 0.16 estimated for bulk Mars [2]. 
The parent/daughter ratios of the martian meteorite 
sources correlate with one another [Fig 4a-b], La/Yb 
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Figure 3. Two stage model isochron constraining the time of 
martian differentiation at ~4.51 Ga. Figure, model calcula- 
tions, and data sources from Borg et al. [4]. 

ratios, and estimates of magmatic f0 2 [5], but not with 
major element abundances of the meteorites. This sug- 
gests that much of the compositional variation ob- 
served in the meteorites are inherited from their ancient 
sources. 

Meteorite sources: The compositions of the mar- 
tian meteorite sources are very similar to the composi- 
tions estimated for the sources of lunar mare basalts. 
For example, mafic lunar cumulate source regions for 
the mare basalts are estimated to have 147 Sm/ ,44 Nd ra- 
tios of 0.26-0.33 and ,76 Lu/ I77 Hf ratios of 0.046-0.072 
[6-7] and are comparable to ratios estimated for the 
QUE94201 source region ( l47 Sm/ l44 Nd = 0.285; 
176 Lu/ 177 Hf - 0.048). Likewise, KREEP is estimated to 
have ,47 Sm/ 144 Nd « 0.172 and 176 Lu/ 177 Hf = 0.014 [8] 
and is comparable to the ratios calculated for Shergotty 
source region ( I47 Sm/ ,44 Nd = 0.182; l76 Lu/ ,77 Hf = 
0.028). The 87 Rb/ 86 Sr ratios of lunar and martian 
source are difficult to compare because of different 
volatile element abundances on the two bodies. How- 
ever, the 87 Rb/ 86 Sr ratios of the lunar sources are within 
~2 of the calculated martian sources. 

Mixing models produced using the compositions of 
lunar mafic cumulates and KREEP as end-members are 
presented in Figure [4a-b]. The 87 Rb/ 86 Sr ratios of the 
lunar sources are increased by a factor of 2 to account 
for different volatile element abundances on the Moon 
and Mars. The mixing models are able to reproduce 
the compositions of the martian meteorite sources sug- 
gesting that the isotopic variations observed in the mar- 
tian meteorites reflect mixing of components formed by 
crystallization of a magma ocean. Positive correlation 



87 Rb/® 6 Sr source 
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Figure 4a-b. Parent/daughter ratios of martian meteorite 
sources calculated assuming differentiation at 4.51 Ga. Rb-Sr 
and Sm-Nd data summarized in [4], Lu-Hf data from [9]. 

between the proportion of the KREEP-Iike component 
estimated in the meteorite source regions and estimates 
of fD 2 suggests that the KREEP-like component may 
be able to oxidize reduced martian magmas, and there- 
fore could be hydrous. In this scenario, variations in 
trace element abundances, isotopic compositions, and 
magmatic f0 2 are produced as a result of melting com- 
positionally diverse magma ocean cumulates. The lack 
of correlation of major element compositions with trace 
element and isotopic compositions is likely to reflect 
the fact that the variations in major element abun- 
dances are primarily controlled by crystallization proc- 
esses after the magmas have left their source regions. 

References: [1] Shih et al. (1982) GCA 46, 2323- 
44. [2] Borg et al. (1997) GCA 61, 4915-31. [3] Jones 
(1989) Proc. 19th LPSC 465-74. [4] Borg et al. (2001) 
LPSC XXXII #1 144. [5] Herd et al. (2001) Am. Min . 
86, 1015-24. [6] Paces et al. (1991) GCA 55, 2025-43. 
[7] Beard et al. (1998) GCA 62, 525-44. [8] Warren & 
Wasson (1979) RGSP 17, 73-88. [9] Blichert-Toft et 
al. (1999) EPSL 173, 25-39. 
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Introduction: The “sampling” of Mars, both re- 
motely and by meteorites, provides evidence for a 
range of reducing and oxidizing environments. Proc- 
esses that modified the crust and mantle on Mars are 
clearly as diverse as those on the Earth with regard to 
oxygen fugacity. Increasingly compelling evidence for 
the presence of water in the near surface regime pro- 
vides a potentially oxidizing near-surface environment, 
especially when coupled with increasing hydrogen loss 
from the planet through time. Dewatering of the man- 
tle and the purported absence of global crust-mantle 
recycling perhaps accounts for the apparent reducing 
regime that dominates the source regions of many of 
the meteorite samples when compared to samples of 
the Earth’s mantle. 

The meteorite samples available must however 
have been near the surface of the planet prior to their 
entrainment in impact ejecta that eventually traveled to 
Earth. Evidence for their interaction with that surface 
regime is surprisingly sparse, or perhaps has not been 
generally recognized. (Have we been looking for the 
wrong things?) Magmatic interactions with crustal ma- 
terial such as water and potential oxidative processes 
that might accompany them have not received as much 
attention as the search for mantle signatures of these 
meteorites. Crustal processes include the crystalliza- 
tion of the magmas, hydrothermal interaction of the 
solid igneous rock with its immediate environment, 
metamorphic modification, later weathering and ero- 
sion to form sediments, and of course impact-induced 
modification. All of these crustal processes have the 
potential to alter the redox character of the samples. 
Our ability to document these processes is improving 
rapidly. Systematic studies of elemental redox state, 
particularly Fe, have the promise to extract detail of the 
changes in oxidation history that each sample has un- 
dergone. The methodologies of these studies include 
both direct and indirect approaches and a corollary of 
these studies is the necessity that consistent results be 
produced by different techniques. 

Crystallization Histories: Despite the broad geo- 
chemical interest in characterizing the Martian mantle, 
it remains obvious that the majority of the Martian 
meteorites represent igneous rocks that crystallized at 
or near the surface of the planet, i.e. in the crust. The 
likelihood of material, that crystallized in the Martian 
mantle, being emplaced in sufficient quantity to ensure 
their sampling by impacts capable of delivering the 


samples to Earth is low. In contrast, typical mafic as- 
semblages that may be ubiquitous throughout the Mar- 
tian crust should dominate the igneous meteorite suite. 
It is therefore imperative that we understand the crustal 
processes that produced the rocks we see. The phase 
equilibria as we know them strongly support such shal- 
low formation histories. We must therefore ensure that 
our studies document that crustal history as completely 
as possible. (For example, the absence of identified 
‘sedimentary-rock’ meteorites from Mars demonstrates, 
at best, an unknown selection process that prevents 
their survival as viable meteorites or at worst our in- 
ability to recognize a ‘meteowrong’ as an unfamiliar 
meteorite type.) 

The Martian meteorites differ from most other ig- 
neous meteorites as they are dominated by crystalline 
rocks. Their histories have not been obscured by the 
brecciation and mixing processes that dominate both 
lunar and the eucrite-diogenite continuum. Shock 
modification, despite the general lack of brecciation, 
remains a definitive characteristic of these samples and 
must be taken into account when assessing petrologic 
and geochemical evidence. Their crystallization histo- 
ries, began with magmatic crystallization events in 
some crustal chamber, proceeded through a sequence 
of events in the crust of Mars (about which we have 
little understanding) and ended with a major nearby 
impact event that sent the samples to Earth. The objects 
were recovered as stones with telltale shock indicators 
to alert us of their extraterrestrial provenance. 

Documenting crustal processes in Martian sam- 
ples: Petrological studies remain a fundamental tool 
for documenting the history of a sample. Extensions of 
those studies using powerful instrumental techniques 
on petrologically characterized subsamples have 
proven to be extremely valuable. The study of rock 
interaction with fluid phases that are no longer present 
is of course a classic problem for petrologists and Mar- 
tian petrology and geochemistry is beginning to deal 
with this thorny problem. 

Redox Exchange Processes: The oxidation or re- 
duction of Fe in Martian near surface environments is a 
central issue to studies of the planet and is our focus 
here. Measurement of the relative abundance of ferric 
and ferrous iron in the phases of redox exchange reac- 
tion is one of the most direct tools for documenting the 
role of that reaction in the history of the sample. It is, 
by no means, the only approach that can be used, nor 
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will it document all the interactions that may have 
taken place. Several approaches for documenting re- 
dox exchanges in the meteorites are available. 

Prior to the availability of techniques for direct 
measurement of ferric/ferrous ratios at micro- and 
nano-scales, either bulk sample measurements or very 
labor intensive, mineral separation measurements were 
necessary. Wet chemical techniques and Mossbauer 
spectroscopy define the state of that art. 

The advent of microanalytical tools with both the 
energy and spatial resolution necessary to make fer- 
rous/ferric measurements at petrologically relevant 
scales (microXANES and EELS) does not reduce the 
importance of the bulk techniques, particularly for in- 
dependent tests of instrument calibration. These micro 
and nanoscale methods have great potential but they 
suffer from the same issues that plagued other more 
familiar microanalytical techniques. The application of 
these techniques, to real problems, highlights weak- 
nesses in their theoretical underpinning that compro- 
mise both their usefulness and acceptance of their re- 
sults. In this regard, they do not differ significantly 
from other microbeam techniques such as secondary 
ion mass spectrometry for which the decades long 
search for fully quantitative theoretical models contin- 
ues. For that matter, microanalysis techniques using the 
venerable electron probe continue to benefit from sig- 
nificant theoretical improvements even after a half- 
century of development. Clearly a role remains for 
both empirical and theoretical techniques in all these 
approaches. For this reason, active exchanges of re- 
sults from these new techniques are critical. For exam- 
ple, when microXANES provides results that are the 
“odd man out” the challenge to find out why becomes 
more acute. 

Indirect approaches to redox processes in the Mar- 
tian meteorites cover the gamut of petrological re- 
search. Phase equilibria provide many of the basic rela- 
tionships necessary for calibrating redox exchange re- 
actions. The changes in phase assemblages as a func- 
tion of oxygen fugacity provide the fundamental results 
against which other approaches can be compared. Fine 
detail of the exchanges involved is sometimes difficult 
to extract. Partitioning studies of redox sensitive reac- 
tions provide some of that detail. For major and some 
minor elements such as Fe and Cr, these studies are 
intimately linked to the phase equilibrium. The ex- 
change of Fe between oxides and silicates and its 
change from ferric to ferrous in that context is both a 
phase equilibrium indicator and a sensitive redox indi- 
cator. For trace element partitioning, however, the link 
to phase equilibria may be less simple given the multi- 
ple possible hosts for many trace elements. Neverthe- 
less such studies have successfully provided usable 


indicators for redox exchange. As the number of tech- 
niques available increases, clearly the opportunity for 
inconsistency among those techniques also increases. 

Back to Mars / In samples that saw a variety of 
processes, conflict among datasets becomes more diffi- 
cult to identify. No indicator of, for example, redox 
state will be sensitive to all crustal process and faith- 
fully record them. This is particularly true of the direct 
measurements of iron redox state such as micro- 
XANES. 

Ferric and ferrous contents of co-existing minerals 
in the Martian samples can potentially be modified 
during crystallization of the primary assemblage. 
However, later interaction with hydrothermal fluids can 
radically alter the magmatic signatures. For the mo- 
ment, we have the luxury of ignoring large scale 
weathering effects in the Martian crust, as none of the 
meteorites presently is recognized to have a strong 
weathering overprint. That may change rapidly, how- 
ever. Similarly, strong metamorphic overprints compa- 
rable to those on Earth are not yet part of the gamut. 
However, the final Martian crustal event that all these 
meteorites experienced was impact and ejection from 
Mars toward the Earth. This raises a suite of issues 
that will need to be explored systematically and in de- 
tail. No aspect of a shock event requires that the target 
experience either reduction or oxidation. However, the 
environment in which the shock is applied can pre- 
sumably provide an impetus in either direction. Identi- 
fying and quantifying localized shock effects and relat- 
ing redox modification (if any) to them, is a problem 
that must be addressed. When shock is accompanied, 
or followed by, a thermal event the potential for reset- 
ting redox indicators is high. For some indicators of 
redox, these late processes are irrelevant, for others 
they will be fundamental. At present, it is simply im- 
possible to say which indicators will be most sensitive 
to each stage in the history of each meteorite sample. 
However, both the mantle and crustal history of these 
samples is probably reflected in the results of the dif- 
ferent types of redox study. 

Petrologically constrained studies that search for 
BOTH the effects of the mantle and crustal history of 
each sample remain desirable. As we look beyond the 
dominant magmatic signatures of the Martian meteor- 
ites, their crustal history should become clearer. While 
we focus on redox process here, such arguments gener- 
alize. We must begin to identify specific crustal AND 
explicitly surface processes so that the meteorites can 
continue to be useful samples of Mars when 
field/mission constrained samples of that planet’s out- 
ermost surface begin to return to Earth. 
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Introduction: The Earth is still a dynamic 
planet, with mantle convection sufficient to drive 
the tectonic plates and outer core convection 
strong enough to produce a magnetic field. Mars 
started its geological life almost as the Earth and 
had in the past a liquid core, able to generate a 
strong magnetic field. But today. Mars has lost 
this internal magnetic dynamo and its mantle is 
probably slowly convecting. The study of radio- 
genic isotopes in the Martian meteorites (SNC), 
reveals that their parent body “must be frozen” in 
its planetary evolution. 

129/132 Xe systematic: From the observed 

higher 129 Xe/ I32 Xe ratio in the Martian atmos- 
phere compared to the interior of Mars [1, 2, 3] 
we postulated an efficient extraction of the halo- 
gens from the interior into the crust. The enrich- 
ment of 129 Xe rad in the Martian atmosphere, pro- 
duced by the decay of 129 I with a half-life of 15.7 
Ma, points to a very early extraction of iodine 
into the crust; the 129 Xe ra d was released into the 
present atmosphere over geological time. In con- 
trast, the 129 Xe/ 132 Xe ratio in terrestrial rocks is 
equal or higher than the atmospheric ratio, be- 
cause of a continuous recycling of the crust and 
mantle. 

Rb-Sr isotopic systematic: All whole rock 
data of the Rb-Sr isotopes cluster in 3 groups 
close to the meteoritic isochron of 4.55 Ga. Start- 
ing with the Rb-Sr systematic of only 6 SNCs we 
postulated about 1 0 years ago the existence of 3 
isotopically distinct reservoirs on Mars, which 
remained isolated for a period of 4.3 ± 0.2 Ga [4]. 
Since then, more than 20 SNC meteorites have 
been found in hot and dry deserts. New isotope 
data from these recently recovered SNCs confirm 
the three isotopically distinct reservoirs. 

The first group, the basaltic shergottites Sher- 
gotty, Zagami and Los Angeles have relatively 
high abundances of radiogenic Sr, which might 
originate from a planetary crust enriched in in- 
compatible elements. The in situ measurements of 
the Martian surface by Viking [5], Phobos [6], 
and Mars Pathfinder [7], thousands of kilometres 
apart, reveal similar compositions of the Martian 


soils and the basaltic shergottites. Of course, we 
have to neglect in this comparison the extremely 
high S- and Cl-concentrations of the soil, which 
derived probably from volcanic exhalations [8]. 
The K contents of the soil measured by Phobos 
[6] and Mars Pathfinder [7], with 0.3 % and 
0.55 % respectively, are higher than the 0.18 % K 
content in the basaltic shergottites. This higher K 
concentration on the Martian surface might indi- 
cate crustal reservoirs, which are more radiogenic 
in Sr as found in the shergottites. However, com- 
pared to the Earth’s upper continental crust, with 
2.8 % K [9], the Martian surface has a moderate K 
inventory. 

A second group of Martian meteorites, charac- 
terized by non-radiogenic Sr, consists of the mafic 
cumulates nakhlites and Chassigny, the olivine 
rich shergottites DaG 476, SaU 005, Dhofar 019, 
and the basaltic QUE 94201 and may represent 
the depleted mantle. This and the first group are 
chemically complementary, suggesting that crust 
formation has caused the mantle depletion, which 
must have taken place during a very early proc- 
ess. This can be derived from the primitive Sr iso- 
topes and the excess Nd-142, the daughter prod- 
uct of the extinct Sm-146, found in Chassigny, 
the nakhlites, SaU 005, DaG476, and QUE 94201 
(E42 group). The observed correlation of 142 Nd 
with 182 W in SNC meteorites by Lee and Halliday 

[10] points to the same incompatible behaviour of 
Nd and W during magmatic processes and an 
early core formation (about 30 Ma). Core forma- 
tion must have occurred concurrently with the 
rapid accretion of Mars, which is inferred from 
excess 182 W in Martian meteorites. 

A third meteorite group with intermediate Sr 
isotopic composition, represented by the olivine 
gabbroic LEW 88516, ALHA77005, and Ya- 
mato793605, might originate from a primitive, 
unfractionated mantle. However a mixing process 
of crust and depleted mantle cannot be excluded 

[ 11 ] . 

The lack of plate tectonics can keep the obvious 3 
distinct reservoirs untouched. On Earth, plate tec- 
tonics might be responsible for the homogeniza- 
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tion of the terrestrial Rb-Sr reservoir. Today, only 
one Rb-Sr isotopic reservoir exists on Earth. The 
terrestrial crust and depleted mantle are developed 
from an isotopic homogeneous mantle at about 2 
Ga while on Mars crust and depleted mantle de- 
veloped at 4.3 Ga. 

Sm-Nd systematic: The Sm-Nd systematic of 
terrestrial samples reveals also a crust-mantle dif- 
ferentiation from a chondritic Sm/Nd reservoir 2 
Ga ago. In the terrestrial samples of crust and 
mantle the Sm/Nd ratio changes by a factor of 2. 
The SNC meteorites, however, have strongly frac- 
tionated Sm/Nd ratios and large variation of the 
143 Nd/ 144 Nd ratio particularly among the rocks 
from a depleted mantle reservoir: Nakhla, Chas- 
signy, DaG 487, SaU 005, Dhofar and QUE 
94201. In the Sm-Nd isotope plot only the data of 
the minerals and the whole rock of ALHA 84001 
fall on the 4.55 Ga meteoritic internal isochron. 
All other SNCs are far away from this isochron 
due to their strongly fractionated Sm/Nd ratio. In 
the Sm-Nd system we cannot find the 3 distinct 
initial isotopic reservoirs as found for the case of 
the Rb-Sr systematic. 

Pb-Pb systematic: The Pb isotopes of all 
measured SNCs show a similar pattern as the Sr 
isotopes. The present study of the Pb isotope sys- 
tematic might give the strongest indication for an 
early differentiation of Martian mantle and crust. 
The initial Pb from hand-picked ultra clean pla- 
gioclase separates Los Angeles, Shergotty, Za - 
garni, and the new shergottite Dhofar 378, from 
the enriched crustal reservoir, and of nakhlites 
and SaU 005, from the depleted mantle reservoir, 
plot close to the 4.5 Ga Pb-Pb isochron. The con- 
formity of the U-Pb and Rb-Sr isotopic systemat- 
ics reflects similar magmatic fractionation behav- 
iour of Rb and U during the evolution of curst and 
mantle. 


relation of radiogenic 182 W with radiogenic 142 Nd 
by [10] points also to a close relationship between 
core formation and mantle melting in the first ~30 
Ma after formation of the Mars. Contrary to Mars, 
on Earth a relationship exists between the Sr- and 
Nd-isotope systems. Because of the remixing of 
the crust into the Earth’s mantle and the homog- 
enization of the isotopes through plate tectonics a 
“common isotopic reservoir” exists from which 
about 2 Ga ago crust and depleted mantle differ- 
entiated. 

References: [1] Becker R. H. and Pepin R. O. 
(1984) EPSL, 69, 225-242. [2] Swindle T. D. et 
al. (1986) GCA, 50, 1001-1015. [3] OttU. and 
Begemann F. (1985) Nature 317 , 509-512. 

[4] Jagoutz E. (1991) Space Sci. Rev., 56, 13-22. 

[5] Clark B. C. et al. (1982) J. Geophys . Res., 87, 
10,059-10,067. [6] Surkov Yu. A. et al. (1989) 
Nature, 341, 595-598. [7] Wanke H. et al. (2001) 
Space Sci. Rev., 96, 317-330. [8] Clark B. C. 
(1993) GCA, 57, 4575-4581. [9] Taylor S. R. and 
McLennan S. M. (1995) Rev. Geophys., 33, 241 - 
265. [10] Lee D.-C. and Halliday A.N. (1997) 
Nature 38 , 854-857. [11] Borg L. E. (2002) GCA, 
66, 2037-2053. 


Summary: The 12 V 29 Xe (T 1/2 16 Ma) and the 
l46 Sm- l42 Nd (T 1/2 103 Ma) isotope systems indi- 
cate a rapid accretion and a very early formation 
of the crust with its enrichment of volatile and 
highly incompatible lithophile elements. The ab- 
sence or at most very limited plate tectonic activ- 
ity on early Mars excludes an extensive crustal 
recycling back to the mantle and preserves the 
Rb-Sr and U-Pb isotopes systems derived from 
the early crustal differentiation. The observed cor- 
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MOSSBAUER SPECTROSCOPY OF SNC METEORITES. M. Darby Dyar, Dept, of Earth and Environment, 
Mount Holyoke College, 50 College St., South Hadley, MA 01075. 


Introduction: With the inclusion of Mossbauer 

spectrometers in the instrument payloads of both MER 
landers in the upcoming Athena mission, it is timely to 
re-examine Mossbauer data on SNC meteorites and fill 
the gaps in our state of knowledge. Much of this work 
was done by Roger Bums ten years ago on whole rock 
samples [1-6], but the research effort was cut short in 
1993. In this project, we are summarizing existing 
Mossbauer spectra of SNC meteorites and showcasing 
new results from a state-of-the-art spectrometer cus- 
tomized to accommodate smaller (less than 10 mg) 
samples. The results demonstrate the problems inher- 
ent in drawing conclusions about individual phases 
chemistry on the basis of Mossbauer spectra of whole 
rock samples. 

Previous Work: Previous work on these meteor- 
ites is summarized in Table 1. The only workers to 


Table 1. Mossbauer Spectroscopy of SNC Samples 


Meteorite 

Analysis 

%Fe 3+ 

ALHA 77005 

WRfl, 2, 5] 

1.7-1. 8 


olivine separate [11, 12] 

4.5 

Chassigny 

WR[3, 4, 51 

2.3 

EETA 79001 

WRn,5, 61 

0.5-1. 9 

Govemador 

Valadares 

WR[5, 6] 

n.s. 

Lafayette 

WRr5, 61 

n.s. 

Nakhla 

WRri,5, 6,8, 131 

1.6 


pyroxene separate [14] 

2.0 

Shergotty 

WR[5, 61 

n.s. 

Zagami 

WRr5, 6, 7, 8, 10, 131 

n.s. 


pyroxene separate [13] 

<1.0 


WR = whole rock analysis, n.s. = not specified. 


extract and analyze mineral separates were Ostertag et 
al. (olivine) [11,12], Vieira et al. [13], and Vistisen 
[14] (pyroxene), and these data show only very small 
percentages of Fe 3+ in olivine and pyroxene. However, 
none of the previous workers has related the Fe 3+ /EFe 
of the bulk separates to the Fe contents and modes of 
individual minerals present, so Fe 3+ /£Fe ratios of indi- 
vidual phases could not be assessed. 

Methods: In the current study, roughly 200 mg 

samples of each of nine SNC meteorites were obtained 
from JSC, NMNH, and the Vatican (Table 2). Whole 
rock Mossbauer spectra were acquired on a 10 mg chip 
of each sample, and the remaining mass was crushed to 
allow handpicking of mineral separates for phases con- 
stituting more than 30% of the mode. 


For all Mossbauer analyses, samples were then 
crushed under acetone to a fine powder, mixed with 
sucrose to avoid preferred orientation, and run on a 
WEB Research spectrometer at room temperature. 
Data were processed using the WMOSS software 
package, which allows use of multiple models for line 
shape. 


Table 2. Mineral Modes (%) for Samples Studied 


Meteorite 

Px 

Fsp* 

Ol 

Ox 

# 

Ref. 

ALHA 77005 

35 

8 

50 



[15] 

ALHA 84001 

97 

1 


2 


[161 

Chassigny 

5 

2 

92 

1 


[17] 

LEW 88516 

35 

8 

50 


5 

[18] 

Los Angeles 

40 

45 

3 

3 

9 

[191 

Nahkla 

80 


10 


10 

[20] 

QUE 94201 

44 

46 


2 

8 

[21] 

Shergotty 

70 

24 


2 

4 

[22] 

Zagami 

70 

24 


3 

3 

[22] 


*Fsp or maskelynite, # = other phases such as glass 


Results: Whole rock spectra of ALHA 84001, 

LEW 88516, QUE 94021, and Los Angeles have not 
previously been published, so they are included here as 
Figures 1-4. Spectra of the other five samples in Table 
2 have also been acquired. The QUE 94201 spectrum 
is typical of orthopyroxene, with a small contribution 
(about 5% of the total Fe) as ilmenite. Both Los An- 
geles and ALHA 84001 show spectra typical of clino- 
pyroxene, and the LEW 88516 spectrum shows a mix 
of Fe in olivine and pyroxene. 

These spectra show good examples of how confus- 
ing the interpretation of Mossbauer spectra of mixtures 
can be! For example, in the spectrum of LEW 88516, 
60% of the area is in the outer olivine doublet, while 
the remaining 40% of the area can be assigned to Cpx. 
Normalizing the mineral mode suggests that the sample 
should be 59% olivine and 41% pyroxene. The fact 
that these numbers match is merely coincidence be- 
cause the percentage of Fe atoms in the two minerals is 
not the same. The olivine is 27 wt%FeO (0.60 apfu) 
and pyroxene is roughly 10-16 wt% FeO (0.28-0.51 
apfu) [23]. The mineral modes reflect the variations in 
effective volume of the phases present. Thus, based on 
the molar volumes of the unit cell of olivine (Z=4, vol- 
ume-290 A 3 ) and pyroxene (Z=4, volume= : 449 A 3 ), the 
percentage of Fe atoms based on the mineral modes 
should be 27% in pyroxene and 73% in olivine if all 
the Fe is Fe 2+ . Thus, the predicted areas are not a bad 
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match to the observed peak areas noted above, espe- 
cially considering the variations in pyroxene chemistry 
upon which these calculations are made (and the fact 
that recoil-free fractions for these two minerals are 
similar but slightly different). 

Two other examples are also useful [24]. Nahkla 
has a mode that is 80% cpx and 10% olivine, but these 
data combined with the chemical compositions suggest 
a Mossbauer spectrum with doublet areas of 75% cpx 
and 25% olivine (close to those actually observed in 
[1] and to our own spectrum). This change occurs be- 
cause the olivine is Fe-richer than the coexisting py- 
roxene. ALHA 77005 has a mode with 50% olivine 
and 35% pyroxene, but its predicted Mossbauer spec- 
trum has the proportions 62% olivine and 38% pyrox- 
ene. 

Conclusions: If mineral modes and the Fe contents 
of the component phases are known, Mossbauer dou- 
blet areas can be predicted with some confidence. 
However, this relationship does not work in reverse 
because there are three variables of which two are un- 
known. Doublet areas in the Mossbauer spectrum of a 
mixture of minerals cannot be used to estimate modal 
mineralogy without known chemical compositions In 
cases where the bulk rock is a mixture of, say, hematite 
and magnetite (i.e., phases where the chemical compo- 
sitions are relatively constant), modal estimates can be 
confidently made based on Mossbauer spectra of mix- 
tures. Unfortunately, in cases where the bulk rock 
contains silicates with variable Fe contents (i.e., nearly 
any silicate!) it is impossible to deduce modes from 
Mossbauer doublet areas. Interpretation of MER data 
for which modes will not be known will be (regretably) 
complicated by this simple constraint. 

References: [1] Bums R.G. and Solberg T.C. (1988) LPSC 
XIX ; 146. [2] Burns R.G. (1989) NASA Tech. Mem., Rpts . PG&G 
Prog. 1988 ; 211-212. [3] Bums R.G. and Martinez S.L. (1990) 
LPSC XXI, 147-148. [4] Bums R.G. and Martinez S.L. (1990) LPI 
Tech. Rpt. 90-06, 91. [5] Bums R.G. (1991) NASA Tech. Mem., 
Rpts. PG&G Prog., 4300, 259. [6] Bums R.G. (1991) LPSC XXII, 
157. [7] Agervist D.P. and Vistisen L. (1993) LPSC XXIV, 1. [8] 
Agerkvist D.P. et al. (1994) LPSC XXV, 1. [9] Morris R.V. (1997) 
Wrkshp. on Early Mars, #3047. [10] Morris R.V. et al. (1998) LPSC 
XXIX, #1326. [11] Ostertag R. et al. (1983) Meteoritics, 18, 368. 
[12] Ostertag R. et al. (1984) EPSL, 67, 162. [13] Vieira V.W.A. et 
al. (1986) Physica Scripta , 33, 180. [14] Vistisen L. et al. (1990) 
Meteoritics, 25, 417. [15] Mason B. (1981) Ant. Meteor. Newsl, 4, 
12. [16] Mason B. et al. (1992) Smithson. Contrib. Earth Sci., 30, 
17. [17] Prinz et al. (1974) Meteoritics, 9, 393. [18] Satterwhite and 
Mason (1991) Ant. Meteor. Newsl . , 14, 19. [19] Rubin et al. (2000) 
Geology, 28, 1011. [20] Nakamura et al. (1982) GCA, 46, 1555. 
[21] Harvey et al. (1996) LPSC XXVII, 497. [22] Stolper and 

McSween (1979) GCA 43, 1475. [23] Harvey et al. (1993) GCA, 51, 
4769. [24] Calculations are based on unpublished EMPA data 
courtesy of J.S. Delaney. 

Figures 1-4 (at right). Room temperature Moss- 
bauer spectra of SNC meteorites. 
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PETROGENESIS OF OLIVINE-PHYRIC SHERGOTTTTES SAYH AL UHAYMIR 005 AND ELEPHANT MORAINE 
A79001 LITHOLOGY A. Cyrena Anne Goodrich, Hawaii Institute of Geophysics and Planetology, University of 
Hawaii at Manoa, Honolulu, HI 96822, USA. cyrena@higp.hawaii.edu. 


Introduction: Martian meteorites Sayh al Uhaymir 
(SaU) 005 and lithology A of EETA79001 (EET-A) 
belong to a newly emerging group of olivine-phyric 
shergottites [1]. Previous models for the origin of such 
shergottites have focused on mixing between basaltic 
shergottite-like magmas and lherzolitic shergottite-like 
material [2-7]. Results of this work, however, suggest 
that SaU 005 and EET -A formed from olivine-saturated 
magmas that may have been parental to basaltic 
shergottites. 

Description: SaU 005 and EET-A have porphyritic 

textures of large (up to -3 mm) olivine crystals (-25% in 
SaU 005; -13% in EET-A) in finer-grained 

groundmasses consisting principally of pigeonite 
(-50% in SaU 005; -60% in EET-A) and plagioclase 
(maskelynite), with <7% augite. Low-Ti chromite 

occurs as inclusions in the more magnesian olivine, and 
with chromian ulvospinel rims in the more ferroan 
olivine and the groundmass. Crystallization histories 
for both rocks were determined from petrographic 
features (textures, crystal shapes and size distributions, 
phase associations, and modal abundances), mineral 
compositions, and melt compositions reconstructed 
from magmatic inclusions in olivine and chromite. 
Results: The following observations indicate that the 
chromite and most magnesian olivine (Fo 74-70 in SaU 
005; Fo 81-77 in EET-A) and pyroxenes (low-Ca 
pyroxene [Wo 4-6] of mg 77-7 4 and augite of mg 78 in 
SaU 005; orthopyroxene [Wo 3-5] of mg 84-80 in EET- 
A) in these rocks are xenolithic. 1) Olivine crystal size 
distribution (CSD) functions show excesses of the 
largest crystals (whose cores comprise the most 
magnesian compositions), indicating addition of 
phenocrysts or xenocrysts [8]. 2) The most magnesian 
low-Ca pyroxenes show near-vertical trends of mg vs. 
A1 2 0 3 (Fig. 1) and Cr 2 0 3 , which suggest reaction with a 
magma. 3) In SaU 005, there is a gap in augite 
composition between mg 78 and mg 73. 4) Chromite 
cores of composite spinel grains are riddled with 
cracks, indicating that they experienced some physical 
stress before being overgrown with ulvospinel. 5) 
Magmatic inclusions are absent in the most magnesian 
olivine, but abundant in the more ferroan, indicating 
slower growth rates for the former. 6) The predicted 
early crystallization sequence of the melt trapped in 
chromite (the earliest phase) in each rock produces its 
most magnesian olivine-pyroxene assemblage (Fig. 
2b, d). However, in neither case is the total 

crystallization sequence of this melt consistent with the 
overall crystallization history of the rock or its bulk 
modal mineralogy (Fig. 2a, c). 


Further, the following observations indicate that in 
both SaU 005 and EET-A the fraction of solid xenolithic 
material is small (in contrast to previous models for 
EET-A [2-5]), and the bulk of the rock formed by 
continuous crystallization of a single magma (possibly 
mixed). 1) CSD functions and correlations of crystal 
size with composition show that most of the olivine (Fo 
69-62 in SaU 005; Fo 76-53 in EET-A) formed by 
continuous nucleation and growth. 2) Groundmass 
pigeonites are in equilibrium with this olivine, and show 
continuous compositional trends that are typical for 
basalts (Fig. 1). 3) The CSD function for groundmass 
pigeonite in EET-A indicates continuous nucleation 
and growth [9]. 4) The melt trapped in olivine of Fo 76- 
67 in EET-A has a predicted crystallization sequence 
similar to that inferred for the bulk of the rock and 
produces an assemblage similar to its modal 
mineralogy. 5) Melt trapped in late olivine (Fo -64) in 
SaU 005 has a composition consistent with the inferred 
late crystallization history of the rock. 

The conclusion that only a small fraction of either 
SaU 005 or EET-A is xenolithic implies that both rocks 
lost fractionated liquids late in their crystallization. 
This is supported by: 1) olivine CSD functions, which 
show a drop in nucleation rate at high degrees of 
crystallization [8]; 2) high pigeonite/plagioclase ratios; 
and 3) low augite contents. For EET-A, this 
fractionated liquid may be represented by EET -B. 

All magmas involved in the petrogenesis of SaU 005 
and EET-A were olivine-saturated. These olivine- 
phyric shergottites did not form by mixing of basaltic 
shergottite-like magmas and lherzolitic shergottite-like 
material. The magmas from which they formed may, 
however, have been parental to basaltic shergottites. 
References: [1] Goodrich C.A. (2002) MAPS , in press. 
[2] Steele I.M. and Smith J.V. (1982) Proc. 13 th LP SC, 
JGR 87, A375. [3] McSween H.Y., Jr and Jarosewich E. 
(1983) GCA 47, 1501. [4] Wadhwa M. et al. (1994) GCA 
58, 4213. [5] Mittlefehldt D.W. et al. (1999) MAPS 34, 
357. [6] Mikouchi T. et al. (2001) MAPS 36, 531. [7] 
Wadhwa M. et al. (2001) MAPS 36, 195. [8] Marsh B.D. 
(1988) J. Petrol. 39, 553. [9] Lentz R.C.F. and McSween 
H.Y., Jr (2000) MAPS 35, 919. [10] Berkley J.L. and 
Treiman A.H. (2000) LPS 31, #1729. [1 1] Berkley J.L. et 
al. (1999) LPS 30, 1588. 
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Fig. 1. AI 2 0 3 content (wt.%) vs. mg for low-Ca pyroxenes in SaU 005 (a) and EET-A (b). For SaU 005 open red squares = opx 
inclusions in the most magnesian (Fo 74-70) olivine; black squares = groundmass pigeonite. For EET -A open red squares = opx cores 
of megacrysts; open green diamonds = opx cores in the exceptional megacryst X-14 [10, 11]; black squares = pigeonite in coronas 
around opx cores and groundmass. Most magnesian low-Ca pyroxenes in both rocks show near-vertical trends, which cannot be 
explained by closed system evolution and suggest reaction with a magma. 


Melt in Chromite in SaU 005 



Modal Abuadancei 
in SatJ QOS 





Fig. 2. Predicted crystallization sequences for melt trapped in low-Ti chromite in SaU 005 (a,b) and EET-A (c,d). In both cases, this 
melt produces the most magnesian phase assemblage in the rock (b,d), but is not consistent with its overall crystallization history or 
modal mineralogy (a,c). Areas indicated by boxes in (a,c) are shown at larger scale in (b,d). 
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MARS METEORITE STATISTICS AND THE MARTIAN UPLANDS. W. K. Hartmann, 
Planetary Science Institute, 620 N. 6 th Avenue, Tucson AZ 85705-8331. 

A significant fraction of Mars, perhaps 30-50%, is covered by heavily cratered uplands. The high 
crater densities constrain their ages to be >2 Gy or even 3 Gy. A somewhat smaller fraction is 
covered by young lava plains with crater retention ages <1.5 Gy, and even < 0.5 Gy over wide 
areas of Tharsis, Amazonis, and Elysium. These young ages for basaltic lava plains agree with 
ages of Martian basaltic meteorites. 

The problem is that all but one of the 4 to 8 sample sites of SNCs lie in the age range of 0.2 to 
1 .3 Gy, and the other one is a sample of primordial crust of 4.5 Gy age. These data suggest that 
the Martian uplands do not efficiently launch meteorites. 

MGS and Odyssey data imply that the upper latitude Hesperian/Noachian Martian uplands, older 
than 2 to 3 Gy, are impregnated with near-surface ice. Furthermore, direct geometric 
considerations of crater areal coverage show that they have been pulverized into megaregolith, 
although that material itself may be recemented by salts, carbonates, and ice. 

Thus, the Martian uplands probably do not launch meteorites as effectively as the young basaltic 
plains, and those that are launched may be rich in rock types that are unfamiliar as meteorites. 
These considerations should affect future meteorite collection strategies. 
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MARTIAN BASALT OXYGEN FUGACITY AND GEOCHEMISTRY: IMPLICATIONS FOR A HETEROGENEOUS 
MARTIAN MANTLE. C. D. K. Herd, Lunar and Planetaiy Institute, 3600 Bay Area Blvd., Houston, TX 77058, 
h e r d@ipiusra.edu . 


Introduction: The oxygen fugacity (J0 2 ) of the mar- 
tian basalts (a.k.a. basaltic shergottites) has recently 
been recognized as important in understanding their 
petrogeneses [1]. The correlation between oxygen fu~ 
gacity and certain geochemical parameters (initial Sr 
and Nd isotopic comp ositions, La/Yb ratios) indicates 
the presence of long-term incompatible -element en- 
riched and depleted reservoirs, which are oxidized and 
reduced, respectively. Oxygen fugacity variations in 
the martian basalts have been modeled based on as- 
similation of oxidized crustal (“crust-like”) material by 
reduced, mantle -derived magmas [1]. However, there 
remains the possibility that the oxygen fugacity of the 
martian basalts is controlled by the rature of their 
source regions [1,2]. Here I review existing oxygen fu- 
gacity data for the martian meteorites, present some 
speculative estimates for newer martian basalts, and 
explore the heterogeneous mantle model in more detail. 

Oxygen fugacity estimates: a review. Estimates 
based on compositions of Fe-Ti oxides are presented 
by [3], and an estimate for Dar al Gani 476 using olivine- 
pyroxene-spinel (Ol-Px-Sp) models is given by [1]. Oxy- 
gen fugacity has also been estimated by [4], using 
£WD G d ratios in pyroxene. Since this method utilizes 
the Eu and Gd concentrations of earliest-formed pyrox- 
ene (augite or pigeonite), it represents a better means of 
estimating magmatic oxygen fugacity than Fe-Ti oxides, 
which are subject to subsolidus equilibration. Initial 
estimates using this method differed greatly from the 
Fe-Ti oxide methods [1,4]. However, [1] and [3] demon- 
strated that Fe-Ti oxides in the martian basalts recorded 
near-magmatic J0 2 conditions, and [1] argued that the 
discrepancy between the methods was due to the lack 
of proper D Eu /D Gd calibration. Since then, other authors 
[5,6] have calibrated this model with martian basalt 
compositions. The estimates based on these calibra- 
tions are in better agreement with Fe-Ti oxide estimates, 
supporting minimal subsolidus equilibration of Fe-Ti 
oxides. 

The observation that J0 2 in the martian basalts var- 
ies by 3 log units or more is supported by all of the 
available methods. General agreement between the 
methods suggests that the oxygen fugacity of a given 
martian basalt was set before crystallization, and did 
not change significantly over the temperature history 
of the rock [6]. 

Oxygen fugacity estimates: new. At least 7 basaltic 
shergottites have been described since [1,3] estimated 
oxygen fugacity from Fe-Ti oxides and Ol-Px-Sp. Of the 


newly discovered basalts, SaU 005 (and pairs), NWA 
1068/1 1 10, NWA 1 195 and Dhofar 019 are texturally the 
most similar to DaG 476 (and pairs) and EET 79001 A; 
that is, characterized by olivine megacrysts in a fine- 
grained groundmass [7,8,9,10]. The Ol-Px-Sp method 
can be used to estimate the/^ of these samples if equi- 
librium between olivine, low-Ca pyroxene and Cr-spinel 
can be established, as was done for DaG 476 [I]. The 
other new basalts, namely NWA 480, NWA 856 and 
Dhofar 378, are texturally similar to Shergotty, Zagami, 
Los Angeles and QUE 94201, as noted by an absence 
of olivine megacrysts [11,12,13]. These contain Fe-Ti 
oxides, but in some cases, significant oxyexsolution has 
occurred (e.g., NWA 480 [11]). Regardless, the pros- 
pects are good for estimating J0 2 using the range of 
methods now available. 

Preliminary estimates of j0 2 can be made using pub- 
lished data for the new martian basalts. The composi- 
tions of oxides in the SaU 094 basalt [7] (a SaU 005 pair) 
indicate very reduced crystallization conditions. Using 
the composition of the most Mg-rich pigeonite and 
olivine cores, and a low-Ti chromite after [7] yields (a 
rather rough) estimate of IW - 1 .7. The compositions of 
oxides in Dhofar 019 suggest J0 2 near the IW buffer 
[10]. In contrast, the J0 2 of NWA 856 is greater than or 
equal to Shergotty, based on Eu/Gd ratios of pyroxene 
cores [14]. 

Oxygen fugacity and LREE enrichment: SaU 005 
and Dhofar 019 appear to be reduced, like QUE 94201 or 
DaG 476. The LREE-depleted bulk REE pattern of SaU 
005 is nearly identical to that of DaG 476 [15], and Dho- 
far 019 is thought to have been derived from a LREE- 
depleted magma [10]. In contrast, NWA 856 is oxidized, 
and has an accordingly LREE-enriched pattern, similar 
to that of Shergotty [14]. 

From these three examples, it appears that the rela- 
tionship between LREE enrichment and oxygen fuga- 
city [1] holds for some of the newer basalts, with more 
reduced basalts being LREE-depleted and more oxi- 
dized basalts being LREE-enriched. If the same rela- 
tionship holds for the remaining basalts, then their rela- 
tive redox conditions can be predicted based on chon- 
drite-normalized La/Yb ratios. For example, Dhofar 378 
has La/Yb ~ 0.9 [13], and would be expected to be oxi- 
dized like Shergotty or Zagami. Likewise, NWA 1068 
has a REE pattern similar to Shergotty and Zagami [14], 
and therefore would be expected to have similar oxygen 
fugacity. 
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NWA 480 does not appear to follow this trend. 
La/Yb is ~ 0.8 [1 1], and so the J0 2 would be expected to 
lie between that of EET 79001 and Shergotty/Zagami. 
However, Ihe reported mineralogy of NWA 480 does 
not support an oxidized basalt. Oxide and Ol-Px-Sp 
compositions after [11] suggest J0 2 ~ IW, like QUE 
94201. The bulk incompatible trace elements of NWA 
480 suggest it is more like the lherzolitic shergottites, 
whereas the compatible trace elements suggest an af- 
finity to Shergotty and Zagami [11]. 

Heterogeneous mantle model: Oxygen fugacity 
does not correlate with bulk composition [1], or with 
texture (e.g., presence/absence of olivine megacrysts). 
The latter apparently holds true for the new martian 
basalts as well as the old: the “reduced” group includes 
QUE 94201, as well as the olivine-bearing DaG 476, SaU 
005 and Dhofar 019. The “oxidized” group hcludes 
Shergotty, Zagami, Los Angeles, and NWA 856. None 
of these are olivine-bearing, but NWA 1068 may belong 
to this group, based on its REE pattern. This decoup- 
ling of the major element compositions from trace ele- 
ment and isotopic compositions, along with the evi- 
dence that the oxygen fugacity of a given basalt is set 
before crystallization, argues for mantle source mixing 
to explain the variations in oxygen fugacity and geo- 
chemistry of the martian basalts. 

Borg [2] argues that the Rb-Sr and Sm-Nd isotopic 
compositions of the martian meteorites can be used to 
calculate parent/daughter ratios of their sources. He 
demonstrates that these ratios correlate with oxygen 
fugacity, supporting the idea that the J0 2 of a given 
martian basalt is inherited from its source region. 

Implicit in this model is that the long-term incom- 
patible element enriched component in the martian ba- 
salts (the “crust-like” component of [1]) is located in a 
source region (possibly the mantle). Indeed, [2] mo dels 
the range of parent/daughter source ratios of the mar- 
tian meteorites using isotopic compositions of lunar 
mafic mantle and KREEP. Therefore, it is appropriate to 
consider the oxidizing characteristics of martian 
KREEP-like material, formed possibly (by analogy with 
the Moon) as the result of the crystallization of a mar- 
tian magma ocean. 

As outlined by [1], the long-term incompatible ele- 
ment enriched component in the martian basalts must 
be enriched in LREE relative to HREE, be more radio- 
genic in terms of Sr and less radiogenic in terms of Nd. 
It must have either an approximately basaltic composi- 
tion, or sufficiently small volume so as not to shift the 
bulk composition away from basalt. Also, the Rb-Sr and 
Sm-Nd isotopic systematics are best e«plained if the 
enriched component is ~4.5 Ga old. All of these charac- 
teristics can be achieved by crystallization of the last 
dregs of a magma ocean early in Mars’ history. By 


analogy with lunar KREEP, the bulk major element com- 
position of such material would not be expected to be 
significantly different from basaltic. 

The additional observation of increasing oxygen 
fugacity with increasing involvement of the KREEP-like 
component suggests a significant difference between 
the martian KREEP-like component and lunar KREEP: 
that the martian component has higher Fe 3+ /EFe, or is 
hydrous. 

It is difficult to envision the crystallization of highly 
oxidized, KREEP-like material from a magma ocean, 
since that material would have to have an oxygen fuga- 
city close to that of the hematite-magnetite buffer [1], 
which is three orders of magnitude more oxidizing than 
the most oxidized basalts on Earth (Earth has the ad- 
vantage of subduct ion). Alternatively, if the KREEP- 
like material is comprised of hydrous material such as 
amphibole or phlogopite, then it may have formed as a 
result of crystallization of a hydrous magma ocean. 
This is possible, although it is a concept not yet con- 
sidered (nor modeled) for Mars. Very little water is re- 
quired to oxidize a basalt from - IW to IW + 3 [1]. How- 
ever, it is unlikely that degassing is consistent and effi- 
cient enough to cause the same amount of oxidation for 
a given degree of involvement of the KREEP-like mate- 
rial. Whether degassing occurred in the more oxidized 
basalts is the subject of ongoing debate [16, 17]. A 
combination of the two mechanisms is more likely for 
the heterogeneous mantle model. 

Future work: Determining the J0 2 and isotopic 
compositions of the new martian basalts (especially the 
enigmatic NWA 480), as well as attempting the same for 
the cumulate martian meteorites will help determine if 
there are two or more main mantle sources. Detailed 
examination of correlations with stable isotopes and 
noble gases will further aid in discerning the multiple 
components that mix to make the SNCs. 

References: [1] Herd C. D. K. et al. (2002) GCA , 66, 
2025-2036. [2] Borg L.E. (2002) This meeting . [3] Herd C. 
D. K. et al. (2001) Am. Mineral, 86, 1015-1024. [4] 
Wadhwa M. (2001) Science, 291 , 1527-1530. [5] Mus- 
selwhite D. S. et al. (2002) This meeting. [6] McCanta 
M. C. and Rutherford M. J. (2002) This meeting. [7] 
Gnos E. (2002) M&PS, 37, 835-854. [8] Barrat J. A. et al. 
(2002) IPS XXXIII , #1538. [9] Irving A. J. (2002) M&PS, 
37, A69. [10] Taylor L. A. (2002) M&PS, 37, in press. 
[11] Barrat J. A. et al. (2002) M&PS, 37, 487-500. [12] 
Jambon A. (2001) M&PS, 36, A90. [13] Dreibus G. et al. 
(2002) M&PS, 37, A43. [14] Wadhwa M. and Crozaz G. 
(2002) M&PS, 37, A145. [15] Dreibus G. et al. (2001) 
M&PS, 36, A50. [16] Lentz R. C. F. et al. (2001) GCA, 
65, 4551-4565. [17] Herd C. D. K. et al. (2002) M&PS, 
37, A62. 
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The existence of a planet-wide early magma ocean on 
Mars is supported by a growing base of petrochemical and 
geophysical observations 1) The parent liquids to the SNC 
meteorites are significantly depleted in A1 2 0 3 and CaO 
relative to terrestrial basalts [1]. Only terrestrial komatiites, 
the products of more than 30% melting of the Archean 
mantle and boninites, wet melts of the mantle wedge in 
island arc regions, have similar low A1 2 0 3 and CaO con- 
tents. Mare basalts and picrite glasses on the Moon have 
similar geochemical depletions, and the major element 
compositions of very low Ti mare basalts bear a striking 
resemblance to the Shergotty parent magmas. What these 
terrestrial and lunar magmas have in common is that the 
parent magmas last equilibrated with a mantle severely 
depleted in magmaphile elements. The boninites and mare 
basalts, in particular, last coexisted with a mantle residue of 
olivine and orthopyroxene. In the lunar case the mantle was 
a product of crystallization from a magma ocean whereas 
the harzburgite parent mantle for boninites was a residuum 
to previous melting events that eliminated diopside from 
the mantle. 2) 182 W and 142 Nd anomalies date the fractiona- 
tion of the core and mantle, respectively, within about 50- 
100 million years of the origin of the solar system [2] [3]. 
The large heavily crated Martian crust and the absence of 
large scale recycling suggests strongly that the crust was 
also a product of this ancient global differentiation and has 
experienced only modest volcanic activity, particularly in 
the southern hemisphere, in subsequent epochs. Whole rock 
Rb-Sr systematics appear to record this planet wide differ- 
entiation at about 4.5 Ga [4] 3) The l43 Nd composition of 
the Martian mantle is significantly more depleted than the 
terrestrial mantle and even the cumulate source regions of 
mare basalts on the Moon [5]. Only Archean lithosphere on 
earth has the extreme positive and negative epsilon values 
so characteristic of the Martian mantle [6). Continental 
lithosphere, by definition, is stable and has withstood the 
homogenizing effects of mantle convection. The extreme 
epsilon values reflect ancient depletion events and subse- 
quent metasomatic perturbations. The data is consistent 
with the early differentiation of a Martian magma ocean 
producing a buoyant crust, dense core and a complemen- 
tary stratified cumulate mantle. The stratified cumulate is 
likely to be gravitationally unstable, at least, in the shal- 
lowmost stratigraphic levels where more iron-rich cumu- 
lates overlie dense magnesian cumulates [7]. Under these 
unstable conditions, solid state differentiation would have 
carried dense, iron-rich and relatively cool cumulates into 
the Martian interior ultimately resulting in a lower mantle 
that is denser and compositional more evolved than the 
upper mantle. This lower mantle would also contain vary- 
ing amounts of heat producing radioactive elements. 

A gravitationally stable compositionally stratified 
Martian mantle explains some important characteristics of 
Mars. 1) Mars, like the Moon, is isotopically much more 
heterogeneous than the Earth [3,5]. While one might argue 
that the small size of the Moon or the absence of water to 
promote solid state creep inhibits thermal convection, such 
arguments do not apply to Mars. 2) The presence of 142 Nd 


anomalies in some SNCs require that their source regions 
preserved their initial compositional heterogeneities [5]. 3) 
The preservation of whole rock Rb-Sr isochrons and the 
ancient ,82 W and 142 Nd anomalies argue against crustal 
recycling subsequent to stabilization of the magma ocean 
cumulates [3]. 4) The isotopic evidence of early differentia- 
tion in the Martian mantle, in contrast to the evidence from 
the terrestrial mantle, indicates that the Martian mantle has 
successfully resisted homogenization by thermal convec- 
tion. Thermal convection is inhibited by a stably composi- 
tionally zoned mantle [7], 

An alternative hypothesis is that the SNC meteorites 
are derived from a thick and ancient lithosphere akin to the 
continental lithosphere on Earth. The terrestrial continental 
lithospheres are up to 200-300 km thick; their formation 
ages vary widely but generally reflect the ages of the subja- 
cent Archean and Proterozoic cratons. The lithospheres 
were stabilized against thermal convection, in part, by their 
compositional buoyancy resulting from the extraction of 
basaltic and komatitic magmas in the Archean. The litho- 
sphere is depleted in A1 2 0 3 and CaO but has elevated MgO 
and Mg* values. The lithosphere is variably metasomatized 
by kimberlitic and carbonatitic magmas so that incompati- 
ble trace elements are decoupled from the more refractory 
major elements. Because of these ancient and continuing 
metasomatic events, the isotopic systematics are widely 
varying. The e Nd values, for example, range from e Nd +42 to 
e Nd -55[6]l 

At first glance a “lithosphere model” would appear to 
explain the compositions of the SNC meteorites and the 
thickness of the Martian crust. A number of observations 
argue against this hypothesis. On earth, there is no Hadian 
(~4.5Ga) lithosphere as indeed there is no pervasive conti- 
nental crust that dates back to the formation of the Earth 
[6]. Second, the existing Archean crust probably never 
covered the entire planet [8]. But if we allow for a planet- 
wide Martian crust and complementary lithosphere, the 
energy associated with radioactive decay in the primitive 
astheno sphere would gradually increase the internal tem- 
perature of the planet until thermal buoyancy of the hot 
interior would cause a catastrophic delamination of ail or 
part of the crust [9]. Crustal -lithosphere foundering and 
widespread volcanic resurfacing of Mars would ensue. 
Perhaps this is the process that explains the northern hemi- 
sphere and the crustal dichotomy. The timescales of the 
process, however, should be over 100’s of millions of years 
[9]. The process, moreover, probably would not be a one- 
time event but would repeat in some cyclic pattern as pos- 
tulated for Venus [9]. While the “lithosphere” hypothesis 
cannot yet be totally discarded, the focus of the rest of this 
abstract will deal with the implications of a magma ocean 
on Mars. 

Fractionation of a magma ocean on Mars would have 
important similarities and differences with the Moon. On 
the Moon, the crystallization sequence of the magma ocean 
is dominated by olivine and orthopyroxene; the progressive 
enrichment of FeO over MgO in these phases makes the 
density of the cumulates increase upward with progressive 
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fractionation leading to gravitational instability [7]. The 
concomitant FeO enrichment in the liquid and its saturation 
in normative plagioclase component leads to the crystalli- 
zation and flotation of plagioclase and the stabilization of 
an anorthosite Lunar crust. 

In contrast, the liquidus phase in the Martian magma 
ocean is mainly gamet-majorite at depths greater than about 
1000 km [10,1 1]. For a core radius of about 1400 km, this 
region will extend over the lower 1000 km of the mantle. 
At these depths, majorite is predicted to be only slightly 
negatively buoyant [12]. It is likely that crystal-liquid frac- 
tionation will be inefficient in the strongly convecting 
magma ocean so that the lower mantle may initially ap- 
proximate the whole Mars composition. Once convection 
ceases the remaining intergranular liquid with the majority 
of the highly incompatible elements, including the heat 
producing ones, will migrate upwards to the still liquid 
magma ocean. The net effect is to sequester A1 2 0 3 (as gar- 
net-majorite) in the lower mantle. Gamet-majorite will 
transform to eclogite in the upper mantle so that it should 
remain as a stable layer in the lower mantle being less 
dense than only dense FeO-enriched sinking cumulates 
from below the crust. 

This scenario has several appealing features. The frac- 
tionation of majorite will decrease the Al 2 0 3 /Ca0 ratio in 
the remaining cumulates, a feature characteristic of Shergo- 
tite parent liquids [1]. The remaining liquids of the magma 
ocean will have less normative plagioclase resulting in 
delayed crystallization of plagioclase. The buildup of H 2 0 
in the residual magma ocean will also delay the onset of 
plagioclase crystallization. The net effect is that the anor- 
thosite crust on Mars should be relatively thinner than that 
on the Moon. The progressive build up of H 2 0 and other 
fluids in the residual liquids of the magma ocean may cause 
plagioclase to become negatively buoyant. The anorthosite 
crust would no longer float resulting in the emption of 
volatile-rich residual melts onto the Martian surface. The 
upper Martian crust, in this model, would be formed by 
products of the late stage magmas. The lower crust would 
be anorthosite-rich and would have acted as the platform on 
which subsequent magmatism was deposited. 

Besides preserving the isotopic and compositional het- 
erogeneity, overturning of an initially unstably stratified 
cumulate mantle may explain additional aspects of the evo- 
lution of Mars. After the overturn, the mantle would be 
stably stratified. Since thermal convection would then be 
inhibited by the compositional stratification, the primordial 
and radiogenically generated heat in the Martian interior 
would be slow to dissipate only by conduction. Similarly, 
the infertility of the cumulate mantle and the inhibition of 
thermal convection could limit magma production by adia- 
batic decompression, perhaps explaining in part the preser- 
vation of ancient crust. 

The crustal dichotomy may be the product of the over- 
turn of initially unstable magma oven cumulates. Horizon- 
tal convergence would thicken the crust above long wave- 


length (spherical harmonic degree one) downwelling of 
dense late stage cumulates. Thinning due to extension or 
delamination of the lithosphere in the opposing (northern) 
hemisphere could induce pressure-release melting and re- 
surfacing. A similar catastrophic resurfacing model has also 
been suggested for the Venus [9]. 

Regardless of the origin of the crustal dichotomy, 
variations in the thickness of the Martian crust are well [13] 
and it is at least plausible that these variations have per- 
sisted from the very early evolution. To preserve this 
crustal thickness variation over a substantial fraction of 
Mars’ evolution requires low temperatures at the base of 
the crust. Evolution models including heat transfer by 
thermal convection suggest established by gravity and to- 
pography data lower crustal temperatures high enough that 
crustal thickness variations may not be preserved [14]. The 
absence of thermal convective heat transfer due to stable 
compositional stratification would allow conductive cool- 
ing to result in lower crustal temperatures. 

Remnant magnetism [ 1 5] in a portion of the ancient, 
heavily cratered Martian crust argues for the existence of 
an internally generated magnetic field early in the evolu- 
tion. The generation of a magnetic field requires suffi- 
ciently rapid cooling of a molten metallic core. It has been 
argued that thermal convection beneath a cool, stagnant, 
conducting lithospheric lid would not cool the core rapidly 
enough, prompting the suggestion that the more rapid heat 
transfer associated with plate recycling may be required 
during the evolution of Mars [16]. Mantle overturn that 
transports relatively cool, dense cumulates to the mantle- 
core boundary may provide an alternative way to cool the 
core rapidly enough to sustain a magnetic dynamo. This 
cooling might induce a transient magnetic field for a part of 
early in Martian history, consistent with the observation 
that not all ancient cratered terrain is magnetized; stable 
mantle stratification after overturn could prevent subse- 
quent convective heat transfer and suppress continued 
magnetic field generation. 
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Jagoutz, E. (1991) Space Sci. Rev. 56, 13. [5] Borg, L.E. et 
al., (1997) GCA 61,4915. [6] Pearson, D. (1999) Lithos 48, 
171. [7] Hess, P.C. and Parmentier, E.M. (1995) EPSL, 
134, 501 [8] deWit, M. (1998) Pre Camb. Res. 91, 181. [9] 
Parmentier, E.M. and Hess, P.C. (1992) GRL 19, 2015. [10] 
Bertka, C.M. and Fei, Y. (1997) JGR 102, 5251. [11] 
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790. [16] Nimo, F., and Stevenson, D.J. (2000) JGR 105, 
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Introduction: Ten yeas ago we had 7 SNC 
meteorites while today we have 26 SNC meteorites. 
The newly found SNC's are essentially duplicates of 
the previously ones with some slight variations in their 
chemical composition. It seems that there are only a 
limited number of lithologies present on the SNC 
parent body, and we have already sampled the 
prominent lithologies. 

For many years the “ages” of SNC meteorites have 
been under discussion. The 150 Ma to 500 Ma mineral 
ages for shergotites will be called the young or 180 Ma 
event, and will be the main focus h this contribution. 
Nyquist argued that the 1 80 Ma event might represent 
the age of crystallization. Indeed we have good 
evidence that some mineral equilibration took place at 
this time since in some cases Rb-Sr and Sm-Nd 
systematics give the same age. However, the 
conclusion that shergotites are magmatic rocks that 
were erupted at this time is not substantiated by data. 
To the contrary there is good evidence that a chemical 
fractionation, which is typical for magmatic rocks (Rb- 
Sr fractionation), was not associated with the young 
disturbance. Instead we find chemical fractionation 
patterns (Sm-Nd fractionation) that are not observed in 
any basalts. 

Rb-Sr systematics: 

We have good insight about the geochemical 
changes in the Rb/Sr during partial melting and 
fractional crystallization. If a reservoir is partially 
melted, it will be depleted in Rb compared to Sr, while 
the partial melt will be enriched in Rb compared to Sr. 
The SNC's, however, plot on the 4.5 Ga reference line 
in a Rb-Sr isochron diagram, implying that the Rb/Sr 
was unchanged for 4.5 Ga. 

Partitioning ofRb in pyroxene: 

Recently in a significant contribution of L. Borg et 
al (2002) new data on Alha77005 and Lev are 
reported. These authors claim that both meteorites 
have the same age of 175 Ma. Shih et al (1982) and 
Jagoutz (1989) also reported isotope data on Alha, and 
they found age's of 180 Ma and 165 Ma, respectively. 
It is interesting to compare the results of these three 
studies. From an in-depth investigation of mafic 
minerals we know that pyroxene does not 
accommodate Rb in the crystal structure. Clean 
mineral separates of pyroxene revealed that the 
concentration of Rb in pyroxene might be less than 
lppb (Zindler and Jagoutz 1986). Even ultra -high 
pressure pyroxene containing up to 1% of K20 does 


not have more Rb (Owen, unpublished data on 
diamond inclusions). Therefore, we use Rb 
concentration as a parameter of how clean a pyroxene 
mineral separate is. 

A Sr versus Rb plot has been made for all 
published mafic minerals and whole rocks from Alha 
77005. The resultant linear array is likely a mixing line 
caused by impact because impact glasses commonly 
have about whole rock chemical composition. If true, 
all the pyroxenes in these three studies would be 
dominated by impact glasses in their Rb-Sr systematic. 
The main difference is that Nyquist did not wash the 
pyroxene at all, while Borg leached the pyroxene using 
2N HC1, and Jagoutz additionally removed shock 
glasses from the pyroxenes using 5% HF. Because the 
pyroxene shows a fine network of cleavages, it is 
unlikely that all the shock glasses was removed by 
leaching. Although considerable time was spend on the 
ion-probe to measure the Rb, efficiency of one count 
per second per ppm cannot measure ppb level 
concentrations of Rb even if all the interferences on the 
Rb masses are filtered out. Furthermore it is hardly 
possible to find areas of the pyroxene without 
cleavages in order to measure the intrinsic Rb 
concentration. 

This linear array does not go thru zero and it seems 
there is a bias of Rb in this minerals. A similar pattern 
is also found in other (shocked) shergotites 
experiencing the young disturbance, but not in the 
much less shocked nakhlites where the array intersects 
at the intrinsic Sr concentration of pyroxene. It seems 
that there are two components of contamination. One 
component is the shock glass, which is mobilized and 
contains Rb and Sr in the proportions of the whole 
rock. The other component does not contain Sr and 
might be carried by a shock-related vapor phase 
containing only highly volatile elements like Rb. 

Age significance of the Rb/Sr data: 

We can assume that prior to the shock event the 
shergotites are well behaved magmatic rocks. Here 
plagioclase, the main host phase of Sr, and pyroxene 
both plot close to the Sr initial isotopic composition as 
we observe in Nakhla. Olivine might contain neither 
Rb nor Sr in appreciable amounts. The low alkalis in 
these rocks might not form an alkali rich mineral. The 
whole rock might have a more radiogenic Sr than any 
of the major phases, requiring that it was produced by 
interstitial Rb or the Rb in minor phases. This 
radiogenic Sr might then easily be equilibrated with 
the plagioclase during a shock event. The observation 
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explaining that the plagioclase and the phosphate have 
different Sr isotopic composition while neither contain 
much Rb. 

The shock essentially mobilizes Rb and produce 
conditions for dating the shock at least in a first order 
approach. We cannot expect perfect isochrones under 
these conditions, but there could be some age 
significance in these “scatter-chrones”. All the Rb-Sr 
systematics do provide a good indication that the “180 
Manage event may possibly represent a shock event. 

Sm-Nd systematics: 

Nd isotopes are much more complicated than the 
Sr isotopes. There are three points, which have to be 
considered in a genetic model. 

1) The Sm-Nd fractionation of the whole rocks is 
very strong in shergotites, especially relative to the 
minor Nd-La and Sm*Yb fractionation. Even Sr-Nd 
and Lu-Hf are fractionated less than Sm-Nd. It seems 
there is a “step-function” in the REE, which cannot be 
explained by a simple magmatic process. 

2) The high Sm/Nd in the pyroxenes and their 
relatively unradigenic ,43 Nd/ ,44 Nd may represent the 
same problem as observed for Sr, and be caused by 
incorporated shock-glass. However, even if all shock- 
glass is eliminated, these pyroxenes still have 
unradiogenic Nd and high Sm/Nd, indicating that the 
depletion of the LREE was caused by the young event. 

3) Terrestrial weathering is a serious obstacle 
especially in LREE depleted shergotites. In some cases 
(Dhofar 019) even clean pyroxene might have been 
contaminated by terrestrial weathering. While we 
appear able to remove terrestrial contamination from 
pyroxene by acid leaching, we cannot use the same 
technique on whitlockites. In fact whitlockites are not 
directly measured but rather it is the acid leaches in 
which the whitlockites are dissolved that are taken as 
whitlockites. Possibly at the same time all interstitial 
terrestrial contamination dissolves in the same leach. 
This is unfortunate because the whitlocktes contain up 
to 90% of the Nd in shergotites. Therefore the ages 
and initial Nd isotopic composition obtained by the 
whitlockite - pyroxene isochrones are highly 
questionable. Leaches appear to form a mixing line 
with terrestrial sediments on an isochron plot. 

U-Pb Systematics: 

In SNC meteorites the U resides almost 
quantitatively in the whitlockites. Maskelinites, the 
host phase for Pb, has a very low U/Pb. We performed 
experiments in order to find the intrinsic Pb from the 
maskelinites and avoid terrestrial contamination. The 
Pb isotopic composition measured in clean 
maskelinites lines up nicely close to the Geochron. 


Presumably this implies that the SNC meteorites were 
formed by a magmatic process at 4.55 Ga, and 
remained undisturbed since then up to the 180 Ma 
event. This young disturbance increased the U/Pb 
drastically, which is in contrast to other meteorites, 
which have had their U/Pb decreased recently. 

Clean pyroxene separates from Nakhla have a low 
U/Pb and, therefore, their Pb isotopic composition is 
similar to the plagioclase. Pyroxene extracted from 
LA, however, have U/Pb of >100 and we must assume 
that the U/Pb was changed by the shock event. 
Terrestrial magmatic pyroxenes have a low U/Pb while 
pyroxenes of chondrites and CAI inclusions have a 
very high U/Pb. Also pyroxenes from angrites have a 
very high U/Pb. Some pyroxene from HED meteorites 
have a high U/Pb and some have a low U/Pb. Our 
current explanation for this observation is that CAI 
inclusions and chondrites condensed from a hot solar 
nebulus; additionally, angrites may have also 
condensed from a hot solar nebulus. In contrast, HED 
meteorites and SNC meteorites are of magmatic origin 
and, therefore, have a low U/Pb in their pyroxenes. 
Some of these meteorites have increased the U/Pb of 
their pyroxene by a superimposed later shock event. 



Ref: L. Borg et al (2002) in press. Shih et al 
(1982) GCA 46, 2323. Jagoutz (1989) GCA 53, 2429. 
Zindler and Jagoutz (1986)GCA 52, 319. 
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Ancient History. Jones [1] first suggested that the 
inverse covariation of initial e( ,43 Nd) and 87 Sr/ 86 Sr of the 
shergottites could be explained by interaction between 
mantle -derived magmas with another isotopic reser- 
voirs) (i.e., assimilation or contamination). In that 
model, magmas were generated in a source region that 
was isotopically similar to the Nakhla source and the 
second reservoir(s) was presumed to be crust. The text 
of [1] also permitted the second reservoir to be another 
type of mantle, but I can confirm that a second mantle 
reservoir was never seriously considered by that au- 
thor. Other features of this model were that (i) it oc- 
curred at a particular time, 180 m.y. ago, and (ii) the in- 
teracting reservoirs had been separated at ~4.5 b.y. In 
a later paper [2], Jones explored this mixing model more 
quantitatively and concluded that magmas from a Nak- 
hla-like source region at 1 80 m.y. would fall on or near 
an isotopic Nd-Sr-Pb hyperplane defined by the sher- 
gottites. This criterion was a necessary prerequisite for 
the parent magma(s) of the shergottites to have initially 
been Nakhla -like isotopically. 

At this juncture, it is perhaps worthwhile to note 
that this mixing model was not presented to explain 
geochemical variations but as a justification for a 1 80 
m.y. crystallization age for the shergottites and a 1 .3 
b.y. crystallization age for the nakhlites [1,2]. In the 
mid-1980’s crystallization ages estimated for Nakhla 
ranged from ~1.3 b.y [3] to 4.5 b.y. [4]. Similarly, pre- 
ferred crystallization ages for the shergottites ranged 
from 360 m.y. [5], to 1.3 b.y. [6], to 4.5 b.y. [7]. In all 
these models, the 180 m.y. event seen in the shergot- 
tites was deemed to be metamorphic [e.g., 6]. The fit 
between the Nakhla -like source region and the shergot- 
tite hyperplane was a validation both of the 1.3 b.y. 
igneous age of Nakhla and the 180 m.y. igneous age of 
the shergottites. 

Less Ancient History. To some degree the original 
mixing model survives today [e.g., 8], but many 
changes have been required in the intervening decade. 
The original model envisioned a homogeneous mantle 
with a e( l43 Nd) of ~+20 and 87 Sr/ 86 Sr of -0.7036 at 180 
m.y. The crustal contaminant had radiogenic Sr 
( 87 Sr/ 86 Sr > 0.722) and e( ,43 Nd) of - -20 [9], We now 
know, for example, that the mantle cannot be homoge- 
neous because of l42 Nd considerations [10]. 

Insights from I42 Nd. 146 Sm decays to 142 Nd with a 
half life of 1 03 m.y. and thus records early planetary 
events. Borg et al. [10] have shown that, although the 
shergottites appear to form a mixing line on a e( 142 Nd) 
vs. e( I43 Nd) diagram, in agreement with the original 


model, the nakhlites do not fall on that trend. So, ironi- 
cally, one of the main premises of the orginal model 
cannot be true. Shergottite parent magmas must come 
from a distinct source region from that of the nakhlites. 
This, in turn, implies that the martian mantle is not ho- 
mogeneous, but the scale of this heterogeneity is un- 
known. Another implication is that the nakhlite source 
region was re-fertilized (i.e., Sm/Nd decreased) early in 
martian history. Since the nakhlites have a very signifi- 
cant amount of radiogenic 142 Nd, this re -fertilization 
must be done carefully (or else the anomaly would have 
been diminished). 

Insights from QUE94201. Another shock to the 
original model came with the finding of QUE94201. At 
180 m.y. the QUE source region w>uld have had an 
e( 143 Nd) of ~+50 and 87 Sr/ 86 Sr of ~0.7014, significantly 
more depleted than the nakhlite source [11]. This raises 
the issue of how depleted a martian mantle can be and 
still produce basaltic magma at 1 80 m.y. 

Insights from Hf-Nd isotopes . The Lu-Hf and Sm- 
Nd chronometer systems both indicate that large ele- 
mental fractionations occurred during SNC petrogene- 
sis [11,12]. These are indicative of complex igneous 
processes that occurred during or immediately preced- 
ing the melting events that produced the SNC’s (e.g., 
QUE94201 [11]). Even so, initial SNC Hf and Nd iso- 
topic ratios, allow reconstruction of the Sm/Nd, Lu/Hf, 
and Ho/Tm ratios of the SNC source regions [13]. 
Comparison of lunar and martian Lu/Hf vs. Sm/Nd sys- 
tematics show strong similarities [14,15]. The 
QUE94210 source appears similar to that of high-Ti 
lunar basalts; and the martian “crust” inferred from 
shergottites appears very similar to KREEP. 

How these observations should be interpreted is 
not clear, but certainly, if the Moon passed through a 
magma ocean stage, then perhaps Mars did as well. 
There are other similarities between lunar and martian 
basalts. Both appear to come from depleted source 
regions (e.g., super-chondritic Sm/Nd and low AI2O3). 
The above discussion indicated that the martian mantle 
has not been homogenized and the same appears true 
for the lunar mantle. 

Insights from experiments . Martian magmas tend to 
be hot by terrestrial standards. The liquidus tempera- 
ture of the groundmass of EET79001 A, with an Mg# of 
~54, is ~1290°C [16]. The highly evolved composition 
QUE94201, with Mg# of ~38, has a liquidus temperature 
of~1170°C [17]. For comparison, the Kilauea tholeiitic 
lava that formed Makaopuhi Lava Lake, with an Mg# of 
~57 had a liquidus temperature of ~1190°C [18]. These 
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elevated liquidus temperatures are probably a conse- 
quence of the depleted nature of the martian mantle. 

Insights from oxygen fugacity. Two recent papers 
have emphasized the correlation of isotopic (Nd, Sr) 
and chemical (La/Yb) parameters with oxygen fugacity 
within the shergottite suite [19,20]. The greater the 
amount of “assimilant” or “contaminant” the higher the 
oxygen fugacity recorded by Eu anomalies and Fe-Ti 
oxide assemblages. Therefore, there is an inference 
that the assimilated (or contaminating) component was 
more oxidized than pristine basalts emanating from the 
martian mantle. This inference, coupled with the lack of 
observational evidence for subduction zones on Mars 
[21], suggests that the martian mantle is reduced (~IW- 
1) and that assimilation of crustal materials oxidizes the 
assimilating magmas to varying degrees. Because 
Shergotty and Zagami experienced the most “contami- 
nation,” it was in fact an accident that early studies of 
the shergottites pronounced them “Earthlike” in terms 
of their oxygen fugacity [22]. 

This new presumption that the martian mantle is ac- 
tually rather reduced makes sense in terms of phase 
equilibria. The FeO content of SNC basalts is not that 
different from typical lunar basalts and eucrites (~18-19 
wt.%). As Stolper [23] has shown, the fo 2 at which 
such basalts are in equilibrium with iron metal is ~IW-1, 
in agreement with the fo 2 inferred for the most reduced 
SNC, QUE94201 [19]. Apparently, the martian mantle 
has remained relatively isolated over geologic time; and 
it maintains the redox state it had at the time of core 
formation, ~4.5 b.y. ago. 

Discussion: From these insights an incomplete pic- 
ture emerges that may nevertheless serve as a useful 
guide to the igneous evolution of Mars. First, if the 
Moon had an early magma ocean, it is likely that Mars 
did too. The contaminant of the shergottites is likely 
the solidified dregs of that magma ocean — martian 
KREEP. The physical location of martian KREEP is 
likely to be one of three regions: (i) the crust; (ii) the 
crust-mantle interface; or (iii) deep in the mantle (due to 
convective overturn). But the correlation between de- 
gree of contamination and redox state argues more for a 
crustal contaminant than for mixing between mantle 
reservoirs. 

Second, it is highly unlikely that the source region 
of QUE94201 could represent typical martian mantle. 
This basalt contains so little in the way of heat- 
producing elements that recent (< 1 80 m.y.) volcanism 
on Mars would be highly improbable, if not impossible 
[24]. A more likely scenario is that QUE mantle was 
melted by an external heat source, perhaps at a thermal 
boundary layer. Because QUE mantle is so depleted, it 
is likely to be buoyant relative to more fertile mantle. 
This conjures an image of a layered martian mantle with 


fertile material below and convectively-isolated, less- 
fertile material above. Conceivably, this layering was 
established during or after the solidification of the 
magma ocean. 

Third, the more complex history of the nakhlite 
source region is most easily achieved by an internal 
differentiation to avoid resetting 542 Nd [10]. A small 
degree partial melt, moving from one portion of the 
nakhlite source to another could lower the Sm/Nd ratio 
of the latter at the expense of the former without chang- 
ing the 142 Nd/ 144 Nd of either. This differentiation pre- 
sumably occurred fairly late (~4 b.y.), after the decay of 
146 Sm. In principle, this differentiation event could also 
have resulted in additional layering. 

Summary: The martian mantle is apparently het- 
erogeneous, which opens the possibility that it is lay- 
ered, with each layer convectively isolated. If this is 
correct, melt generation should occur either at thermal 
boundary layers or in plumes generated at those 
boundaries. Mantle layering may be a good means of 
slowing the planet’s cooling rate, allowing young vol- 
canism. Layering may also provide a means for keeping 
the crust and upper mantle cool, allowing the preserva- 
tion of ancient variations in crustal thickness [25]. 

References: [1] Jones J.H. (1986) GCA, 50, 969-977. 
[2] Jones J.H. (1989) Proc. Lunar Planet. Sci. Conf 
19th, 465-474. [3] Gale N.H. et al. (1975) EPSL 26, 195- 
206. [4] Papanastasiou D.A. and Wasserburg G.J. (1974) 
GRL 1, 23-26. [5] Jagoutz E. and Wanke H. (1986) GCA 
50, 939-953. [6] Shih C.-Y. et al. (1982) GCA 46, 2323- 
2344. [7] Chen J. H. and Wasserburg G.J. (1986) GCA 
50, 955-968. [8] Borg L.E. et al. (2002) GCA 66, 2037- 
2053. [9] Longhi J. (1991) Proc. Lunar Planet. Sci. Conf. 
21 st . 695-709. [10] Borg et al. (2001) Lunar Planet. Sci. 
XXXII. #1144. [11] Borg et al. (1997) GCA 61,4915-4931. 
[12] Blichert-Toft J. et al. (1999) EPSL 173, 25-39. [13] 
Bostok D. (1973) Trans. Phil. Soc. Agua Blanca 97, 892- 
897. [14] Jones J.H. (2000) Lunar Planet. Sci. XXXI. 
#1474. [15] Beard B.L. et al. (2002) Lunar Planet. Sci. 
XXXIII. #1933. [16] Wasylenki L.E. et al. (1993) Lunar 
Planet. Sci. XXIV, 1491-1492. [17] McKay G. et al. (2002) 
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Herd C.D.K. et al. (2002) GCA 66, 2025-2036. [21] Phil- 
lips R.J. and Ivina E.R. (1979) Phys. Earth Planet. Inter. 
19, 107-148. [22] Stolper E.M. and McSween H.Y. Jr. 
(1979) GCA 43, 1475-1498. [23] Stolper E.M. (1977) GCA 
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The shergottites are basaltic meteorites from Mars, many 
of which have igneous crystallization ages of about 1 80 million 
years ago [1,2]. The observed density of small impact craters in 
MOC images of Olympus Mons and Elysium Planitia indicates 
that some volcanic flows in these regions are just 10 to 30 
million years old [3]. The existence of young volcanism on 
Mars implies that adiabatic decompression melting and hence 
mantle convection has remained an important process on Mars 
to the present day. Numerical simulations of mantle convection 
and magma production set important constraints on the current 
properties of the mantle of Mars, particularly the minimum 
abundance of radioactive elements that has been retained in 
the mantle. 

Mantle Convection Model 

The mantle convection simulations are performed using 
finite element methods in spherical axisymmetric geometry 
[4,5]. The best observational constraint on the pattern of den- 
sity heterogeneities in the martian mantle is the non-hydrostatic 
geoid, which shows a high degree of axi symmetry about central 
Tharsis [6] . This indicates that the spherical axisymmetric flow 
model is a very good approximation for modeling present-day 
convection on Mars. The model includes both internal heating 
from mantle radioactivity and basal heating due to the heat 
flux out of the core. The internal heating is related to the broad 
Tharsis plateau, and the basal heating produces mantle plumes 
that feed individual shield volcanos such as Olympus Mons. 

The finite element grid has a vertical resolution of 13 km 
and an assumed core radius of 1700 km. The model includes 
a 210 km thick high viscosity lid at the top of the model, with 
a viscosity contrast of 10 5 between the surface and the base of 
the lid. The characteristic flow velocity at the top of the high 
viscosity layer is just 0.01 mm/year, confirming that the model 
is in the stagnant lid convective regime [7]. In the preferred 
model, the heat flux is 1 6 mWm~ 2 at the base of the crust and 
22 mWm~ 2 at the surface. This is consistent with flexure 
studies using Mars Global Surveyor gravity and topography 
data, which indicate that the near-surface heat flux at large, 
Amazonian age volcanos is 16-28 mWm~ 2 [8]. 

Radioactivity Models 

Models for the silicate composition of Mars are con- 
strained by the chemistry of the SNC meteorites and imply 
present-day radioactive heating rates of 4. 1 -6.2 x 10" 12 W kg~ l 
[9,10]. In the convection simulations, 30 to 90% of the total 
radioactivity is partioned into the crust. The crust is assumed 
to be 50 km thick [11], with radioactivity uniformly distributed 
within the crust. The remaining radioactive heating is assumed 
to be uniformly distributed throughout the mantle. Models 
that retain large amounts of radioactivity within the mantle 
have higher present-day mantle temperatures and thus pro- 
duce greater amounts of magmatic activity. An important goal 
of this study is to determine the minimum amount of radioac- 
tivity that must have been retained in the mantle in order to 


permit pressure release melting to occur on present-day Mars. 

Magma Production Calculations 

Melting relationships are taken from the experimental 
work of Bertka and colleagues on a Mars-analog composi- 
tion [12, 13]. These experiments cover the pressure range 
1-3 GPa and 12-25 GPa. In my calculations, most melting 
occurs near a pressure of 5 GPa, so some interpolation of the 
solidus location is necessary. Recent work by Draper et al. 
[14] on a slightly different composition covers the pressure 
range 5-9 GPa. Melting calculations using this solidus will 
be presented at the conference. I have assumed a dry solidus 
based on studies indicating a low oxygen fugacity in the man- 
tle [15-17]. Existing experimental results are for a primitive 
martian mantle composition. Because the present-day mantle 
has had its lowest melting point components removed via pre- 
vious melting, the actual solidus should be somewhat higher 
than indicated by these experiments. A 50 km thick crust [11] 
constitutes only about 5% of the mantle volume, so this ef- 
fect is likely to be small, but experimental melting studies of 
an evolved mantle composition would be helpful. Melt pro- 
ductivity above the solidus is assumed to be 0.4% per degree 
based on studies of KLB-1 [18]. Melt volumes are determined 
by tracing streamlines through the melting region [19]. Latent 
heat of melting and specific heat are taken from the compilation 
of Navrotsky [20]. 

Results 

In these models, magma production is restricted to mantle 
plume heads, which is consistent with the localized distribution 
of young volcanism on Mars (Figure 1). In order for the man- 
tle to remain warm enough to permit on-going pressure-release 
melting, the minimum average radioactive heating rate in the 
mantle must be at least 1.6 x 10" 12 Wkg~ x . This corresponds 
to 40% of the total radioactivity in the Wanke and Dreibus 
model [9] or 25% of the total radioactivity in the Lodders 
and Fegley model [10]. The corresponding average mantle 
temperature and viscosity corresponds to a thermal Rayleigh 
number for the martian mantle that exceeds 10 6 at present. 
This implies that mantle convection on Mars remains relatively 
vigorous at present. As a corollary, the density anomalies and 
dynamic topography associated with this convective flow must 
be considered in any complete model for the long-wavelength 
geoid and topography of Mars. Surface volcanic loads obvi- 
ously make an important contribution to the long- wavelength 
geoid [e.g, 21], but models that consider only surface loads are 
incomplete descriptions of the planet. 

The best current convection simulation has a present-day, 
average magma production rate of 7 x 10~ 6 km 3 y ear ~ x . As- 
suming a mean lava flow thickness of 5 meters and that 10% 
of total magma production is extruded on the surface, this cor- 
responds to resurfacing 140 km 2 per million years. Transient 
periods with resurfacing rates a factor of 1 0 larger than this 
sometimes occur in the models. For comparison, the Hawaiian 
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Figure 1: A simulation of mantle convection on Mars. The color scale shows potential temperature, which is the physical temper- 
ature with the adiabat removed. The white region in the mantle plume head is above the solidus. 


hotspot on Earth has a long-term average magma production 
rate of about 0.01 km 3 year" 1 [22]. The mean melt fraction in 
the numerical models ranges up to 8%, in very good agreement 
with estimates of 2-8% partial melting in the shergottites based 
on rare earth element concentrations [23]. 

The melt production rate is extremely sensitive to the pre- 
cise value of the assumed solidus. Raising or lowering the 
solidus by 50 K can change the melt production rate by a fac- 
tor of 5. The high pressure solidus measured by Draper et al. 
[14] is about 50 K hotter than the solidus that I have interpo- 
lated from the results of Bertka and collaborators [12, 13]. This 
new constraint on the high pressure solidus implies that the re- 
quired amount of radioactivity retained in the mantle of Mars is 
somewhat higher than the value cited above. Detailed calcula- 
tions of this are presently underway and will be reported at the 
conference. Alternatively, if a wet solidus is appropriate, the 
required abundance of mantle radioactivity would be reduced. 
The mantle radioactivity value derived here is a mantle-wide 
average. Some of the shergottites (notably QUE94201 and 
SaU005) have very low radioactivity abundances and imply 
source regions that are highly depleted in radioactivity. The 
results of this study indicate that the average martian mantle 
must have a much higher present-day concentration of ra- 
dioactive elements. One possibility is that the martian mantle 
consists of blobs of both depleted and primitive compositions. 
At least in some cases, the low radioactivity blobs must occur 
in close association with high radioactivity blobs to account 


for the formation of meteorites such as QUE94201. 
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Introduction: Shergottites from Mars, like basalts 
on Earth, are important because they provide an oppor- 
tunity to determine the chemical and thermal evolution 
of the planet’s mantle, relationships between the crust 
and mantle, and the role of crustal assimilation in the 
petrogenetic processes that produce surface lavas. At- 
tempts to utilize these rocks, however, have been frus- 
trated because they are cumulate rocks rather than bulk 
melt samples that lie on liquid lines of descent. To 
overcome this problem, estimates of bulk melt compo- 
sitions have been made by trying to identify and sub- 
tract the compositions of cumulate minerals [1-5], ana- 
lyzing trapped melt inclusions [6-9], interpreting phase 
equilibria [10], and using minor element compositions 
in pyroxene with equilibrium coefficients [11]. Unfor- 
tunately, these attempts have produced contradictory 
results. 

QUE 94201: When Martian meteorite QUE 94201 
was first described, it appeared to be a bulk melt rather 
than a cumulate fraction [12-14], thus providing our 
first unambiguous melt composition for modeling pur- 
poses. The evidence seemed compelling. Typically, 
basaltic and lherzolitic shergottites contain cumulus 
olivine and/or pyroxene, estimated to range from 14% 
to 63% of the rocks [1,5]. In contast, there is no evi- 
dence of any xenocrystic or cumulus minerals in QUE 
94201. Indeed, the mineral assemblage, mineral com- 
positions, textures, trace elements in pyroxene, and the 
bulk composition of QUE 94201, suggested it fraction- 
ally crystallized in a closed system [12-15]. QUE 
94201 contains numerous melt inclusions and abundant 
interstitial melt pockets, suggesting that most intersti- 
tial melt was trapped in the rock rather than being re- 
moved during the crystallization process. Likewise, 
and in contrast to other shergottites, the rock contains 
an unusually large amount of accessory minerals, again 
suggesting that lower temperature melt fractions were 
not separated from the higher temperature crystalliza- 
tion products. This can be illustrated quantitatively by 
comparing the compositions of pyroxene cores (the 
highest temperature, and thus earliest, crystalline phase 
in the rock) with the calculated pyroxene composition 
that would be in equilibrium with the bulk rock (i.e., 
when it was initially wholly molten). Figure 1 shows 
the molar Mg/(Mg+Fe) across a zoned pyroxene crys- 
tal in QUE 94201. Also shown is the calculated py- 
roxene composition that would be in equilibrium with a 
melt of the bulk composition of the meteorite. Because 
the core pyroxene composition is identical to the equi- 


librium composition, it is clear the pyroxene is not a 
cumulus mineral and that a fraction of the melt was not 
lost after pyroxene began to crystallize. In contrast, 
similar agreement is completely lacking in other sher- 
gottites, as illustrated for the case of Elephant Moraine 
(EET) 79001 (Fig. 1, right panel). Consequently, the 
bulk composition of QUE 94201 seems to represent a 
bulk melt rather than a cumulate fraction. The minera- 
logical evidence is similarly expressed in the chemistry 
of the rock. QUE 94201 is more aluminous than any 
other shergottite, all of which have compositions that 
are skewed from melt compositions by Al-poor cumu- 
lus phases. In addition, QUE 94201 has higher abun- 
dances of many other elements which are incompatible 
during the earliest stages of crystallization (Ti, P, and 
REE). 



Fig. 1. Molar Mg/(Mg+Fe) in a linescan across a zoned 
pyroxene grain in QUE 94201 . The core composition is 
similar to the calculated composition for a pyroxene in equi- 
librium with a melt with the bulk composition of QUE 
94201. In contrast , core compositions of pyroxene in EET 
79001 lithology B are not similar to the calculated composi- 
tion for pyroxene in equilibrium with a melt with the bulk 
composition of that rock. 

Petrogenetic Contradictions: If QUE 94201 is a 
bulk melt, then it can be used to investigate the petro- 
genetic evolution of magmas on Mars using liquidus 
phase diagrams. When one does so, however, it is 
clear that something is amiss. To illustrate the prob- 
lem, Fig. 2 shows the projection of the bulk composi- 
tion of QUE 94201 into an Olivine-Plagioclase-Quartz 
ternary in the psuedoquartemary Olivine-Plagioclase- 
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Wollastonite-Quartz system often used to investigate 
basaltic melts [16]. Mineralogically, QUE 94201 is 
dominated by plagioclase (now maskelynite) and py- 
roxene [e.g., 12]. Plagioclase is zoned from An 57 to 
An 68 . Pigeonite is also continuously zoned, from 
Wo 9 En 6 oFs 31 to Wo^En^s^. Augite compositions are 
bimodal, averaging Wo 3 oEn 45 Fs 25 and Wo 40 En 2 Fs 5 8 . 
This is coherent set of compositions that suggests 
closed system fractional crystallization. Consequently, 
one would expect the bulk composition of QUE 94201 
to plot close to the pigeonite-augite-plagioclase bound- 
ary. Instead, the bulk composition plots in the olivine 
stability field, even though olivine does not appear to 
be an early crystallizing principal phase in QUE 94201. 
(The meteorite contains traces of fayalite, Fa 97 . 9 9 , 7 [12], 
but it was produced from the last vestiges of melt and 
is not relevant to the projected phase composition.) 
Consequently, while the phase diagram suggests the 
rock contains olivine, it does not. 


Qtz 



Fig. 2. Composition of QUE 94201 projected into a 
quartz-olivine-plagioclase liquidus ternary phase diagram in 
the psuedoquartemary (Qu artz- Olivine- PI agiocl as e- 

Wollastonite) system (after [15] and modifications after J. 
Longhi, personal comm., 1996]. Compositions Ql, Q2, Q4, 
and Q5 represent the uncertainties among different analyses 
of the QUE 94201 ([12] and unpublished data; [17]). Stan- 
dard abbreviations are used: Ol=olivine, Pl=plagioclase, 
Qtz-quartz, sil~silica, aug-augite, lpyx~low-Ca pyroxene. 

A Porphyritic Olivine Basalt?: There are two 

ways to resolve this conundrum. The chemical compo- 
sition of QUE 94201 used in the projection could be in 
error. If the rock contained several more weight per- 
cent Si0 2 , then the composition would not plot in the 


olivine stability field. The data seem robust, however, 
having been replicated with multiple splits in the same 
laboratory and in different laboratories. 

The conundrum can also be resolved if the meteor- 
ite is not representative of its parent lithology. While 
QUE 94201 is in chemical equilibrium with olivine, it 
may have been physically separated from olivine. It is 
perhaps worth remembering that QUE 94201 is a very 
tiny sample, weighing only 12 grams, for a medium- 
grained (1 to 5 mm) rock. It is conceivable that it is 
not representative of the outcrop or near surface lithol- 
ogy from which it was excavated from Mars. In par- 
ticular, it could potentially be the matrix of a porphy- 
ritic olivine basalt, rather than a representative bulk 
sample of the basalt. This accounts for olivine being 
part of the residuum in the phase space topology, but 
not actually part of the meteorite. While other sher- 
gottites are plagued with too many crystals (xenocrysts 
and cumulate phases), QUE 94201 may be plagued 
with having too few (the missing olivine). Unfortu- 
nately, this problem suggests QUE 94201 is not 
representative of a true bulk melt, nor can a true bulk 
melt be calculated easily because the proportion of 
olivine in the bulk basalt from which it came is un- 
known. If olivine was physically separated from QUE 
94201, the meteorite’s bulk composition will not lie on 
an equilibrium liquid line of descent. Rather, it will lie 
on an open system crystal fractionation line that is 
controlled by the loss of olivine. 

References: [1] Stolper E. and McSween H.Y. Jr. 
(1979) GCA, 43, 1475-1498. [2] Smith J.V. and Hervig 
R.L. (1979) Meteoritics, 14, 121-142. [3] Treiman 
A.H. (1986) GCA, 50, 1061-1070. [4] McSween H.Y. 
Jr. et al. (1988) EPS XIX, 766-767. [5] McCoy T.J. et 
al. (1992) GCA, 56, 3571-3582. [6] Johnson M.C. et 
al. (1991) GCA, 55, 349-366. [7] Harvey R.P. and 
McSween H.Y. Jr. (1992) GCA, 56, 1655-1663. [8] 

Harvey R.P. et al. (1993) GCA, 57, 4769-4783. [9] 
Treiman A.H. (1993) GCA, 57, 4753-4767. [10] 
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464. [11] Treiman A.H. and Sutton S.R. (1992) GCA, 
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Detailed chemical and isotopic measure- 
ments reveal complexities in planetary dif- 
ferentiation that will make any reconstruc- 
tion of the SNC source region difficult. Re- 
cent age determinations of shergottites have 
resolved a long-standing uncertainty in their 
crystallization ages in favor of ages in the 
range of 200 to 400 Ma [1], These young 
ages, together with 1.3 Ga ages of the 
Nahklites and Chassigny, require that the 
SNC source regions maintained more pro- 
nounced depletions of incompatible ele- 
ments than the terrestrial MORB source 
since primordial differentiation (4.525 Ga) 
and that the parent magmas of the classical 
shergottites (Shergotty, Zagami) assimilated 
a long term light-REE enriched component 
(crust) [2], The isotopic evidence indicates 
little or no mixing between differentiated 
and primitive mantle or between differenti- 
ated mantle and crust during martian history. 

Estimated SNC parent magma composi- 
tions [e.g., 3] have lower A1 2 0 3 concentra- 
tions than terrestrial MORB and OIB [2], 
These low AI 2 O 3 contents are consistent 
with highly depleted source regions as in- 
ferred from the isotopic data. New modeling 
of polybaric melting reveals that the extent 
of the depletion in terms of major elements 
is greater than expected for basalt extraction. 
There are two general scenarios conistent 
with the major elements: one is that the por- 
tion of the martian mantle from which the 
SNC magmas are derived is a magma ocean 
cumulate; the other is that initial melts in a 
polbaric fractional fusion event were never 
aggregated with later melts (possibly the 
inial melts were too dense). Loss of the ini- 
tial melts is also consistent with the ex- 
tremely high Sm/Nd ratios inferred for the 


parent magmas of the Antarctic shergottites. 
Figure 1 contrasts the results of batch and 
fractional fusion calculations for a martian 
mantle (dDW — open diamond) that is de- 
pleted relative to primitive Mars (DW [4] — 
open square) in a manner analogous to the 
MORB source (dPUM) is related to primi- 
tive upper mantle (cross — PUM [5]). source 
compositions with the compositions of the 
SNC magmas. The string of blue circles 
represents 10% batch melts of dDW from 
35 kb to 6 kb. The string of open suares are 
polybaric fractional melts of the dDW 
source also beginning at 35 kb. The strings 
of triangles represents the same melting be- 
ginning at 35 kb, but from which the 35 to 
26 kb melts have been removed. The string 
of cicles is analogous to the triangles, but in 
tis caseThe AI 2 O 3 content of the source was 
arbitarily reduced by 50%. These results 
suggest a way to have source regions with 
apparently low A1 without a magma ocean. 
However, at least two sources appear neces- 
sary to yield the nahklite and shergottite 
parent magmas. 
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Figure 1. Olivine projection in the system lamite (La)-olivine(01)-nepheline + CaAl 2 0 4 (NeCA)- 
quartz (Qtz). Liquidus boundaries and field of SNC parent magma compositions modified after [2,3]. 
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Introduction: Oxygen fugacity (t0 2 ) is a key physi- 
cal parameter which influences the crystallization se- 
quences of magmas, as well as the composition of the 
resulting minerals [e.g. 1,2]. On a planet-wide basis, the 
distribution of elements between the metallic core and 
the silicate portion of the planet is influenced by the 
f0 2 at the time of differentiation. Additionally, if the 
intrinsic f0 2 of a planet is known, the processes which 
may change that value can be investigated. Although 
magmatic f0 2 can reflect the degree of oxidation of the 
magma source region, it can also provide insight into 
processes, such as metasomatism, tfegassing, or as- 
similation, which may have operated on and in the 
magma as it moved to the surface. 

Recent work on the basaltic shergottites has pro- 
vided new measurements of the f0 2 recorded in the 
pyroxenes of these meteorites [3,4,5]. These pyroxenes 
exhibit ai f0 2 range of ~3-4 orders of magnitude [5]. 
Several processes have been proposed to explain the 
origin of these differences, the majority of which rely 
on assimilation [e.g. 6]. New data, however, is consis- 
tent with intrinsic fD 2 differences in the magma source 
region being responsible for the measured SNC f0 2 
variations. More detailed analyses of natural SNC py- 
roxenes are needed in order to better understand the 
source of the f0 2 range recorded in the basaltic sher- 
gottites. 

Determination of f0 2 : Because of the cumulate na- 
ture and inferred slow cooling of the SNC’s, there are 
difficulties in measuring f0 2 with traditional methods. 
Typically, f0 2 measurements are made using the 
Fe 2+ /Fe 3+ ratios of oxides. However, oxides must cool 
rapidly to lock in their crystallization f0 2 . It has been 
demonstrated that oxides can reequilibrate in the sub- 
solidus [7], leading to non-magmatic f0 2 estimates. For 
example, Fe-Ti oxides in the meteorite Shergotty record 
temperatures of 682°C, with a log f0 2 of -20 [8]. This 
temperature is not magmatic, implying the oxides ex- 
perienced subsolidus reequilibration thereby discount- 
ing the f0 2 calculation. A new technique utilizing elec- 
tron probe oxygen measurements and stoichiometric 
determinations has resulted in more accurate f0 2 meas- 
urements [9], although the compositions still appear to 
record post-magmatic conditions. Furthermore, oxides, 
except chromite, come on the liquidus late in most sher- 
gottite magma crystallization sequences, thereby not 
recording the early f0 2 history of the melt. Thus, for 
the currently available SNCs, Fe-Ti oxide f0 2 calcu- 
lations may not reflect magmatic conditions. 


Eu-oxybarometer. The f0 2 of a melt can also be 
calculated using the partitioning of Eu, relative to adja- 
cent REEs, into the first crystallized minerals [10,11]. 
The Eu distribution coefficient between melt and min- 
eral depends on the ratio of Eu 3+ /Eu 2+ in the melt, which 
is a function of the f0 2 conditions during crystallization 
[10]. Previous experimental studies have investigated 
the partitioning of Eu between augite and melt as a 
function of f0 2 in order to use this relationship to un- 
derstand natural systems [10,11]. We have applied this 
approach to a synthetic shergottite melt and low-Ca 
pyroxene assemblage in order to determine the f0 2 of 
the resulting magmas (Figure 1) [3,5]. Since pigeonite is 
known to be one of the first minerals to crystallize from 
a shergottite melt, the pigeonite Eu-oxybarometer pro- 
vides f0 2 estimates of the earliest magmatic history. 



Figure 1 . Distribution coefficient of Eu/Gd and Eu/Sm in pi- 
geonite in a dry shergottite melt vs. f0 2 . Error bars shown are 
calculated from the standard deviations between ion probe 
measurements of the same sample. 

The Eu-oxybarometer has an important advantage 
over other methods which rely solely on Fe measure- 
ments . Although the Fe 2+ /Fe 3+ ratio of a melt or mineral 
phase is a function of f0 2 , this ratio is known to be re- 
set through processes such as mineral reequilibration, 
metasomatism, or shock [7,12]. One or more of these 
processes have certainly affected the SNC meteorites. 
This resetting of ferrous/ferric ratios can lead to erro- 
neous f0 2 measurements as discussed above. The 
REEs, acting as trace elements, on the other hand, are 
essentially locked into a mineral during crystallization 
and are unlikely to be reset by later igneous or shock 
events. This leads to more robust crystallization t0 2 
measurements made using the Eu-oxybarometer tech- 
nique. Due to the small number and special nature of 
the SNCs (slowly cooled cumulate rocks subjected to 
high shock pressures), the measurement of f0 2 as re- 
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corded in the pyroxenes provides the most accurate 
measurements of magmatic oxidation state. 

Results: The new Eu-oxybarometer calibration for 
pigeonite yields refined f0 2 values for the basaltic 
shergottites, ranging from IW 3.5 for Shergotty to 
IW+0.2 for QUE94201. These numbers differ in value, 
but not trend, from those published previously [e.g. 
9,11]. These P0 2 values were calculated using natural 
pigeonite core data from [13]. 

In order to better constrain small variations in P0 2 
recorded in the pigeonites, new, more detailed analyses 
of natural phenocryst REE contents are needed. This 
data will be collected on the Woods Hole Oceano- 
graphic Institution 3F and 1270 ion probes. Detailed 
traverses of pigeonite cores will be carried out using 
the smallest beam diameter feasible. This will allow for 
a closer look into potential small-scale f0 2 variations in 
single meteorites and may shed light on the origin of 
the fD 2 differences exhibited by the SNCs . 

ft> 2 differences: The new results from the pigeonite 
Eu-oxybarometer detailed above show ~3-4 orders of 
magnitude difference in f0 2 between the basaltic sher- 
gottite meteorites. This variation is less than the ob- 
served terrestrial variation of ~9 orders of magnitude 
[14,15], but much greater than that seen on the moon. 
Several possibilities exist that could explain these dif- 
ferences. These could reflect intrinsic f0 2 differences 
in the magma source regions on Mars (i.e. the martian 
mantle or lithosphere) or they could result from interac- 
tions between the magma and an assimilant, that is as 
yet unsampled, either at depth or near the surface. This 
new f0 2 data may provide insight into which of these is 
the more likely explanation. 

Experimental studies have shown that pigeonite is 
one of the first minerals to crystallize from a shergottite 
melt, whereas Fe-Ti oxides come on the liquidus late in 
the melts crystallization sequence (after 30 and 60% 
crystallization in hydrous and dry experiments, respec- 
tively, to possibly 95% in Ios Angeles) [e.g. 16,17]. 
Therefore pigeonite core f0 2 values should record a 
different portion of magmatic history than the oxides. 
However, when the pigeonite f0 2 values are compared 
with the most recent determination of Fe-Ti oxide f0 2 
values a close correlation is observed (Figure 2). This 
implies that the melt fD 2 for each of these meteorites 
was buffered internally throughout much of the crystal- 
lization sequence. That is, the early crystallizing miner- 
als were exposed to the same relative f0 2 environment 
as the late crystallizing minerals. These results are 
consistent with the argument that the variations in f0 2 
recorded in the basaltic shergottites are the result of 
intrinsic f0 2 differences in the magma source regions of 
these rocks. However, assimilation cannot be com- 
pletely ruled out on this basis alone. If the assimilation 


occurred before any crystallization had taken place, the 
primaiy fC 2 recorded in the phenocrysts would be that 
of the source region modified by the assimilant. 



Figure 2. Comparison of f0 2 data from pigeonite Eu- 
oxybarometer and Fe-Ti oxide equilibria. Oxide data from [6]. 

Conclusions: We have presented a calibration of 
the pigeonite Eu-oxybarometer for the basaltic shergot- 
tites. This tool allows for the measurement of ft) 2 dur- 
ing the earliest crystallization history of the melt. The 
new f0 2 values in the pigeonites correlate well with 
recent Fe-Ti oxide data. This apparent f0 2 buffering 
during crystallization is consistent with the f0 2 varia- 
tions between the SNCs being the result of intrinsic fD 2 
differences in the magmas ’ source regions, the martian 
mantle or lithosphere. However, an assimilation event 
which took place before crystallization commenced 
cannot be ruled out. New, more detailed REE analyses 
of natural pigeonite cores are being collected now in 
order to explore the lull implications of these results. 
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Introduction: Mafic igneous rocks serve as probes 
of the interiors of their parent bodies - the composi- 
tions of the magmas contain an imprint of the source 
region composition and mineralogy, the melting and 
crystallization processes, and mixing and assimilation. 
Although complicated by their multifarious history, it 
is possible to constrain the petrologic evolution of an 
igneous province through compositional study of the 
rocks. Incompatible trace elements provide one means 
of doing this. I will use incompatible element ratios of 
martian meteorites [1] to constrain the early petrologic 
evolution of Mars. Incompatible elements are strongly 
partitioned into the melt phase during igneous proc- 
esses. The degree of incompatibility will differ de- 
pending on the mineral phases in equilibrium with the 
melt. Most martian meteorites contain some cumulus 
grains, but nevertheless, incompatible element ratios of 
bulk meteorites will be close to those of their parent 
magmas. ALH 84001 is an exception, and it will not 
be discussed. The martian meteorites will be consid- 
ered in two groups; a 1.3 Ga group composed of the 
clinopyroxenites and dunite, and a younger group 
composed of all others. 

Planetary Comparisons: On Earth, P and La are 
highly correlated in igneous rock suites, resulting in a 
small range in P/La ratios for terrestrial mafic magmas 
(Fig. 1). Magmas with higher incompatible element 
contents have lower P/La, indicating La is more in- 
compatible than P. This is also supported by various 
estimates of the composition of the bulk continental 
crust [2], which is more highly enriched in La than P 
compared to an estimated bulk silicate Earth [2]. Ma- 
fic rocks from the Moon follow the terrestrial example 
- P/La ratios exhibit a narrow range (Fig. 1). In con- 
trast, martian meteorites show a wide range in P/La 
(Fig. 1). The planetary distinctions are also clearly 
shown on a P vs. Yb diagram (Fig. 1). The terrestrial 
mafic rocks show no correlation between P and Yb, 
lunar rocks occupy a band of increasing P and Yb, 
while in martian meteorites these elements are strongly 
correlated. This suggests a fundamental difference in 
the petrologic evolution of Mars as compared to the 
Earth or Moon. The constancy of Yb/P implies that the 
partitioning of these elements is governed by a single 
phase, or two phases either in constant proportions or 
with similar P/Yb partition coefficient ratios. The 
phase or phases would have a much lower partition 
coefficient for La, such that P and La are decoupled. 


• M an • Earth A Moon 



• Mars •Earth A Moon 
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Figure 1 . P-La and P-Yb plots for mafic igneous rocks 
from Mars, the Earth and Moon. 

Crustal Assimilation: Recently, a number of geo- 
chemical characteristics of a subset of martian meteor- 
ites have been ascribed to assimilation of crustal mate- 
rial by mafic magmas [e.g. 3]. The P/La ratios for mar- 
tian meteorites <1 Ga in age are positively correlated 
with their e 143 Nd [1] (Fig. 2), a presumed measure of 
crustal assimilation/contamination, indicating that this 
process may have contributed to the scatter in P/La. 
This is in general in accord with terrestrial geochemis- 
try - the P/La ratio for the bulk continental crust is 
lower than that estimated for the bulk silicate Earth [2]. 
However, P and La are enriched in the continental crust 
by roughly lOx and 30x, respectively [2]. QUE 94201 
has the highest e 143 Nd (Fig. 2) and the highest P con- 
tent (Fig. 1). This rock is considered to have suffered 
minimal crustal contamination [see 3]. Hence, if ba- 
salts like Shergotty were formed by contamination of 
primary or evolved magmas, then these magmas, and 
possibly their source regions, must have had incom- 
patible element contents different in detail from those 
of QUE 94201. Thus, for example, their Sm/Nd ratios 
may have been lower, and the source region would then 
have had lower e 143 Nd. Regardless of complications 
from possible assimilation of crust, the P/La ratios of 
martian mafic igneous rocks show wide ranges, while 
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the P/Yb ratios are roughly constant, quite unlike the 
cases for the Earth or Moon. 


• <1 G* 01.3G* 



Figure 2. Correlation of P/La with e l43 Nd in martian 
meteorites. Q = QUE 94201; b = martian basalts Los 
Angeles, Shergotty and Zagami. 

Incompatible Element Ratios: Martian meteorites 
1.3 Ga in age do not follow the trend exhibited by the 
younger rocks (Fig. 2). These rocks have by far the 
lowest P/La ratios, yet intermediate e 143 Nd. They are 
also distinct in other trace element ratios (Fig. 3). 
These distinctions are not due to crustal assimilation - 
the terrestrial continental crust has Cl-normalized 
Hf/Sm of ~1.3 and Ta/La of -1.0 [2] and it is likely 
that the martian crust is also unfractionated in these 
ratios from the primitive Mars composition. Thus, it is 
unlikely that the low ratios observed for the 1.3 Ga 
martian meteorites could be a result of contamination 
by martian crustal rocks. The 1 .3 Ga martian meteor- 
ites and QUE 94201 have the lowest and highest P/La 
ratios, respectively. 
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Figure 3. Ta/La vs. Hf/Sm for martian meteorites 
showing the distinction between the 1.3 Ga and 
younger rocks. Q = QUE 94201 . 

Modeling Mars: Using the experimental high 

pressure mineralogy [4] of a Dreibus-Wanke Mars 
composition and a variety of partition coefficients [1], I 
have attempted to model the martian meteorite trace 
element signatures as arising from a sequentially de- 
pleted primitive mantle (i.e. fractional fusion). Models 
were done for pressures in the garnet lherzolite and the 
majorite-clinopyroxene-P-phase stability regimes. In 


no case did either QUE 94201 or the 1.3 Ga martian 
meteorites fall within the envelop of possible melt 
compositions in any of the tested models. In short, 
sequential melting of a Dreibus-Wanke martian mantle 
composition under any P-T regime will not yield melts 
with incompatible element ratios that come close to 
those of either the 1.3 Ga martian meteorites, nor of the 
most pristine younger volcanic rock, QUE 94201 . Al- 
ternative models must be sought. 

High Pressure Cumulates: I have begun evaluat- 
ing high pressure magma ocean cumulates as alterna- 
tive source regions for the martian meteorites. The 
garnet-like REE pattern of QUE 94201 suggests that 
cumulus garnet may be required. However, garnet 
sensu stricto is not a near-liquidus phase in ultramafic 
melt compositions, and thus is not likely to be a cumu- 
lus phase. At P >15 GPa, majorite does become a liq- 
uidus or near liquidus phase in ultramafic melts [5]. 
Limited partition coefficient data [1] indicate that ma- 
jorite crystallizing from a martian magma ocean should 
have a garnet-like REE pattern, and thus could impose 
a gamet-like signature to the cumulates. In addition, 
very limited P data suggest that majorite/melt K d s for P 
and Yb are similar at -1. Thus, cumulus majorite 
could define the P/Yb of the martian meteorite source 
regions, and P and La would be decoupled. Finally, 
measured Hf/Sm partition coefficient ratios for ma- 
jorite are >2, although Hf is incompatible in majorite 
(kas -0.5-0.9). Hence, it appears plausible that the 
parent magma of QUE 94201 originated by remelting a 
magma-ocean-cumulate source region in which ma- 
jorite was a cumulus phase. However, a serious stum- 
bling block might be Al. Martian basalts are noted for 
being Al-poor [see 3], while a majorite-rich source 
region should be Al-rich. 

The 1.3 Ga martian meteorites must have come 
from a different source in order to explain the different 
P/La, Hf/Sm and Ta/La ratios (Figs 1, 3). Other possi- 
ble high P phases are Mg-perovskite, y-spinel, P-phase 
and clinopyroxene. I will be evaluating these as possi- 
ble cumulus phases in addition to majorite to explain 
the geochemical characteristics of the 1 .3 Ga martian 
meteorites. 

References: [1] References for martian meteorite 
data, and high P partition coefficients can be found at 
http://ares.jsc.nasa.gov/mittlefehldt/. 
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102, 5251. [5] Ohtani et al (1989) Contrib. Mineral 
Petrol, 103, 263. 
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Introduction: Chemical and (oxygen) isotopic com- 
positions of SNC meteorites have been used by a num- 
ber of workers to infer the nature of precursor materials 
for the accretion of Mars [1-4]. The idea that chondritic 
materials played a key role in the formation of Mars has 
been the central assumption in these works. Wanke 
and Dreibus [1] have proposed a mixture of two types 
of chondritic materials, differing in oxygen fugacity but 
having Cl type bulk chemical composition for the non- 
volatile elements, for Mars’ precursor. But a number of 
studies based on high pressure and temperature melt- 
ing experiments do not fivor a Cl type bulk planet 
composition for Mars, as it predicts a bulk planet Fe/Si 
ratio much higher than that reported from the recent 
Pathfinder data (e.g., [5]). Oxygen forms the bulk of 
Mars (~40 % by wt., e.g., [1]), and might provide clues 
to the type of materials that formed Mars. But models 
based on the oxygen isotopic compositions of SNC 
meteorites predict three different mixtures of precursor 
materials for Mars: 90% H + 10% CM [2], 85% H + 1 1% 
CV + 4% Cl [3] and 45% EH + 55% H [4]. As each of 
these models has been shown to be consistent with the 
bulk geophysical properties (such as mean density, and 
moment of inertia factor) of Mars, the nature of the ma- 
terial that accreted to form Mars remains ambiguous. 

The N and O isotopic model: We have adopted a 
new approach to this study, which involves the iso- 
topic compositions of nitrogen (a 15 N) and oxygen 
(A I7 0) in SNC meteorites. As we have shown earlier [6], 
the nitrogen and oxygen isotopic systematics of Mars 
mantle derived from the SNC data (5 15 N = -30%o [6]; 
A 1 7 0 ~ 0.321 ±0.01 3%o [7]) are consistent with a mixture 
of EC and OC type precursor materials for the planet. A 
similar result is obtained if one considers A 17 0 and ra- 
diogenic chromium (e 53 Cr) data for the SNCs [8]. More- 
over, such a mixture is also consistent with the bulk 
planet Fe/Si ratio for Mars as predicted from the mo- 
ment of inertia factor obtained from the Mars Pathfinder 
data [5]. We would like to point out here that the rela- 
tive contributions of EC and OC in our earlier calcula- 
tions [6] were based the A 17 0 of EH and H. But as we 
have already pointed out, one cannot actually resolve 
such a mixing proportion from those involving the 
other types of EC and OC by the N and O isotopic sys- 
tematics. Therefore we repeated the calculations by 
using the mean values for the EC [9] and OC [10], which 
suggested the EC:OC ratio to be 74:26 (Fig. 1). We also 
tried to further resolve the possible mixing proportions 
between EC and OC by using the bulk planet Fe/Si for 


Mars predicted by Bertka and Fei [5] as a constraint. 
But such an exercise did not prove encouraging. 



Fig. L The nitrogen and oxygen isotopic systematics 
of Mars, in relation to those of different types (inset) of 
chondrites (after [6]). Mars is consistent with a mixture 
of EC and OC in a ratio of 74:26 (main plot). 

Therefore we derived the bulk chemical composition 
for Mars by mixing 74% EC and 26% OC, and investi- 
gated the consequences of such a model. The bulk 
planet Fe of 24.7 wt % and Fe/Si ratio (by wt.) of 1.39 
obtained from the above mixture are consistent with the 
ranges given by Bertka and Fei [5]. It is possible to 
reproduce the bulk chemical composition (major oxides) 
of the martian mantle (Table 1) and its mean density as 
suggested by the SNCs [1] if one accommodates a 
small fraction (~0.08) of the planet’s Si in the core. The 
core predicted by the present model thus has 6.7 wt % 
of Si, a density of 6.5 g/cc and accounts for 23 % of the 
planet’s mass. The moment of inertia factor of 0.368 ± 
0.003 calculated for Mars, by assuming a crust with 
thickness of 20 to 50 km and density of 2.7 to 3 g/cc, is 
consistent with the value predicted from the Mars Path- 
finder data [11]. 

Si in martian core: Si has commonly not been in- 
voked as a light element in the core of Mars [12] al- 
though its presence in the terrestrial core has often 
been suggested from geochemical arguments [13]. 
Wanke and Dreibus [14] have pointed out that the de- 
pletion of Si in Earth’s mantle -derived rocks as com- 
pared to samples from Mars can be explained by the 
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presence of 14.4 wt % of Si in Earth’s core. Although, 
the presence of significant Si in Earth below core- 
mantle boundary is a matter of debate at the present, 
recent results from high pressure and temperature melt- 
ing studies [15] and theoretical studies based on ab 
initio calculations and seismic data [16] suggest as low 
as 8 mol % of Si in Earth’s outer core. Iron sillicide as a 
possible host phase for Si in the core has thus been the 
focus of a number of experimental and theoretical stud- 
ies. The results from Guyot et al. [17], based on both 
experimental and theoretical observations, show that 
significant amounts of Si can be dissolved in Fe be- 
tween the pressures 8 and 15 GPa. A number of similar 
studies have also shown that at pressures below 25 
GPa liquid iron can accommodate significant amounts 
of Si in solution, while at pressures above 25 GPa it (Si) 
may exsolve from the solution [13]. The expected pres- 
sure at the core -mantle boundary of Mars for the pre- 
sent model is about 20 GPa, which falls within the re- 
gime of significant Si solubility indicated by the above 
studies. Therefore, we cannot neglect the possibility of 
Si being present as a light element in the martian core. 
The partitioning of Si into the martian core would re- 
quire a reducing environment to prevail during the 
core -mantle separation, which is supported by the re- 
ducing nature of the martian mantle as has been re- 
cently pointed out [18-19]. The possible partitioning of 
Si into the core is also compatible with the chemical 
data on the opaque phases (metal bearing) in enstatite 
chondrites. A recent study [20], for example, has shown 
up to 12.5 wt % of Si in such phases. 

References: [1] Wanke H. and Dreibus G. (1988) 
Phil Trans . R. Soc. Lond., A3 25, 545-557. [2] Delaney J. 
S. (1994) Meteoritics , 29, A459. [3] Lodders K. and Fe- 
gley B. Jr. (1997) Icarus, 126, 373-394. [4] Sanloup C. et 
al. (1999) Phys. Of Earth Planet. Interior., 112, 43-54. 
[5] Bertka C. M. and Fei Y. (1998) Science, 281, 1838- 
1840. [6] Mohapatra R. K. and Murty S. V. S. (2001) 
Meteoritics & Planet. Sci., 36, A138. [7] Franchi I. A. 
et al., (1999) Meteoritics Planet. Sci. 34, 657-661. [8] 
Lugmair G. W. and Shukolyukov A. (1998) Geochim. 
Cosmochim. Acta, 62, 2863-2886. [9] Newton J. et al. 
(2000) Meteoritics & Planet. Sci., 35, 689-698. [10] 
Clayton R. N. et al. (1991) Geochim. Cosmochim. Acta, 
55, 2317-2337. [11] Folkner W. M. (1997) Science, 278, 
1749-1752. [12] Stevenson D. J. (2001) Nature, 412, 214- 
219. [13] Hilgren V. J. et al. (2000) In Origin Earth and 
Moon (eds. R. M. Canup and K. Righter), Univ. Arizona 
Press, Tuscon, pp. 245-263. [14] Wanke H. and 
Dreibus G. (1997) IPS XXVIII, 1495-1496. [15] Dobson 
D. P. et al. (2002) American Miner., 87, 784-787. [16] 
Alfe D. et al. (2002) EPSL, 195, 91-98. [17] Guyot F. et 
al. (1997) European Jour. Miner., 9, 277-285. [18] Wad- 


hwa M. (2001) Science, 291, 1527-1530. [19] Herd C. D. 
K. et al. (2002) Geochim. Cosmochim. Acta, 66, 2025- 
2036. [20] Lin Y. and El Goresy A. (2002) Meteoritics. 
Planet. Sci., 37, 577-599. 

Table. 1. Bulk features of the mantle and core of Mars 
obtained from the present model, by partitioning ~8% 
of the planet’s Si into the core. Shown for comparison 
is the Wanke and Dreibus (1988) model (WD). 



This work 

WD [1] 

Mantle+crust 



Major Oxides (wt %) 



MgO 

29.1 

30.2 

ai 2 o 3 

2.5 

3.02 

Si02 

47.1 

44.4 

CaO 

1.9 

2.45 

Na 2 0 

1.2 

0.50 

Ti02 

0.1 

0.14 

Cr 2 0 3 

0.6 

0.76 

MnO 

0.4 

0.46 

FeO 

16.9 

17.9 

Density (g/cc) 

3.55* 

3.52 

Core 



Fe 

68.0 

77.8 

S 

17.8 

14.2 

Ni 

6.6 

7.6 

Co 

0.30 

0.36 

Si 

6.7 

0 

P 

0.6 

0 

Mass (%) 

23.0 

21.7 

Density (g/cc) 

6.5 

6.8 


* Refers to bulk density of mantle only. 
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Introduction: Redox variations have been reported 
among the shergottites [1,2,3]. [1] used the Eu and Gd 
partitioning experiments of [4], designed for the 
LEW86010 angrite, to infer a range of fo 2 for the sher- 
gottites. [2,3] inferred fo 2 using equilibria between Fe-Ti 
oxides. There is fairly good agreement between the Fe- 
Ti oxide determinations [2,3] and the estimates from Eu 
anomalies [1] in terms of which meteorites are more or 
less oxidized. The Eu anomaly technique and the Fe-Ti 
oxide technique both essentially show the same trend, 
with Shergotty and Zagami being the most oxidized and 
QUE94201 and DaG 476 being the most reduced. Thus, 
the variation in fo 2 appears to be both real and substan- 
tive. However, although the redox trends indicated by 
the two techniques are similar, there is as much as a 
two log unit offset between the results of [2,3] and [1]. 
One explanation for this offset is that the Eu calibration 
used for the shergottites was actually designed for the 
LEW86010 angrite, a silica-undersaturated basalt 
whose pyroxene (diopside) compositions are rather 
extreme. To correct this, [5] conducted a set of experi- 
ments on the redox relationship of Eu partitioning rela- 
tive to Sm and Gd for pyroxene and melt compositions 
more relevant to Martian meteorites. We have taken 
the experimentally determined relationship between 
D Aug ite /md. ( Sm Eu Gd ) ^ fo 2 of [5] and applied it to 
augite Eu and Gd data for the basaltic Shergottites of 
[1] as well as previously unreported augite Sn data 
(analytical methods identical to [1]). 

Figure 1: 


Data Set The results of the controlled f0 2 partition- 
ing experiments of [5] are shown in figure 1 . Table 1 


shows the Eu and Gd augite core data of [1 ] along with 
previously unreported Sm augite core data. Iking 
using these data and whole -rock Sm, Eu, and Gd data 
[6], we have calculated the D^Dod and IhJD Sm ratios 
for selected basaltic Shergottites (also Table 1). Ap- 
plying these D ratios to the curves in Figure 1, we de- 
termined the f0 2 for the meteorites. These are given in 
Table 2 along with f0 2 ’s determined from Fe-Ti oxides 
[2,3] for comparison. 



Discussion and Conclusions: Figure 2 shows a 
plot of f0 2 for Shergottites determined from Efc u /D Gd 
both from this study and [1] versus f0 2 as determined 
from Fe-Ti oxides [2,3]. The f0 2 ’s determined from 
DE U /D Gd using the calibration of [5] are offset to higher 
values at higher f0 2 ’s compared to those determined 
from E^u/Dcd [1] using the calibration of [4] but to 
slightly lower values at the lowest f0 2 ’s. It appears 
that for some of the Shergottites (Shergotty, Zagami 
and EET 79,1 A), some of the discrepancy between 
fC 2 ’s determined using the Fe-Ti method and that using 
the D Eu /D Gd method was due to the need to calibrate the 
relationship between D Eu /D Gd and f0 2 for Martian basalt 
and pyroxene compositions. However, for all of the 
Martian meteorites considered here except Shergotty, a 
discrepancy remains between the f0 2 determined from 
Deu/D od and that determined from Fe-Ti oxides. That 
discrepancy may in fact be real in the sense that the 
two methods are recording ft) 2 at different stages in the 
crystallization of the magma. However, before we ac- 
cept this conclusion, other factors must also be con- 



f02 (log ( AIW)) 
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sidered. In particular, two important considerations 
may be exemplified by Zagami and EET79,1 A. 

(1) How representative are the whole rocks?: In 
the case of Zagami, whole rock data used to estimate 
Deu/Dree ratios [6] are highly variable. This variability 
is evident in separate analyses conducted in the same 
laboratories as well as between those performed in 
separate labs. These analyses are mostly INAA meas- 
urements conducted on small chips. The variability of 
the modal compositions of the chips (in particular, the 
relative phosphate and pyroxene abundances) are likely 
the cause of the variability of the in the REE ratios. 
Therefore, care needs to be taken in selecting the ap- 
propriate analyses for estimating the whole -rock REE 
ratios. 

(2) Are the augite cores in equilibrium with the melt 
as represented by the whole -rock REE ratios?: This 
could be particularly problematic for EET79, 1 A. Augite 
may have formed quite late in the crystallization of its 
parent magma. In fact, equilibrium crystallization ex- 
periments conducted on the EET79,1A composition at 
successive temperatures between its liquidus and 
solidus, produce no pyroxenes of augite composition 

[7]. 

Figure 3: 



f02 from DEu/DSm Aug Res (log (DIW)) 


Figure 3 shows a plot of f0 2 for shergottites de- 
termined using Efc u /D Gd versus that determined using 
CWDsm both from this study. Agreement is not found 
for the f0 2 ’s determined with these two ratios for Za- 
gami and EET79,1A. The f0 2 ’s determined by the 


£WDree method may be suspect for these meteorites 
for the reasons stated above. Agreement is seen within 
error for the f0 2 ’s determined with these two ratios for 
Shergotty, Los Angeles, EET79,1B and QUE94201 . We 
might, therefore, conclude that the fCVs determined by 
the EWDree method for these rocks is the correct one 
for their parent magmas at the early stages of magma 
crystallization. However, before accepting this conclu- 
sion, further corroboration is necessary. Specifically, 
employment of Deu/Dree for pigeonite and melt is sug- 
gested. [8] have begun to calibrate the relationship 
between E^/Dree and f0 2 for pigeonite/melt at Martian 
compositions. Up until now, however, REE abun- 
dances on pigeonite cores have not been of suitable 
precision to employ this method with any certainty. 

References [l]Wadhwa, M. (2001) Science 
291, 1527-1530. [2] Herd C. and Papike J.J. (2000) MAPS 
35, A70. [3] Herd C. et al. (2001) Geochim. Cosmochim. 
Acta 66, 2025-2036. [4]McKay G.A et al. (1994) Geo- 
chim. Cosmochim. Acta 58,2911-2919. [5] Musselwhite, 
D.S. and Jones, J.H. (2002) MAPS 37, A101 [6] Meyer, 

Table 1: SIMS analyses of Augite cores and D 

Sm ppb Eu ppb Gd ppb DEu/DSm DEu/DGd 


SHERGOTTY164 ±14 

58 ±3 

315 ±16 

0.88 ±0.11 

0.75 ±0.0 

ZAGAMI 

369 ±31 

103 ±8 

670 ±59 

0.58 ±0.15 

0.60 ±0.01 

LOS ANGE 

456 ±35 

140 ±11 

946 ±79 

0.75 ±0.11 

0.57 ±0.0 

EET 79,10 

214 ±22 

56 ±4 

374 ±28 

0.53 ±0.07 

0.59 ±0.0 

EET 79,1(E 

240 ±27 

73 ±5 

630 ±47 

0.65 ±0.09 

0.47 ±0.0 

DaG 476 


62 ±5 

607 ±46 


0.50 ±0.0 

QUE 94201 

239 ±18 

56 ±4 

552 ±28 

0.52 ±0.06 

0.44 ±0.01 


C. (1998? NASA Johnson Space Center Publication 
27672, Revision A, 88. [7] Jones, J.H. personal commu- 


Table 2: fp for selected Shergottites determined by various methods 

from Fe-Ti 

from Eu/Gd 

from Eu/Gd 

from Eu/Sm 

[2,3] 


[i] 


This Study 


Shergotty 2.6 ±o.i 

1.9 

+0.9/-0.7 

2.5 

+0.9/-0.7 

1.8 

+ 1.2/-1.0 

Zagami 2.4 ±o.3 

0.4 

+0.V-0.7 

0.8 

+0.8/-0.7 

-1 

+0.9/-1.0 

Los Angel 2.4 ±0.3 

0.2 

+0.6/-0.7 

0.4 

+0.6/-0.7 

0.6 

±1.0 

EET 79,1/ 1.7 ±0.2 

0.3 

+0.9/-0.8 

0.7 

+Q.9/-0.8 

- 1.9 

±1.1 

EET 79, IB 1.7 ±0.2 

- 0.7 

+0.6/-0.8 

- 0.7 

+0.6/-0.8 

- 0.3 

+0,8/- 1.1 

DaG 476 1.2 ±0.4 

- 0.5 

+0.7/-0.9 

- 0.5 

+0.7/-0.9 



QUE942C 1 ±0.8 

o 

> 

+0.7/-0.9 

- 1.5 

+0.7/-0.9 

-2 

+1.0/- 1.2 


nication. [8] McCanta, M.C. et al. (2002) MAPS 37, A97. 
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Introduction: The shergottite Martian meteorites 
present a variety of oft -confusing petrologic features. 
In particular, represented among this subgroup are ba- 
salts with very depleted LREE abundances, as well as 
those with nearly chondritic overall REE abundances. 
The LREE-depleted basalts appear to more closely re- 
cord the REE and isotopic features of their mantle 
source regions. Those basalts with more nearly chon- 
dritic REE abunandances appear to contain an extra 
component often referred to as a “crustal” component. 
The addition of the crustal component tends to restore 
the overall REE abundance pattern towards chondritic 
relative abundances. Here we suggest that the crustal 
component could derive from “andesitic” rocks ob- 
served remotely to occur on the Martian surface, and 
which were analysed at the Pathfinder site. 


Martian Meteorites 
with Pathfinder Soils and Rocks 



Figure 1. Mg/(Mg+Fe) vs. K/(Ca+Na+K) for Mar- 
tian meteorites and Pathfinder soils and rocks. 

Comparison of Martian meteorites to Pathfinder 
rocks and soils: Figure 1 compares Mg/(Mg+Fe), i.e., 
mg-values, and K/(Ca+Na+K) for Martian meteorites to 
those in Pathfinder rocks and soils [1]. The mg-values 
are related to olivine abundances in the shergottites. 
“Picritic basalts” and “peridotitic shergottites” are en- 
riched in olivine compared to the “basaltic shergot- 
tites”. However, it is clear that the shergottites cannot 
be adequately characterized by a single parameter. 
Also, K/(Ca+Na+K) cannot be identified with a single 
mineralogical feature of the samples, but was chosen 
partly to reflect the feldspar composition of the sam- 
ples, partly to reflect the feldspar composition in possi- 
bly assimilated precursor materials, and partly because 


K has geochemical affinities to both La and Rb. Thus, 
this parameter might be expected to correlate to the Sm- 
Nd and Rb -Sr isotopic systematics of the meteorites. 

The meteorite and Pathfinder data fields are com- 
pletely separated in Figure 1 . However, the peridotitic 
shergottites and picritic basalts lie on the extension of 
the trend of Pathfinder data, and thus could be inter- 
preted as capable of providing a mafic component to 
the Pathfinder soils. The need for such a comp onent, 
not represented among the analysed Pathfinder rocks, 
was recognized by Wanke et al. [1]. 

Preliminary Nd and Sr isotopic data for NWA 
1068: In Figure 1, the picritic basalt displaced most 
strongly towards the pathfinder data is NWA 1068. It 
contains large olivine megacrysts characteristic of 
strongly LREE-depleted picritic basalts like DaG 476 
and 489, SaU 005, and EET79001, but it has nearly 
chondritic relative REE abundances like basaltic sher- 
gottites, Shergotty, Zagami, and Los Angeles [2]. We 
have begun Sr and Nd isotopic investigations of NWA 
1068. Comparing preliminary Sm-Nd data to that of Los 
Angeles (LA), another hot-desert meteorite, shows the 
effect of terrestrial contamination on analyses of an 
unleached whole rock (WR) and a leachate from a sec- 
ond whole rock (WR(1)). Sm-Nd of the leach residue, 
WR(r), are very similar to leach-residue data for LA, 
and lie on the ~172 Ma isochron determined for LA [3]. 
The WR and WR(1) data are displaced from WR(r) 
along a mixing line terminating at the isotopic data for 
caliche physcially removed from LA. The slope of this 
mixing line corresponds to a maximum apparent age of 
245 Ma. The initial Rb-Sr isotopic data show analo- 
gous relationships, with the WR(r) datum falling in the 
range of values for ~170 Ma isochrons for Shergotty, 
Zagami, and Los Angeles. Again, the WR and WR(1) 
data are displaced along a mixing line terminating at the 
caliche datum. For this discussion, we assume an age 
of 173 Ma for NWA 1068, characteristic of shergottites 
with similar REE abundances [4], and calculate initial 
I43 Nd/ 144 Nd and 87 Sr/ 86 Sr ratios from the WR(r) data. 
Additonal isotopic data will be presented at the work- 
shop. 

Olivine accumulation and fractional crystallization 
in shergottites: Figure 2 shows that magma -chamber 
and/or near-surface processes of olivine accumulation 
and fractional crystallization have been operative on 
shergottite magmas of differing isotopic compositions. 
These effects have operated independently of whether 
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Figure 2. Illustrating the effect of olivine 
accumulation and fractional crystallization, 
resp., on shergottite magmas. 

the isotopic systematics were inherited from the Mar- 
tian mantle directly (meteorites with large, positive val- 
ues of £N d ), or likely represent mixing of mantle and 
crustal components (meteorites with negative ^ d ). Oli- 
vine accumulation is manifest in the presence of large 
megacry sts of olivine in the picritic basalts NWA 1068, 
EET79001 A, Dhofar 019, SaU 005, and DaG 476, 489. 
The effect of fractional crystallization is manifested in 
low Mg/(Mg+Fe) for several of the basaltic shergot- 
tites, including especially Los Angeles and QUE94201 
(Figure 1). Los Aigeles has “crustal” isotopic sys- 
tematics, whereas QUE94201 has “mantle” isotopic 
systematics (Figure 3). Olivine fractionation clearly 
occurred after the isotopic and REE characteristics of 
the magmas were established. 

Isotopic mixing relationships: We presented two 
component isotopic mixing models for the initial Nd- 
and Sr-isotopic composition in EET79001 and Los An- 
geles [5]. The simplest scenario, binary mixing, is illus- 
trated in Figure 3, updated from [5]. 

Because late-stage magma -chamber and near- 
surface fractionation processes affect the abundances 
of Sr and Nd in the basalts, but not their isotopic com- 
positions, we allow the abundances of Sr and Nd in the 
mixing endmembers to vary within reasonable limits. 
The assumed isotopic characteristics of the crust are 
fixed by requiring the crustal component in LA to be 
<20%. To match the isotopic characteristics of LA, a 
highly fractionated basalt, we use the relatively high Sr 
and Nd abundances in QUE94201 for the “mantle- 
derived” basalt mixed with a crustal component with 
-30X higher Nd, and 4.4X higher Sr abundances. The 
parameter K/(Ca+Na+K) = 0.0033 for QUE, -30X less 
than for the Pathfinder soils and rocks, so the required 
level of Nd enrichment (-1 17 X Cl) may be readily avail- 
able in Martian crustal rocks. 

Figure 4 shows the results of an AFC model using 
the same endmembers as for Figure 3. A high ratio of 


Basaltic and Picritic Shergottites 
(Mantle-"Crust M Binary Mixing Model) 



Figure 3. Isotopic binary mixing models for man- 
tle and crustal components of fixed isotopic com- 
position present in differing abundances. 

the rate of assimilation to the rate of crystallization, r - 
0.3, is shown for illustration. The amount of assimilant 
is -3X higher than for the binary mixing model. The 
resultant is enriched in Nd and Sr by -3X for LA, but 
by only -20% for EET79001B. Thus, the models are 
equivalent for EET79001B, not for LA. 


Basaltic and Picritic Shergottites 

(Mantle-"Crust" AFC Model) 



Figure 4. AFC models for the same mantle and 
crustal components as for Figure 3. 

References: [1] Wanke H. et al. (2001), Chron. & 
Ev. Mars , ISSI, 317-330. [2] Barrat J. A. et al. (2002) Geo- 
chim. Cosmochim. Acta (submitted). [3] Nyquist L. E 
and Reese Y. D. (2001) Eos Trans . AGU, 82 (47), Fall 
Mtg. abs. P51A-02. [4] Nyquist L. E. et al. (2001) Chron . 
& Ev. Mars , ISSI, 105-164. [5] Nyquist L. E. et al. (2001) 
LPS XXXII, abs. #1407. 
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Introduction: As part of an ongoing study 
of the shergottites [1-3], we have measured the 
isotopic signature and concentration of xenon in 
mineral separates of Shergotty and EET79001 
Lithology-B. Similarities between the meteor- 
ites include the enhanced concentration of mar- 
tian atmospheric xenon in opaque and maske- 
lynite minerals relative to pyroxene [2]. The 
interior component appears to be present in all 
minerals but is best defined in pyroxene sepa- 
rates. This component consists of isotopically 
solar xenon with a significant contribution from 
fission xenon and is similar to those observed in 
Chassigny, ALH84001 and Nakhla [4]. The 
proportion of fission xenon in the pyroxene 
separate of EET79001 is lower than that in 
Shergotty. Pyroxene separates from Shergotty 
also have a well defined 129 Xe/ 132 Xe ratio higher 
than solar. 

One explanation for these variations and 
which we were seeking the effects <f in this 
study, is that they reflect the degree of incorpo- 
ration of ‘crustal’ material into the parent melt 
of the shergottites. The covariation of oxygen 
fugacity and e N d among the Shergottites has 
been argued to reflect different degrees of as- 
similation of crustal naterial (or other material 
with similar geochemical properties) in the par- 
ent melt. Geochemical, both plutonium (the 
most strongly implicated precursor of fission 
xenon in the martian neteorites) and uranium 
would be expected to be associated with the 
high eNd endmember, and so 136 Xe/ I32 Xe might 
be expected to be higher in Shergotty than in 
EET79001, as observed. 

Here we present preliminary results from a 
simple model tracking xenon isotopically 


through differentiation, outgassing and atmos- 
pheric loss of Mars. We require an atmosphere 
with elevated 129 Xe xs / 130 Xe and low 
I36 Xe*/ 129 Xe xs ; in addition, we attempt to pro- 
duce two interior reservoirs, one consisting of 
solar xenon with no radiogenic xenon and one 
with a high ratio of fissiogenic I36 Xe*/ 130 Xe and 
a low 129 Xe*s/ 136 Xe* ratio to account for the 
variable interior component. 

Model of Martian Xenon Reservoir 
Evolution; The model uses a simple rate equa- 
tion approach to trace the passage of *non 
and its parent radiogenic isotopes 129 I (t ]/2 = 17 
Myr), 244 Pu (ti/ 2 = 82 Myr), and 238 U (t 1/2 = 
4.468 Gyr) between reservoirs and radioactive 
decay. Rates of differentiation, degassing and 
atmospheric loss can be varied as free parame- 
ters to investigate under what circumstances 
reservoirs are qualitatively similar to those nec- 
essary to explain martian meteorite xenon com- 
ponents. Initial concentrations and ratios are 
similar to those of Swindle and Jones [5]. Xe, I 
and Pu are assumed to have behaved as n- 
compatible elements during differentiation, with 
the exception that xenon could also be de- 
gassed from the crust - xenon appears to have 
behaved as an incompatible element during the 
formation of the nakhlites [6, 7]. 

Results and Discussion As expected, the 
requirement to produce a crust with high 
136 Xe/ 130 Xe and low 129 Xe/ 130 Xe, and with sig- 
nificant concentrations of 136 Xe* imposes quite 
stringent constraints on its degassing history. 
Early major loss of 130 Xe is required to allow 
the 136 Xe/ 130 Xe ratio to evolve upwards subse- 
quently. This loss must have continued on time- 
scales long compared to that characteristic of 
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I29 I decay (>~100Ma) to allow high 
136 Xe*/ 129 Xexs ratios. In order to preserve sig- 
nificant concentrations of 136 Xe*, this degassing 
rate must have declined drastically on time- 
scales short compared to that characteristic of 
244 Pu decay (<~800 Ma). The broad charac- 
teristics of the martian atmospheric xenon iso- 
topic composition can also be reproduced 
within this model, using the atmospheric loss 
history proposed by Pepin [8, 9]. 

The identification of a pure solar component 
in the Chassigny meteorite argues that, if any- 
thing, our model underestimates the initial xenon 
to iodine/plutonium ratio. Since degassing of the 
source region rich in plutonium is thus crucial to 
the model’s success (decay of Pu can only ki- 
crease 136 Xe/ 130 Xe once the original comple- 
ment of solar xenon has been effectively b- 
moved), should our interpretation withstand 
further tests constraints are imposed on the lo- 
cation of this reservoir. Although labeled 
‘crustal’, its identification with the martian crust 
remains controversial - a mantle reservoir rich 
in incompatible elements is also a candidate 
source region. However, the requirement for 
the source region to be degassed argues in fa- 
vor of a location close to the surface, in effect 
strengthening the case for identification with the 
crust itself. 

Extending the model to the other major ra- 
diogenic rare gas component has presented 
some problems. In our model, the same prop- 
erties that lead to elevated 136 Xe* in the crust 
lead to elevated 40 Ar* excesses. As yet, no 
interior component enriched in 40 Ar from 40 K 
decay has been identified. Notably, K-Ar (to- 
tal 40 Ar/K) ages of the nakhlites are virtually 
identical to Ar-Ar ages and crystallization ages 
in other systems (1.3 Ga) [10] arguing for com- 
plete degassing of argon from the parent melt or 
its source shortly before formation of the mete- 
orites. This is consistent with the proposal that 
they formed in a near- surface flow, but is diffi- 


cult to reconcile with the presence of inherited 
136 Xe*. We note in passing that this must cast 
some doubt on models relying on the total re- 
tention of xenon during formation of the 
nakhlites [7]. The resolution of this conundrum 
remains unclear. 

Summary: Although work is at an early 

stage, our model seems capable of qualitatively 
accounting for the evolution of xenon isotopes 
in the martian interior and atmosphere provided 
degassing of the crust is limited to an early pe- 
riod of martian history. It produces 
,36 Xe*/ 129 Xe xs ratios low in the atmosphere and 
high in the crust and a source of gas close to 
solar in the mantle, as required by our present 
understanding of die components. This can be 
seen as a consequence of the concentration of 
incompatible elements in the crust during differ- 
entiation, supporting the argument for crustal 
assimilation. 

We plan further work on the basaltic sher- 
gottites to establish whether the variation be- 
tween Shergotty and EET79001 Lithology-B 
represents a wider trend, and hope to present 
data from DaG489 at the meeting. 

References: [1] Ocker, K. D. and Gil- 
mour, J. D., (2001) XXXII, CD#1782. [2] 
Ocker, K. D. and Gilmour, J. D., (2001) 36 (9 
abstracts), A 152. [ 3 ] Ocker, K. D. and Gil- 
mour, J. D., (2002) XXXIII, CD# 1603. [4] 
Mathew, K. J. and Marti, K., (2001) JGR 
106, 1401-1422. [ 5 ] Swindle, T. D. and 

Jones, J. H., (1997) J. Geophys. Res.-Planets 
102 (El), 1671-1678. [6] Gilmour, J. D., et 
al., (2001) GCA 65, 343-354. [ 7 ] Marty, B. 
and Marti, K., (2002) EPSL 196 (3-4), 251- 
263. [8] Pepin, R. O., (1991) Icarus 92, 2- 
79. [91 Pepin, R. 0., (1994) Icarus 111, 

289-304. [ 10 ] Podosek, F. A., (1973) EPSL 
19, 135-144. 
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Introduction: Noble gases in the solar wind and 
giant-planet atmospheres are generally assumed to be 
the best — perhaps the only — available proxies for 
isotopic distributions in the early solar nebula. Both, 
however, may be isotopically fractionated to some de- 
gree from their source composition, the wind in proc- 
esses that transport solar plasma to and release it from 
the corona [1], and giant-planet Ar, Kr and Xe by trap- 
ping in icy planetesimals if these were indeed the prin- 
cipal suppliers of heavy noble gases to their present 
atmospheres [2]. Noble gas isotopic ratios in the solar 
wind are reasonably well established from lunar and 
asteroidal regolith studies [3,4] and the Galileo Probe 
mass spectrometer has given us our first look at com- 
positions in Jupiter’s atmosphere [2]. 

Modem theories of atmospheric evolution on the 
terrestrial planets are focused on the nebula as the pri- 
mary supplier of primordial planetary volatiles, either 
directly — i.e., in ways that do not fractionate isotopes 
(by adsorption on planetary embryos or "ingassing" 
from dense, gravitationally condensed atmospheres) — 
or by accretion of cometary ices carrying noble gases 
that could have been either isotopically solar or mildly 
fractionated during trapping from the ambient nebula 
[5], The observation that nonradiogenic Ne, Ar, Kr and 
Xe in Earth's current atmosphere are all isotopically 
heavier than their solar counterparts is an important 
clue to the nature of the processes that subsequently 
acted on these primordial planetary reservoirs. This is 
also the case on Mars, except for Kr which is isotopi- 
cally near-solar (and therefore an interesting challenge 
for evolutionary modeling), and on Venus as well, as 
far as one can tell from the limited data on hand [1,5]. 

Much of the evolutionary modeling over the past 1 5 
years or so has focused on fractionation from primor- 
dial source compositions during thermally-driven, hy- 
drodynamic atmospheric escape (reviewed in [5]). Hy- 
drogen-rich primordial atmospheres on partially or 
fully accreted planets are heated at high altitudes by 
intense ultraviolet radiation from the young sun or al- 
ternatively, in the case of the Earth, by energy depos- 
ited in a large Moon-forming impact. Under these con- 
ditions hydrogen escape fluxes can be large enough to 
exert upward drag forces on heavier atmospheric con- 
stituents sufficient to lift them out of the atmosphere. 
Lighter species are entrained with the outflowing hy- 
drogen and lost more readily than heavier ones, leading 
to mass fractionation of the residual atmosphere. 

Models in which a planetary atmosphere is driven 
from isotopically solar-like progenitors to its present 
compositional state by energetic hydrodynamic escape 
of the primary atmosphere, followed by mixing of the 


isotopically fractionated atmospheric residue with solar 
gases evolved from the interior, are able to account for 
many of the details of contemporary noble gas isotopic 
distributions on Earth [5-7]. The one glaring exception 
is Xe. Isotope ratios of terrestrial nonradiogenic Xe 
cannot be generated by hydrodynamic escape 
fractionation of the measured solar-wind Xe composi- 
tion, An initial composition called U-Xe, which ap- 
pears to be isotopically identical to solar-wind-Xe ex- 
cept for sharply lower abundances of the two heaviest 
isotopes, is required [8]. There is some evidence that a 
U-Xe-like composition may have been present in the 
presolar molecular cloud and at some stage in the his- 
tory of the solar nebula. But if so, we are still left with 
the puzzling question of why it differs isotopically 
from the Xe composition in the solar wind and pre- 
sumably also in the sun. 

Applied to Venus, the hydrodynamic escape model 
makes specific predictions for the currently unknown 
Kr and Xe isotopic distributions in its contemporary 
atmosphere [7]. This opportunity for an observational 
test of the modeling assumptions is one of the primary 
scientific drivers in current proposals for a suitably 
instrumented Discovery mission to the planet. 

Mars: Mars, as usual, is interestingly different. 
Atmospheric Kr is isotopically similar to solar-wind 
Kr, as noted above, and so is evidently largely unfrac- 
tionated. Severe fractionation, however, is needed to 
establish present Xe isotope ratios, and here the com- 
position of its primordial progenitor is completely con- 
sistent with solar-wind Xe, not U-Xe. Current Ar and 
Ne isotopic compositions appear to result from inter- 
play of planetary degassing and fractionating sputtering 
loss from the exosphere over the past several Ga. 

Current modeling divides Martian atmospheric his- 
tory into early and late evolutionary periods (e.g. [9- 
14]). The first is characterized by an initial episode of 
hydrodynamic escape that sets the Martian Xe inven- 
tory, followed by high C0 2 pressures and a possible 
greenhouse, and the second by either a sudden transi- 
tion to a low pressure environment similar to present- 
day conditions on the planet — perhaps initiated by 
abrupt polar C0 2 condensation ~3.7 Ga ago [9] — or by 
a more uniform decline to present pressure. 

The early atmosphere . Early in the pre-3.7 Ga ep- 
och, nonradiogenic Xe isotopes are assumed to be hy- 
drodynamically fractionated to their present composi- 
tion, with corresponding depletions and fractionations 
of lighter primordial atmospheric species [11]. Subse- 
quent C0 2 pressure and isotopic history is dictated by 
the interplay of estimated losses to impact erosion, 
sputtering, and carbonate precipitation, additions by 
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outgassing and carbonate recycling, and perhaps also 
by feedback stabilization under greenhouse conditions 
[11,14]. It should be stressed, however, that little is 
actually known about the values of the parameters gov- 
erning these various processes. Models of this epoch 
are no more than qualitative illustrations of how they 
might have driven early atmospheric behavior. 

Later atmospheric evolution: Sputtering loss. Prob- 
able operation, on Mars, of a sputtering loss mecha- 
nism that results in fractionation of elements and iso- 
topes in the residual atmosphere was demonstrated by 
[15]. Oxygen atoms in the Martian exosphere, ionized 
by solar ultraviolet radiation and accelerated in the 
motional vxB electric field of the solar wind, can im- 
pact species near the top of the atmosphere with 
enough energy transfer to eject them from the planet's 
gravitational field. Fractionations are governed primar- 
ily by diffusive mass separation in the upper atmos- 
phere, which enhances the relative abundances of light 
to heavy species. Sputtering loss therefore preferen- 
tially removes lighter masses, leaving the atmosphere 
enriched in heavier constituents [5,10,1 1]. 

C0 2 pressure collapse on Mars near 3.7 Ga would 
have triggered the rapid sputtering removal of pre- 
existing Ne and Ar from the atmosphere (Kr and Xe 
are too heavy to be affected) [11]. This has the interest- 
ing consequence that no isotopic memory of their ear- 
lier processing survives. Current abundances and iso- 
topic compositions are entirely determined by the ac- 
tion of sputtering on gases supplied by outgassing dur- 
ing the late evolutionary epoch, and final distributions 
of the light noble gases are therefore decoupled from 
whatever their elemental and isotopic inventories might 
have been in the pre-3.7 Ga atmosphere. The issue of 
how much Ne and Ar degassing was required to pre- 
vent their complete sputtering annihilation from the 
atmosphere during this epoch is unresolved. One study 
concluded that the Ne-Ar outgassing flux attributable 
to degassing during major episodes of volcanic activity 
would have been ~l-3 orders of magnitude too low to 
balance sputtering losses, and thus another major 
source of juvenile volatiles must have contributed to 
the atmosphere over geologic time, perhaps via input 
from gas-enriched hydrothermal systems [12]. 

Sputtering is greatly attenuated by the presence of a 
planetary magnetic field [13]. Efficient operation of 
this fractionating loss mechanism over time on Mars is 
thus centrally linked to the timing of the disappearance 
of the Martian paleomagnetic field [13,16]. 

Atmospheric Xe and the 244 Pu problem. The iso- 
topic composition of Martian atmospheric Xe most 
recently derived from Xe trapped in the shock-glass 
phases of shergottites [17] is almost perfectly consis- 
tent (except for l29 Xe) with a composition generated 
entirely by fractionation of solar-wind Xe [8]. This has 
the interesting consequence that additional contribu- 
tions to the Xe inventory from degassed fission Xe are 
either very minor or absent altogether. And yet the pre- 


sent atmosphere is clearly heavily enriched in 129 Xe 
from extinct ,29 I decay. The presence of radiogenic 
129 Xe but apparently little if any 244 Pu fission Xe, a 
situation quite unlike that on Earth, has been awkward 
to reconcile with models of Martian geochemical evo- 
lution and degassing history. A less severe fractiona- 
tion of SW-Xe does allow the presence of an additional 
component with Pu-Xe fission yield ratios, amounting 
to ~5% of total 136 Xe [18] — similar to the fraction on 
Earth. At the moment the central question of the pres- 
ence or absence of Pu-Xe is plagued by apparent coin- 
cidences. The excellent match of fractionated SW-Xe 
alone to the atmospheric data is presumably fortuitous 
if Pu-Xe is present; and if it is absent, the fact that a 
weaker fractionation generates a pattern in good accord 
with 244 Pu fission yields must likewise be accidental. 

Interior Xe. The nonradiogenic Xe isotope ratios 
found in Chassigny appear to be indistinguishable from 
solar values [19]. The relative abundances of I29 Xe and 
136 Xe are also close to solar, and so there is little scope 
for radiogenic additions. Data from Nakhla [20] indi- 
cate that there are other interior Martian reservoirs that 
contain solar Xe but with resolvable fissiogenic contri- 
butions, and so have had lower Xe/Pu ratios. It seems 
likely that the interior is relatively undegassed, with 
high ratios of primordial Xe to parent elements. 

Degassing history. Viking measurements of atmos- 
pheric 40 Ar and estimates of 40 K in the mantle from the 
SNCs suggest that only ~2% of the 40 Ar generated in 
the planet has degassed. Atmospheric 129 Xe from ex- 
tinct 129 I decay is consistent with ~2% degassing if the 
closure age of Mars against gas loss during accretion 
was ~70 Ma; the same closure age and degassed frac- 
tion would be indicated by the atmospheric 136 Xe in- 
ventory if ~5% of it derives from 244 Pu decay [5] — 
currently an open issue, as noted above. 

References: [1] Wieler R. (2002) in Noble Gases in 
Geochemistry and Cosmochemistry, Rev. Mineral Geochem. 
47, pp. 20-. [2] Mahafly P.R. et al. (2000) JGR 105, 15061- 
15071. [3] Palma R.L. et al. (2002) GCA 66, 2929-2958. [4] 
Pepin R.O. et al. (1995) GCA 59, 4997-5022. [5] Pepin R.O. 
and Porcelli D. (2002), in Noble Gases in Geochemistry and 
Cosmochemistry, Rev. Mineral Geochem. 47, pp. 191-246. 
[6] Pepin R.O. (1991 ) Icarus 92, 2-79. [7] Pepin R.O. (1997) 
Icarus 126, 148-156. [8] Pepin R.O. (2000) Space Sci. Rev. 
92, 371-395. [9] Haberle R.M. et al. (1994) Icarus 109, 102- 
120. [10] Jakosky B.M. et al. (1994) Icarus 111, 271-288. 
[11] Pepin R.O. (1994 ) Icarus 111, 289-304. [12] Hutchins 
K.S. and Jakosky B.M. (1996) JGR 101, 14933-14949. [13] 
Hutchins K.S. et al. (1997) JGR 102, 9183-9189. [14] Carr 
M.H. (1999) JGR 104, 21897-21909. [15] Luhmann J.G. et 
al. (1992) GRL 19, 2151-2154. [16] Connemey J.E.P. et al. 
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Introduction: Martian meteorite (shergottite) in- 
pact melt glasses that contain high concentrations of 
martian atmospheric noble gases [1] and show signifi- 
cant variations in 87 Sr/ 86 Sr istopic ratios [2] are likely to 
contain Martian surface fines mixed with coarser e- 
golith materials. The mixed soil constituents were mo 1- 
ten due to shock at the time of meteoroid impact near 
the Martian surface and the molten glass got incorpo- 
rated into the voids and cracks in some shergottite 
meteorites. Earlier, Rao et al. [3] found large anich- 
ments of sulfur (sulfate) during an electron-microprobe 
study of several impact melt glass veins and pods in 
EET79001, LithC thin sections. As sulfur is very abun- 
dant in Martian soil, these S excesses were attributed to 
the mixing of a soil component containg aqueously- 
altered secondary minerals with the LithC precursor 
materials prior to impact melt generation. Recently, we 
studied additional impact melt glasses in two basaltic 
shergottites, Zagami and Shergotty using procedures 
similar to those described in [3]. Significant S enrich- 
ments in Zagami and Shergotty impact melt glass veins 
similar to the EET79001, LithC glasses were found. In 
addition, we noticed the depletion of the mafic comp o- 
nent accompanied by the enrichment of felsic compo- 
nent in these impact melt glass veins relative to the 
bulk host rock in the shergottites. To explain these ob- 
servations, we present a model based on comminution 
of basaltic rocks due to meteroid bombardment on mar- 
tian regolith and mechanical fractionation leading to 
enrichemnt of felsics and depletion of mafics in the fine 
grained dust which is locally mobilized as a result of 
saltation and deflation due to the pervasive aeolian 
activty on Mars. 

Methods and Results: Several impact melt glass 
veins and pods in EET79001, LithC (,77 ,78 ,18, 20A and 
,506) sections ; Zagami (,992 a2 & c2) and Shergotty 
(dbs) glasses were studied using a Cameca SX100 rri- 
crobeam automated electron microprobe equipped with 
standard Cameca instrument controls and PAP data 
reduction software. Standard analytical methods and 
calibration procedures were used. The elemental abun- 
dance data obtained by raster scans along different line 
traverses at homogeneous sites in the veins and pods 
are averaged and normalized to the bulk composition of 
the host rock. Typical results (average) for 
EET79001, LithC are shown in Fig 1. In these shergottite 
glass veins, the maximum enrichment of A1 2 0 3 is upto 
~60%; CaO is upto ~20%; and Na 2 0, upto ~60%. On 
the otherhand, the maximum depletion of FeO is upto 


~15%; MgO is upto ~36% and Ti0 2 upto ~36%. Also, 
the volatile element, S showed large S0 3 excesses (upto 
a maximum of ~300%) in these shergottite glasses over 
the host rock. Further, sulfur was shown to be present 
predominantly as sulfate (oxidized S / reduced S = 1 .83) 
in LithC glass veins by [4]. However, the Cl abundance 
in the LithC glasses is very low (<0.05%) compared to 
the Cl abundance estimated from the CI/S ratio of the 
Viking and Pathfinder soils [5,6]. Note that precautions 
were taken to avoid heating losses during measurement 
of alkalies and halogens by using defocussed beam, 
low current and long counting rates. The results sug- 
gest that the impact melt glass veins in shergottites are 
enriched in felsics (Al & Ca) and depleted in mafics (Fe 
& Mg) compared to the host rock. Further, S which is 
an index element for Martian soil is enriched in these 
glass veins, whereas Cl (also a volatile) is depleted in 
these veins. These results suggest that the martian soil 
in the shergottite glass veins is different from the Vi- 
king and Pathfinder soils. 

Discussion and the Model: During the study of 
glass veins in Zagami, [7] noted “slight enrichments of 
Si0 2 , AI2O3 and Na 2 0 along with depletions of FeO and 
MgO relative to the bulk”. Similar observations were 
made by [8] in the case of EET79001, LithC glasses. 
Further, [9] noted large enrichment of LREE in 
EET79001, LithC glass during an ion-microprobe study. 

Felsic enrichment and mafic depletion. Modest- 
sized rocks in Martian regolith undergo collissional 
destruction by multiple impacts due to meteoroid bom- 
bardment. It was shown by [10] and [11] that small- 
sized meteoroids survive passage through martian at- 
mosphere and encounter the surface with few km/sec 
velocity. As a result, the unconsolidated materials will 
be incresingly fractured and shock-processed with time 
into fine-grained soil. This bombardment causes exten- 
sive regolith production and redistribution of dust and 
rubble on Mars surface[12]. From the impact-gardened 
regolith on Mars, finer material is mobilized and redis- 
tributed locally by deflation and saltation processes 
due to aeolian activity. 

Hoerz and Cintala [13] carried out a series of lab ex- 
periments using the high velocity impactor machine at 
JSC to simulate the conditions of bombardment of ba- 
saltic roks in space using gabbro target. The results 
indicate that in the fine-grained comminution products, 
plagioclase (Al and Ca) is enriched as it comminutes 
easily whereas mafic materials (Fe and Mg) which are 
hard are depleted relative to the starting materials. This 
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mechanical fractionation was found to be most severe 
in fine grain sizes (<20 urn). Further, in thencase of lu- 
nar soils, [14] showed that the <10 um fraction of lunar 
soils is enriched in felsics (A1 2 0 3 , CaO, Na 2 0, K 2 0 and 
LREE) and depleted in mafics ( FeO, MgO, MnO) rela- 
tive to the bulk lunar soil. These results are consistent 
with those obtained for shergottite glasses, indicating 
that their precursor materials were similarly produced 
on Mars surface. 

S and Cl in shergottite glasses: Secondary sulfate 
and chloride ninerals in Martian soil result from the 
volatiles entring into the atmosphere through volcanic 
emnations and subsequently interacting with aqueous 
media on Mars surface. S is enriched in EET79001, Lith 
C glasses compared to Lith A. Elemental correlation 
diagrams indicate that A1 2 0 3 and CaO correlate posi- 
tively with S0 3 whereas MgO correlates negatively 
with S0 3 . In contrast, the Cl abumdance in LithC sam- 
ples is very low. Halite was never found during SEM 
studies in these glasses. The S/Cl ratio in the shergot- 
tite glasses is >50 whereas the same ratio in Viking and 
Pathfinder soils is 5. These results suggest that in- 
soluble A1 & Ca rich sulfates are decoupled and mobi- 
lized from soluble sulfates of Mg and chlorides of Na & 
K by local aqueous activity near the source regions of 
the precursorr materials on Mars regolith. 

Conclusion: The shergottite impact glasses contain 
fine-grained soil, which was exposed to local fluid 
activity in the subsurface regolith associated with 
young volcanic upland regions (Olympus Mons or 
Elysium) on Mars. 

References: [1] Bogard D.D. and Johnson P. (1983) 
Science, 221 , 651-654. [2] Nyquist L.E. et al. (1986) 
LPSC XVII , 624-625 (abstract). [3] Rao M.N. et al. 
(1999) GRL, 26, 3265-3268. [4] Gooding J.L. and Mue- 
now D.W. (1986) GCA, 50, 1049-1059. [5] Clark B.C. et 
al. (1982) JGR, 87 , 10059-10067. [6] Reider R. et al. 
(1998) Science , 278 , 1771-1774. [7] McCoy T.J. et al. 
(1992) GCA, 56 , 3571-3582. [8] McSween H.Y. Jr. and 
Jarosewich E. (1983) GCA, 47, 1501-1513. [9] Shearer 
C.K. et al. (2000) First Astrobiology Conf Ames . [10). 
Hoerz F. et al. (1999) Science, 285 , 2105-2107. [11] 
Hartmann W. K. et al. (2001) Icarus, 149 , 37-53. [12] 
Hoerz F. and Cintala M. (1997) MAPS, 32, 179-209. [13] 
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761-826. 
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Several quantitative studies on thermal stories of SNC meteorites have been carried out to 
infer burial depth. Two approaches are generally used: (i) indirect estimation based on 
textural observation compared with products from controlled cooling crystallization 
experiment (ii) direct quantitative cooling rate calculation from diffusion zoning in mineral. 
The first method, applied on basaltic shergottite, is largely based on textures and compositions 
of pyroxene and/ or maskelinite. The second one, performed on picritic shergotite (EET79001 
and Dar al Gani 476) and the nakhlite NWA 817, uses zoning of olivine to extract cooling 
rate. At the onset of cooling, olivine was homogeneous and was then modified by atomic 
diffusion due to contact with Fe-rich surrounding melt when olivine was incorporated as a 
xenocryst. 

It had long been argued that olivine speedometer, giving absolute cooling rate, was the 
reliable method compared to the first one which bears large uncertainties. Unfortunately any 
contribution using olivine as a speedometer in meteorites are negligent (or at least optimistic) 
about the assignment of uncertainties to cooling rate estimate. In such studies numerical 
procedures, where diffusion equation is solved by finite difference approximation, are 
generally preferred to analytical one. However error propagation becomes obvious when 
using analytical methods in as much they allow the dependence of cooling rate results on 
input parameter to be clearly seen. For instance temperature To and diffusion data (E the 
activation energy and Do the preexponential term of diffusion coefficient D expressed 
through the Arrhenius law) are first order parameters in cooling rate equation. Small errors on 
temperature and diffusion data are magnified due to the logarithmic relation ship between 
input T and cooling rate S. 

Reading the literature on SNC meteorites, uncertainties on To, temperature at which cooling 
started, are at least of ± 50°C. In NWA 817, for a given D (Buning & Buseck 74 for T< 1 125 
°C), S is 0.41 K/h for To = 1 100°C and 0.043 K/h for To = 1200°C (B & B 74 T>1 125°C). 
Uncertainties on input diffusion parameters are worst compared to temperature one, as 
diffusivity of a given element (Fe-Mg) in a given mineral (olivine) is highly variable due to 
different laboratory experimental conditions and crystal impurities. For sake of demonstration 
let's used two extreme values of E>Fe-Mg given respectively by B & B 74 and Chakraborty 97. 
In NWA 817, for a given To of 1 100°C, B&B 74 diffusion coefficient gives S of 0.41 K/h 
and C.97 's one gives 2.66 10‘ 3 K/h i.e. 2 orders of magnitude difference in the cooling rate 
estimate. 

To conclude, cooling rate quantification of SNC meteorites using olivine speedometer is a 
perilous exercise as such a method bears itself large uncertainties. Those results should be 
cross-calibrated with data obtained from controlled cooling crystallization experiments. 
However it should be noted that fast crystallization rate does not always relate to fast cooling 
and may be rather due to supersaturation phenomenon. 
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Introduction: Isotopic compositions of nitrogen 

and noble gases have played a key role in understand- 
ing SNC meteorites and Mars. For example nitrogen and 
argon data from the glass of EET 79001 were one of the 
convincing evidences that linked these meteorites to 
Mars [1]. Similarly, the observation of comp onents with 
largely different isotopic compositions (of both nitro- 
gen [2, 3] and noble gases [4]) in these meteorites has 
led to the realization that the martian atmosphere (as 
measured by Viking) is quite different from the martian 
interior (as pictured from data from these meteorites). 
But at the same time we have not yet been able to 
clearly define the (interior) mantle end member and re- 
duce the uncertainties presented by the Viking data for 
the martian atmosphere. Similarly we do not know for 
sure the 20 Ne/ 22 Ne signature of Mars. It is clear that 
many of these ambiguities as to the indigenous martian 
components result from interferences from cosmogenic 
contributions and terrestrial contamination, which are 
impossible to avoid. But they can be minimized by step 
temperature heating, and selection of mineral phases 
rather than bulk samples. We have initiated a simulta- 
neous nitrogen and noble gas isotopic study of mineral 
separates and a few bulk samples from martian meteor- 
ites with the newly set up experimental facility at Mainz 
[5]. Here we report data for maskelynite and pyroxene 
separates from Zagami. 

Experimental: Gases were extracted by vacuum py- 
rolysis in an iridium crucible in temperature steps, while 
simultaneously adsorbing them on a stainless-steel- 
mesh (SSM) finger (at -196 °C) near the crucible. A split 
of ~5 % of the sample gas was processed for the analy- 
sis of nitrogen, while the rest was processed for the 
noble gases using standard procedures [6]. The nitro- 
gen fraction, prior to analysis, was cleaned by oxygen 
generated by an in-line copper oxide finger in a line 
made of quartz and pyrex. Pipettes of Air standard were 
run similar to the sample steps to calibrate the experi- 
ment, which were consistent with a sensitivity of 4xl0' 7 


Table 1 . Nitrogen and selected noble gas data. 


N 

a 15 N 

“Ar/^Ar 

129 Xe/ 132 Xe 

(ppm) 

(%o) 



Mask. 25 

23 

3600 

1.20 

Px. 10 

13 

2160 

1.09 


|Xg/mV (±10 %). The nitrogen isotopic composition 
(a 15 N) of the standards was up to 2 %o heavier than that 
of air. A typical blank for nitrogen at 1600 °C was 60 
picograms, with a a 15 N of +2 %o. The noble gas blanks 


and calibration parameters were similar to those of [6]. 
The nitrogen data have additionally been corrected for 
background contributions from CO using the procedure 

of Pi- 

Discussion: The concentrations of nitrogen and 
a ,5 N in the maskelynite sample are about 2 times higher 
as compared to that in pyroxene (Table 1). Although 
one expects the low temperature extraction (at 400 °C) 
to have significant contribution from terrestrial con- 
tamination (as seen in the 40 Ar/ 36 Ar and 129 Xe/ U2 Xe 
ratios), the a 15 N (+22 %o in maskelynite and +13 %o in 
pyroxene) are much different from that of terrestrial 
atmosphere. The high temperature release (accounting 
for ~85 % of nitrogen) in maskelynite shows a variation 
from +15 to +25 %o in a 15 N (Fig. 1). Nitrogen in pyroxene 
is dominated by the release (accounting for -70 %) at 
1200 °C which has a a ,5 N of +10 %o. a l5 N of the 800 and 
1600 °C steps of this sample are similar to those of 
maskelynite. Marti et al. [8] have shown the presence of 
a martian atmospheric component (for ritrogen and 
noble gases) in glass samples of Zagami. A heavy a 15 N 
of +147 %o observed by Mathew et al. [9] in a bulk 
sample of Zagami also suggests the presence of a 
martian atmospheric component in this meteorite. While 
the present nitrogen concentrations are at least an 
order of magnitude higher than those in the Zagami 
glass [8], the a 15 N show significant contribution from a 
component with lighter a 15 N. 



I m Mask. (14 ms) 

0 4 —t ■ ■ i *-r : * — •H 

0 20 40 60 60 100 


% Release of Nitrogen j 

Fig. 1 . Stepped temperature release of nitrogen. 

Similar to nitrogen, concentrations of trapped Ar, Kr 
and Xe in maskelynite are a factor of two higher as 
compared to those in pyroxene. A considerable fraction 
(Figs. 2 and 3) of these gases released in the 400 °C 
step, could have been contributed by terrestrial con- 
tamination, as suggested by the 40 Ar/ 36 Ar and 
129 Xe/ 132 Xe signatures. But the gases released in the 
high temperature (>400 °C) extractions are clearly dif- 
ferent. Neon and argon in the present samples are 



Workshop on linmixing the SNCs (2002) 


6010.pdf 


N & NOBLE GASES IN MINERAL SEPARATES FROM ZAGAMI: S. P. Schwenzer et al. 


dominated by the cosmogenic component (up to 90 % 
for Ne and 79 % for Ar). We have corrected the nitro- 
gen data for cosmogenic contributions using the 
21 Ne cosm following [10]. The corrections amount to less 
than 1 %o. Similarly, the measured 40 Ar/ 36 Ar ratios have 
been corrected for the cosmogenic contributions fol- 
lowing [4]. 



% Release of M Ar 


Fig. 2. Stepped temperature release of trapped argon. 

Although 40 'At/ 36 Ar of the low temperature extrac- 
tion in these samples shows a contribution from terres- 
trial atmospheric argon, that in the high temperatures 
shows the presence of contributions from an end- 
member similar to martian atmosphere (as shown by the 
release in 800 °C in pyroxene) and radiogenic argon (as 
shown by the elevated signatures in 1200 and 1600 °C 
releases of maskelynite). The higher contribution from 
the radiogenic component in maskelynite compared to 
pyroxene can be explained by the higher K content of 
maskelynite. 40 Ar/ 36 Ar in the 1200 and 1600 °C extrac- 
tions of pyroxene, contrary to maskelynite, indicates 
contributions from a component with low 40 Ar/ 36 Ar 
(<1000), which can be explained by the martian mantle 
(e.g., [11]). I29 Xe/ 132 Xe in the present samples varies 
from an air-like signature (in the 400 °C extractions) to 
1 .3 (in 1 600 °C extraction of maskelynite), and is smaller 
than the signatures observed in the Zagami glass [8]. 



0 20 40 60 SO 100 

% Release of 132 Xe 


Fig. 3. Stepped temperature release of trapped xenon. 

The higher gas (apparently from martian atmos- 
phere) concentrations in maskelynite as comapraed to 
pyroxene may be explained by the fact that the former is 


a result of shock metamosphism which has been shown 
to be efficient in implanting ambient atmo spheric gases 
into shock-produced phases [12]. However, it should 
be noted that maskelynite (+mesostasis) and pyroxene 
from Shergotty [11] do not show any significant differ- 
ences in their trapped gas concentrations (except for 
He). 

References: [1] BeckeT R. H. and Pepin R. O. (1984) 
Earth Planet Sci. Lett, 69, 225-242. [2] Mohapatra R. 
K. et al. (1998) Meteoritics & Planet Set, 33, Al 12. [3] 
Mathew K. J. and Marti K. (2001) JGR, 106, 1401-1422. 
[4] OttU. (1988) Geochim. Cosmochim. Acta, 52, 1937- 
1948. [5] Mohapatra R. K. et al. (2001) Meteoritics & 
Planet . Sci., 36, A139. [6] Schelhaas N. et al. (1990) 
Geochim. Cosmochim. Acta, 54, 2869-2882. [7] Murty 
S.V.S. (1997) Meteoritics & Planet Set, 32, 687-691. [8] 
Marti K. et al. (1995) Science, 267, 1981-1984. [9] 
Mathew K. J. et al. (1998) Meteoritics & Planet . Sci., 
33, 655-664. [10] Mathew K. J. and Murty S. V. S. (1993) 
Proc. Indian Acad. Sci. Earth Planet Sci., 102, 415- 
437. [11] Terribilini et al. (1998) Meteoritics & Planet. 
Sci., 33, 677-684. [12] Wiens R. C. and Pepin R. Q 
(1988) Geochim. Cosmochim. Acta, 52, 295-307. 





Workshop on Unmixing the SNCs (2002) 


6017.pdf 


CHEMICAL COMPOSITION OF NEWLY COLLECTED ANTARCTIC NAKHLITES, Y000593 AND Y000749. N 

Shirai, Y. OURA and M. Ebihara, Department of Chemistry, Tokyo Metropolitan University, hachioji, Tokyo 192-0397, 
Japan (correspondence to ebihara-mitsum@c.metro-u.ac.ipT 


Introduction: Y000593 and a possibly paired mete- 
orite Y000749are grouped into nakhlites, a group of SNC 
meteorites, by their noble gas compositions [1] They are 
the first nakhlite recovered on Antarctica. So far about 30 
SNC meteorites (shergottite, nakhlite and chassignite) 
have been identified, most of which are classified into 
shergottites. Nakhlites are rare, to which 6 meteorites 
including Y000593 and Y000749 belong. In this study, we 
have analyzed Y000593, Y000749 and Nakhla for their 
bulk chemical compositions by nuclear analytical 
methods. At the Antarctic meteorite symp osium at NIPR 
of this year, we 1 reported some preliminary results [2]. 
We have being continuing a suit of analyses and their 
results are presented in this woikshop. 

Experimental: We first performed neutron-induced 
prompt gamma-ray analysis (PGA). For PGA, lump 
samples of Y000593,66 (1.2 g), Y000749,46 (3.6 g) and 
Nakhla (0.88 g) were used. PGA was conducted by using 
cold neutron beam at JRR-3M, Japan Atomic Energy 
Research Institute. The specimens of Y000593 and 
Y000749 were loaned from NIPR and were subjected to 
PGA without any additional sample preparation at our 
side. Because the size of the Y000749 specimen is larger 
than the neutron beam size used in PGA, only content 
ratios normalized to Si content were determined for 
Y000749 and the other samples by the internal 
mono-standard method developed by our group. Ab- 
solute contents were also obtained for Y000593 and 
Nakhla by conventional comparison method using JB-1 
(a basaltic geological reference sample issued by GSJ) 
and chemicals of analytical grade as reference standards 
for quantification. For instrumental neutron activation 
analysis (INAA) and instrumental photon activation 
analysis (IPAA), powdered samples were prepared for 
Y000593 and Nakhla. Six chunks weighing about 2 g each, 
which are located at different position from each other, 
were sampled from the interior of Y000593 and powdered 
at NIPR. About 100 mg of each powdered sample was 
used for INAA and IPAA. 

Results and Discussion: Since Si-normalized con- 
tents for major elements obtained by PGA using the 
internal mono-standard method are almost the same for 
three nakhlites analyzed, it is inferred that their absolute 
elemental abundances are also similar Our absolute 
content values for Y000593 and Nakhla obtained by PGA 
are shown in Table 1 with contents in Nakhla compiled by 
Lodders [3]. All major elements and several minor and 
trace elements can be determined by PGA . It is seen in 
Table 1 thatthere is no apparent difference in contents of 


analyzed elements between Y000593 and Nakhla except 
forH and Cl andthat our PGA results for Nakhla are also 
consistent with compilation data [3]. There appear clear 
boundaries in Mg/Si, Al/Si, Ca/Si, and/or Ti/Si weight 
ratios among SNC groups; for example, Mg/Si - Ca/Si 
diagram is shown in Fig.l . These ratios for Y000593 and 
Y000749 are different from those for Shergotty and 
Chassigny but consistent with those for Nakhla. 
Therefore it isconfirmed based on chemical composition 
that Y000593 is a nakhlite. Our determination of Cl in 
Nahklais about 16 times higher thanthe compilation data. 
As our sample of Nakhla was sampled from the interior of 
a large mass at the lunar sample processing room at JSC, 
NASA and no additional treatment was given for PGA, a 
high content of Cl cannot be attributed to experimental 
contamination. It must be highly probable that Cl is 
heterogeneously distributed in nakhlite, just like the case 
for chondrites [4]. As only one value has been reported 
for Cl content in Nakhla, additional analyses are desired 
for Cl (and other halogens) in future. INAA and IPAA 
data are under examination, and their results and dis- 
cussion will be presented at the woikshop. 


References: [1] Imae N. et al. (2002) LPS XXXIII , 
#1483 (CD-ROM). [2] Oura Y. et al. (2002) Antarctic 
Meteorites XXVII, 143-145 (abstract). [3] Lodders K. 
(1998) Meteoritics Planet Set 33, A183-A190. [4] 
Zolensky M. E. et al. (1999) Science 285, 1377-1379. 
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Chemical composition of Y000593 and Y000749: N. Shirai et al. 


Table 1 . PGA results for Y000593 and Nakhla* 


Element 

Y000593 

Nakhla 

Nakhla $ 

H # , ppm 

541 (14) 

252(10) 


B, ppm 

3.47 (6) 

3.02 (7) 

4.6 

Na, % 

0.432 (3) 

0.447 (25) 

0.34 (5) 

Mg, % 

6.26 (26) 

6.93 (29) 

7.3 (0.2) 

Al, % 

0.999 (30) 

0.974 (33) 

0.89 (11) 

Si, % 

22.2 (4) 

23.5 (4) 

22. 7 (8) 

Cl, ppm 

52.9 (6.1) 

872(13) 

80 

K, % 

0.119(7) 

0.117(7) 

0.107 (19) 

Ca, % 

10.2 (0.3) 

10.8(0.5) 

10.5(0.5) 

Ti, % 

0.218(4) 

0.208 (4) 

0.202 (25) 

Cr, % 

0,179 (6) 

0.208 (7) 

0.177(28) 

Mn, % 

0.398 (27) 

0.388 (30) 

0.382 (31) 

Fe, % 

15.3 (5) 

14.5 (4) 

16.0(1.2) 

Co, ppm 

91(11) 

79 (12) 

48(5) 

Ni,ppm 

179 (48) 

191 (55) 

90 

Sm, ppm 

1.46 (9) 

1.01 (7) 

0.77(8) 

Gd, ppm 

1.17(9) 

1.08 (7) 

0.86(8) 


♦Values in parentheses are errors due to errors (la), 
These values correspond to last digit(s) of the preceding 
figures. 

information values. 

S [3]. 
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Figure 1. Elemental abundance ratios of Mg/Si vs. Ca/Si 
for SNC meteorites. 
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dle@u. arizona.edu) 


Introduction: Nakhlites, the intermediate-aged 

(-1300 Ma) Martian meteorites, are similar to one an- 
other in noble gas characteristics, as well as mineral- 
ogy, but distinct from other Martian meteorites. There 
are also are subtle differences among the well-studied 
nakhlites (Nakhla, Lafayette, and Govemador Val- 
adares), which may be important in understanding the 
processes affecting Martian volatiles. This abstract will 
summarize some noble gases properties of Martian me- 
teorites, particularly nakhlites, and then highlight five 
unanswered questions. 

To first order, Xe in these meteorites appears to be 
a mixture of an interior component, seen mostly clearly 
in Chassigny, and an atmospheric component seen 
mostly clearly in shergottites like EETA 79001 [1, 2]. 
Nakhlites are rich in the atmospheric component (Nak- 
hla being the richest), but the component appears to be 
elementally fractionated, with a higher Xe/Kr ratio than 
the (atmospheric) component seen in EETA 79001. 
Nakhlite Ar has been studied most in the context of Ar- 
Ar age experiments. All three appear to be about 1300 
Ma old, but with some disturbances to the system, 
most noticeable in Nakhla. Shergottites all appear to 
have either incorporated variable amounts of two dif- 
ferent Ar components (interior and atmosphere?), to 
have incorporated different components from one an- 
other, or to have incompletely degassed upon forma- 
tion. 

Question 1: How did atmospheric (?) Xe get 
trapped in the nakhlites? This topic has been inten- 
sively investigated, and several different scenarios 
have been suggested, each of which has advantages 
and drawbacks. The basic problem is that while most of 
the Xe is isotopically indistinguishable from the Mar- 
tian atmosphere [3], the Kr/Xe elemental ratio is very 
distinct: it appears that atmospheric gas has been frac- 
tionated in favor of the heavier gas. ALH84001 shows a 
very similar effect although, unlike Nakhla, there may 
be some isotopic differences [4]. 

a. Incorporation related to liquid water : Liquid 
water will produce elemental fractionation of the right 
magnitude and direction. Since the nakhlites contain 
the strongest evidence for interaction with liquid water 
[5], Drake et al. [6] suggested that nakhlites contain 
atmospheric gas that was elementally fractionated by 
water. A similar scenario was suggested for ALH84001 
[7, 8]. In this scenario, the gas would be in the aqueous 
alteration products, iddingsite for the nakhlites, car- 
bonate for ALH84001. One problem with this scenario 


is that Nakhla has more Xe than Lafayette, although 
Lafayette has been more heavily altered. Most seri- 
ously, two sets of experiments on mineral separates 
have now shown that the bulk of the gas in the 
nakhlites is in the major minerals, not alteration prod- 
ucts [9, 10]. The gas in ALH84001 is also in pyroxene, 
not carbonate [11]. 

b. Magmatic Incorporation : Gilmour et al. [9] pro- 
posed that the gas was incorporated during crystalliza- 
tion, perhaps from soil degassing as the rock crystal- 
lized. It is not clear whether the outgassed material 
could be incorporated that efficiently. More recently, 
the same group suggested the next scenario [12]. 

c. Shock implantation of adsorbed gas : Adsorption 
will also produce elemental fractionation in favor of the 
heavy gas, but little or no isotopic fractionation. If 
shock can implant gas without further fractionation, 
this process could work. Gilmour et al. originally sug- 
gested this for ALH84001 [11], and later suggested it 
for nakhlites as well [12]. This is probably the leading 
candidate at the moment, but it has some problems that 
will be discussed below. 

d. Changing atmospheric elemental composition: 
Musselwhite and Swindle [13] pointed out that Xe 
could be preferentially incorporated into polar clath- 
rates under certain climate conditions, which would 
lead to an elemental fractionation in the atmosphere, 
which could have been implanted without any elemen- 
tal fractionation (as appears to have happened with the 
shergottites). Then, the differences between ALH84001 
and the nakhlites, on the one hand, and the shergot- 
tites on the other could be because of different climate 
conditions. Since ALH84001 and the nakhlites were all 
ejected from Mars >10 Ma ago, and the shergottites 
were all ejected <5 Ma ago, this is at least plausible. 

e. Recent addition of a cometary component to the 
Martian atmosphere: Owen and Bar-Nun [14] also 
suggested that the Martian atmosphere may suffer ele- 
mental changes over millions of years, but suggested 
that it could be the result of input of cometary gas. 

Question 2: Why did the nakhlites and ALH84001 
trap elementally fractionated atmosphere, when sher- 
gottites and Chassigny dich’t? Attributing the effect to 
aqueous alteration, would have explained this naturally 
for the nakhlites, but has been ruled out by further ex- 
periments. Invoking shock mplantation of adsorbed 
atmosphere is logical, but any shock scenario requires 
some special pleading. 
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The nakhlites have experienced the least shock of 
any of the Martian meteorites [e.g,, Table I in [15]], yet 
have abundant elementally fractionated atmosphere. 
One possibility is that different amounts of gas were 
adsorbed, presumably because of temperature. If the 
nakhlites and ALH84001 were ejected when the atmos- 
phere was colder (because of season, latitude, altitude, 
or some combination), they might have had more ad- 
sorbed atmosphere available for incorporation. Mean- 
while, the shergottites experienced higher shock pres- 
sures, high enough to produce melt glass that incorpo- 
rated unfractionated atmosphere. But what about 
Chassigny, which very possibly was ejected in the 
same event as the nakhlites [15]? Also, EETA79001, 
ALH84001, and Chassigny experienced comparable 
shock pressures [15], but the first incorporated copious 
quantities of un fractionated atmosphere, the second 
incorporated elementally fractionated atmosphere, and 
the third incorporated no atmosphere at all. Other dif- 
ferences (e.g., the presence of vesicles or the time of 
incorporation) could explain at least some of the differ- 
ence, but could such differences explain everything? 

Question 3: Why did Chassigny incorporate so 
much more interior Kr and Xe than any other Martian 
meteorite? By far the best sample of Kr and Xe from a 
component that appears to represent the Martian inte- 
rior comes from Chassigny [1,2, 4], although this com- 
ponent is apparently present in Nakhla as well [3]. On 
Earth, the best samples of mantle noble gases come 
from quickly cooled undersea basalts [16]. Although it 
is a dunite, a reasonable choice for a mantle -derived 
rock, Chassigny is definitely not a chilled basalt. 

Question 4: Why do shergottites show signs of in- 
terior-derived Ar, when the nakhlites and Chassigny 
do not? Shergottites are notorious for having one or 
more trapped argon components that make it virtually 
impossible to perform 40 Ar- 39 Ar dating [17]. This would 
be relatively easy to understand if it were merely be- 
cause of the presence of atmospheric Ar associated 
with the elemental unfractionated (i.e., high Ar/Xe) at- 
mosphere trapped in impact glasses. However, a simple 
correction for Martian atmosphere still does not yield 
interpretable Esults [17]. Rather, the “trapped” argon 
seems to be a mixture of Martian amosphere and a 
component with a lower 40 Ar/ 36 Ar Btio, presumably 
from the crust or mantle [17, 18]. Chassigny and the 
nakhlites, on the other hand, seem to have much sim- 
pler Ar-Ar systems [17, 19, 20]. Yet in Xe, the shergot- 
tites seem to have less evidence for any interior com- 
ponent. An obvious source for the interior-derived Ar 
in shergottites is crustal material that may have been 
assimilated. However, shergottite QUE94201 has one of 
the lowest 40 Ar/ 36 Ar ratios [18] (implying the most 


crustal material), but it seems to have assimilated the 
least material, based on other criteria 

Question 5: Why does the Ar-Ar pattern in Nakhla 
the most disturbed of the nakhlites, when Lafayette 
has the most evidence for aqueous alteration? Both 
Lafayette and Govemador Valadares have relatively 
simple Ar-Ar patterns, suggesting crystallization 
~1300Ma ago, with only a little evidence for loss of 
radiogenic 40 Ar since [19, 20]. Nakhla shows much more 
evidence of later disturbance than Lafayette or Gover- 
nador Valadares [19, 20], and also gives an Ar-Ar age 
that is the same or older, even though it gives slightly 
younger ages in other systems [15]. A quick answer to 
part of the question is that the products of aqueous 
alteration must not be affecting the Ar-Ar systematics 
much. While that must be true, it still doesn’t explain 
what the source of disturbance in Nakhla is. The key to 
the answer may be that while the nakhlites are all simi- 
lar, they are not identical. For example, Nakhla’s Sm-Nd 
systematic is also distinct from the other two [21, 22]. 
Remember that in the Kr/Xe system, Nakhla is again 
slightly, but distinctly, different as well. The distinction 
in histories may be sufficient to explain the noble 
gases. 

To solve these noble gas puzzles, we may ultimately 
have to understand more details about the magmatic, 
eruption, and even climate history of these meteorites 
than we currently do. 
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Correlations of elements abundances among the 
Martian meteorites have been used to infer the bulk 
composition of the Martian mantle [1-3], and thence the 
broad petrogenesis of the planet. It has been assumed 
that the correlations represent igneous fractionations. 
However, the correlations are also consistent with mix- 
ing, broadly similar to that inferred from radio-isotope 
data. If the correlations represent mixing of geochemical 
components, it is not clear that they represent Mars’ 
mantle or bulk composition. 

Introduction: Most estimates of the chemical com- 
position of Mars’ mantle and core derive from element- 
element correlations bulk analyses of Martian meteor- 
ites [1-3]. To use the correlations in this manner as- 
sumes implicitly that they reflect mineral/melt fractiona- 
tions (e.g., partial melting, fractional crystallization, 
accumulation) among a cogenetic suite of igneous 
rocks. However, similar correlations can qjpear from 
mixing of chemical components, and need not imply 
that the components were cogenetic nor that they had 
an igneous origin. 

Several aspects of the elemental and isotopic chem- 
istries of the Martian meteorites can be explained by 
mixing of geochemical components (mantle or crustal), 
e.g. [5,6]. Thus, it is worth asking whether mixing rela- 
tions might be seen in bulk chemical data, and (if so) 
whether the mixing components can be used to infer 
average planet or mantle compositions. 

Data: Data are from compilations [2, 7]. Where many 
analyses were available, a single average value was 
used. Where there were few or discrepant analyses, 
judgment favored large samples and recent esults. 
Considerable scatter is expected for elements (S, Se) 
that are concentrated in minor or trace phases. 

Results: For the most part, new analyses confirm 
the correlations noted >20 years ago. Treiman et al. [2] 
found that nominally incompatible elements could be 
divided into two groups that covaried coherently (Fig- 
ure 1). Group 1 elements (highly incompatible) include: 
La, Ag, As, Au, Ba, Br, Cl, Cs, K, Li, Rb, Sb, Ta, Th, U, 
and W. Group 2 elements (moderately hcompatible) 
include: Ti, Al, Cu, Eu, Ga, Hf, the HREE, In, Na, P, S, 
Sc, Se, Tl, and V. Elements in one group did not covary 
cleanly with those of the other group (Figure 1). Most 
of these correlations are as expected for basaltic sys- 
tems; e.g., highly incompatible elements do not enter 
solid phases, and so become equally concentrated in 
progressive melt fractions. Correlations of Ti with Na, 
Al, and Ga are also expected for Martian basalts, where 


plagioclase is not a liquidus phase. A few elements 
were found to behave compatibly (Ni, Cr), and a few 
‘indifferent’ elements had constant abundances across 
the whole sample suite (e.g., Co, Re). 


Figure 1 . Correlations of P with La and Ti in bulk chemical 
analyses of the Martian meteorites. Data from [2,7]. 




Unexpected Results . The behaviors of a few ele- 
ments are unexpected or bizarre. Abundances of phos- 
phorus closely track those of other moderately com- 
patible elements (Fig. lb), whereas in most cases P 
tracks highly incompatible elements like La (Fig. la). 
Second, abundances of the chalcogenides S and Se 
also track those of Ti and other moderately compatible 
elements (Fig. 2). In terrestrial basalts and mantle rocks, 
S and Se are not correlated with any lithophile elements 
but track other chalcogenides (Re, Ag) so to suggest a 
separate sulfide phase in the source [2]. 

Interpretations: The element-element correlations 
listed above have been interpreted as products of ig- 
neous fractionation, but might also represent mixing 
lines. 

The correlation of P with other moderately incom- 
patible elements is a problem - it requires that P the 
average bulk p D so]idsAMsaH be > 0.1, significantly higher 
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Mars its ownself: A.H. Treiman 



than in typical terrestrial petrogenesis. This can per- 
haps be justified by invoking high P-T fractionation 
with olivine fD ol/basa | t -0.1) and/or garnet £ D gar/basalt - 
0.3) [8,9]. Even higher bulk D values may be required by 
the fractionation of HREE from LREE [10]. 

The correlation of S and Se with moderately incom- 
patible lithophile elements is difficult to explain by ig- 
neous processes. Neither sulfide nor sulfate (nor se- 
lenide nor selenate) has is significantly soluble in com- 
mon rock-forming minerals. Igneous sulfide is common 
but minor in the Martian meteorites, and there is no 
reason to expect that it would fractionate coherently 
with silicate minerals (e.g., pyroxene). The Martian me- 
teorites contain no igneous sulfate minerals. They do 
contain rare secondary sulfate minerals but there is no 
obvious reason that the abundance of secondary sul- 
fates should follow those of Ti, HREE, Al, etc. 

Mixing as an Explanation (?): In a formal sense, the 
co-linearities seen in bulk chemical data can be ex- 
plained by mixing of at least three components: a ~ a 
reservoir rich in highly incompatible (Group 1) ele- 
ments; b - a reservoir rich in moderately incompatible 
(Group 2) elements; and c - a reservoir strongly de- 
pleted in all incompatible elements and so presumably 
rich in ferromagnesian and compatible elements. The 
meteorite richest in component a would be Nakhla, 
which has high La and low Ti and P. This assignment is 
stronger for Nakhla ’s parent magma, which is inferred 
to contain Ti=0.46%, P=0.26%, La~6ppm [10,11]. The 


meteorite richest in component b would be QUE94201, 
which in has the highest abundances of P and Ti, but 
low La etc. The meteorites seemingly richest in compo- 
nent c are the olivine-phyric basalts like DaG and SaU, 
which have low abundances of elements in both 
Groups 1 and 2 (e.g., La and Ti). 

These bulk-chemical components have some corre- 
spondence with components derived from isotopic 
systematics [5,6,12]. It may be possible to “fie 
the isotopically defined components into estimates of 
rock compositions, and thus better constrain their ori- 
gins. 

Estimates of the bulk composition of Mars are 
thrown in some doubt if the element-element correla- 
tions are mixing lines. Then, one has to ask if the mixing 
components represent only products of an igneous 
event (e.g., magma ocean differentiation). If so, then the 
element-element correlations may reflect planetary 
compositions. If not, or if the components have been 
affected by other processes (e.g. aqueous alteration of 
crust with introduction of sulfur), then all bets are off. 
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Introduction: Los Angeles is a non-brecciated, highly 
shocked basaltic shergottite. The shock effects in Los Ange- 
les can be classified into two basic types: (1) mechanical 
deformations and transformations that take place essentially 
in the solid state; and (2) localized regions of shock melting 
forming glassy to microcrystalline enclaves. A brief descrip- 
tion of the former type, based largely on optical observations, 
is given by [1]. They describe the complete transformation of 
plagioclase to maskelynite, mosaicism and planar fracturing 
of clinopyroxenes, pyrrhotite veins within clinopyroxene, 
polycrystallinity of some pyrrhotite grains, and a 350 pro- 
long fault with -15-pm displacement. These results are in 
agreement with this study. 

This work presents a detailed description of the shock- 
related microtexture and mineralogy of the localized melt 
pockets which have formed by in situ melting of local min- 
eral phases. The melt pockets range from 3 mm x 3 mm - 
0.07 mm x 1.25 mm, in direct contact with the host matrix. 
Two polished thin sections of stone 1 (1.3 cm x 1 cm, 
2.1 cm x l.2cm) have been investigated to determine the 
mineralogy and microtextures of the shock-induced melt 
pockets using a JEOL 6400 digital scanning electron micro- 
scope equipped with a Link Analytical eXL energy disper- 
sive spectrometer (EDS) fitted with a Si (Li) LZ-4 Pentafet 
detector, and a FEG-SEM equipped with the In-Lens Ther- 
mal FEG which produces high probe current sufficient for 
EBSP, WDS, and EDS. These instruments are capable of 
characterizing compositional variations and microtextures of 
the melt pockets that are unresolvable using optical tech- 
niques. 

Melt pockets — an overview: Los Angeles stone 1 
contains several unusual localized regions of shock melting 
which are texturally very different from typical melt veins 
observed in other SNC meteorites (e.g., Zagami) and ordi- 
nary chondrites. Melt pockets are contained within the two 
dimensions of the thin sections. Optically the melt pockets 
appear as irregularly-shaped dark brown to black coloured 
areas which are often observed to be spatially related to 
large, open fractures infilled with calcite. Three texturally 
distinguishable types have been observed: (1) Small 

(0.5 mm x 0.70 mm; 0.07 mm x 1.25 mm) vesiculated, clast- 
rich melt pockets with abundant crystallites indicative of 
quenching. Entrained clasts have subangular morphologies 
and consist of isotropic whitlockite, isotropic plagioclase 
significantly enriched in K 2 0 (Or 19 - 2 g) compared to maske- 
lynite in the adjacent host rock (An^Or^n), and Si0 2 clasts 
(95.02-98.35 wt% Si0 2 with variable amounts of Ai 2 0 3 ). 
Clast size ranges from 200 pm to 0.5 pm; avg. = 39 pm (n = 
100). The contact with adjacent host rock minerals varies 
from smooth boundaries to intrusive stringers and veins of 
shock-induced melt material One such melt pocket contains 
clasts which display interior quench textures, possibly inher- 
ited from a previous impact event. (2) Small (0.03 x 
0.04 mm) vesiculated, clast-poor (relatively) melt pockets 
with abundant crystallites of hollow Si0 2 rods (with minor 


amounts of A1 2 0 3 ). (3) With increasing volume of melt 
(3.5 mm x 3.0 mm) textures are dominantly vesicular, clast- 
poor (relatively) with flow textures and abundant 
spheres/blebs of immiscible Fe-sulphides. These range in 
size from 20 pm to >0. 1 pm and comprise mono- and poly- 
crystalline pyrrhotite. Most vesicles are spherical, indicative 
of quenching at low confining pressure, however, some are 
elongated with their c-axis parallel to the direction of flow. 
Opaques within these shock-produced glasses are perfectly 
spherical. It is interpreted that the vesicles were elongated by 
fluidal flow in the molten host, followed by quenching to 
produce glass, and spherical globules of sulphides. Vesicles 
range from 0.5 mm to >1 pm apparent diameter. 

Quench Textures: The crystals formed during quench- 
ing have a variety of distinctive shapes (skeletal, dendritic, 
feathery) and comprise a variety of minerals (olivine, 
whitlockite, ulvospinel, plagioclase). Quench crystals are 
invariably microscopic (microlites, crystallites). The 
groundmass consists of glassy material rich in Fe, Mg and Ca 
and is relatively pristine showing no signs of devitrification 
or alteration to aggregates of new phases (e.g., clays, idding- 
site). 

Severely fractured, melt intruded clinopyroxenes adja- 
cent to the melt pockets allow for preferential nucleation of 
dendritic crystals that extend from the contact into glasses. 
These consist of intimate intergrowths of pigeonite and or- 
thopyroxene with dendritic textures. In addition, clasts 
within melt pockets serve as nucleation sites for dendritic 
growth. These textures are consistent with rapid crystalliza- 
tion of the melt. 

Discussion: Remobilization and injection of opaques 

and melt from the melt region into pre-existing fractures in 
adjacent pyroxene in the host rock could have occurred dur- 
ing compression and subsequent stress relaxation. Although 
the contact between the melt pockets and adjacent host rocks 
is observed to be smooth to arcuate, it is also intrusive in 
areas, it is thought that these could have formed by pressure 
generated by a volume increase associated with the melting 
event, and that the melt pockets have formed in situ . With 
increasing volume the melt pockets comprise texturally het- 
erogeneous (flow-structured and meandering) schlieren-rich 
glasses. These show compositional variability related to the 
melting of variable proportions of the major rock-forming 
minerals (plagioclase and pyroxenes) of the host meteorite 
and incomplete mixing of melt. Abundant quench textures 
attest to rapid cooling of the shock-induced melt. The high 
abundance of relatively cold clasts in the smaller melt pock- 
ets probably contributed to the quenching of the melt to gen- 
erate abundant crystallites and an absence of flow textures. 

References: [1] Rubin, A.E. et al. (2000) Los An- 
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